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Ihis  report  describes  a  series  of  FCRTBAH-IV  subroutines  used  to 
prepare  graphical  displays  of  three-dimensional  data. 

Considerable  work  has  been  expended  to  make  these  programs  as 
flexible  as  possible  °  Many  are  based  on  routines  obtained  from  the 
University  of  Michigan  Department  of  Geography  or  the  Harvard  University 
Laboratory  for  Computer  Graphics,  and  are  believed  to  be  among  the  best 
routines  available  <, 

It  is  hoped  to  produce  a  revised  version  of  this  report  during  the 
summer  of  I968  which  will  contain  additional  programs  now  in  the  developmental 
stage o  Comments  from  users  regarding  corrections  or  modifications  to  the 
routines  are  welcome  <> 
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PARTI 

hjtsoixetion 


INTRODUCTION 

The  graphical  display  of  three-dimensional  data  is  important  in  many 
acientific  and  engineering  studies.  Such  data  form  surfaces  which  may  be 
of  many  types  ranging  from  real  surfaces  such  as  the  surface  of  the  earth 
to  various  mathematical  and  functional  relationships „ 

This  manual  describes  an  integrated  series  of  routines  which  can  be 
combined  to  produce  various  types  of  graphical,  displays  of  thrse-dimensional 
surfaces  using  either  the  printer  or  the  CALCOMP  plotter.  Use  of  the 
printer  has  several  advantages  over  the  CALCOMP  plotter,  particularly  in 
the  early  stages  of  a  project,  since  the  output  is  produced  as  part  of  the 
regular  job  output,  the  printer  is  cheap  and  fast  (allowing  many  more  plots 
to  be  produced),  and  its  use  does  not  require  exts>a  tape  mounts.  Due  to 
the  limited  resolution  of  the  printer  which  is  controlled  by  the  line  and 
character  spacing,  printer  displays  may  be  distorted,  although  under  favors- 
able  conditions  distortion  can  be  eliminated.  Since  the  CALCOMP  routines 
are  designed  to  be  easily  compatible  with  the  printer  routines,  the  plotter 
can  be  used  for  final  drafting  of  selected  displays. 

PROGRAM  LANGUAGE 

These  routines  are  written  in  FORTRAN-IV ,  Most  of  these  routines  can 
be  run  on  either  the  IBM  TO9U  or  the  CDC  6500  with  only  minor  changes  to 
accommodate  the  different  word  lengths  on  the  two  machines. 

PROGRAM  AVAILABILITY 

FORTRAN- IV  source  decks  are  available  for  copying  in  the  Airphoto 
Interpretation  and  Photogrammo  Room  B-2U,  Civil  Engineering 

Building  and  in  Room  G-150^  M a.'Hi •  Science**  BuiS^io^. 

METHODS  OF  SURFACE  REPRESENTATION 

The  three  most  common  methods  of  graphically  displaying  three-dimensional 
data  are: 

(1)  perspective  or  isometric  views  of  the  surface 

(2)  «hree-dimensional  histograms 

(3)  contour  maps, 

Routines  presently  under  development  will  provide  the  capability  of  pre- 
paring perspective  and  isometric  views,  and  three-dimensional  histograms. 
It  is  hoped  that  these  new  routines  will  be  available  this  summer. 


This  first  edition  of  the  manual  describes  a  series  of  main  programs 
and  associated  subroutines  to  perform  the  following  functions  ■» 

a)  Production  of  contour  maps  utilizing  polynomial 
equations . 

b)  Interpolation  of  irregularly  spaced  to  gridded  data. 

c)  Production  of  contour  maps  of  the  printer* 

d)  Production  of  contour  maps  on  the  Calcomp  plotter, 

BENEFITS  AMD  LIMITATIONS  OF  MftCHIffE-COMFUTED  CCMTOUB  MAPS 

Contouring  can  be  extremely  literal,  extremely  interpretative,  or  a 
compromise  between  the  two-  The  routines  deve       re  give  literal 
contours.  As  such,  the  maps  produced  are  not  "unbiased,"  as  is  often 
claimed,  but  are  in  fact  strongly  and  consistently  biased  toward  literal 
interpretation  of  the  data. 

As  a  consequence,  computer«generated  contour  maps  developed  by  these 
routines  may  not  be  equally  suitable  for  all  applications.  They  are  ex- 
tremely useful  as  "quick-look*0  maps  to  give  the  user  a  general  idea  of 
his  data,  or  to  check  the  general  validity  of  the  data. 

Furthermore,  since  all  maps  are  made  consistently,  they  can  also 
be  used  for  comparison  purposes .  Contouring  several  data  sets  by  humans 
inevitably  involves  some  learning  processes  which  introduce  an  uncertain 
and  subtle  bias  into  the  results  rendering  comparison  of  maps  difficult . 
In  contrast,  the  computer  will  not  "learn w  from  one  problem  to  the  other, 
but  will  handle  each  data  set  uniformly. 

PROCEDURE 

Computer  contouring  of  data  involves  two  steps 

1)  determination  of  regularly  spaced  (gridded) 
values  of  the  dependent  variable  (elevation) 

2)  linear  interpolation  to  locate  contour  lines 
within  the  grid  squares. 

(Obviously,  if  data  is  already  in  grid  form, 
step  1  may  be  omitted.) 


Determination  of  Gridded  Values 

l)  By  Regression  Analysis 

The  determination  of  regularly  spaced  values  can  be  accomplished  by 
first  developing  an  analytic  equation  to  describe  the  data,  and  then 
evaluating  this  equation  at  a  series  of  grid  coordinates.  Least  squares 
regression  procedures  based  on  the  polynomial  expansion  are  commonly  used, 
although  double  Hurler  series  expansion  procedures  have  also  been 


j  - 

developed  (2).  Statistical  procedures  allow  the  research  to  determine  the 
ability  of  the  various  surfaces  to  conform  to  his  data  (l,  3).  These  techniques, 
termed  Trend  Surface  Analysis,  have  long  been  used  by  geologists  and 
geographers  (3,  k,   fi,  7).  They  have  found  residual  maps,  contour  maps  show- 
ing the  lack  of  fit  to  the  regression  equations,  useful  in  locating  anoiaalous 
areas  (6). 

Program  POLYFIT  and  its  associated  subroutines  will  fit  polynomials  up 
to  the  fifth  power  and  produce  trend  surface  maps.  Program  RSSMA?  will 
produce  residual  maps. 

2)  By  Weighted  Moving-Average  Procedures 

Weighted  moving-averages  or  the  fitting  of  local  polynomials  are  two 
other  techniques  for  interpolating  irregularly  spaced  da  ta  to  a  uniform 
grid.  Weighted  moving  averages  are  generally  computationally  easier,  and 
therefore  faster,  than  the  fitting  of  local  polynomials. 

Subroutine  GRID,  based  on  en  algorithm  developed  by  Professor  W.  R. 
Tobler  at  the  University  of  Michigan,  uses  a  weighted  moving  average  to 
determine  grid  values.  As  in  all  interpolation  procedures  soma  smoothing  of 
the  data  results  from  its  use.  A  limited  amount  of  testing  suggests  that 
the  smoothing  is  normally  less  than  2$  of  the  data  range,  rising  to  nearly 
10$  of  the  range  with  extremely  erratic  data. 

Linear  Interpolation  Within  Grid  Squares 

Subroutines  F^TCOKl  and  P&CC0N2  perform  this  step  on  the  printer; 
program  COKTUR,  performs  a  similar  operation  on  the  CALCOMP  plotter „ 


II  and  Pg|rC0N2  differ  only  in  the  type  of  output  producedo 
P$TC0N1  produces  maps  with  various  characters  representing  different 
contour  intervals.  Options  allow  for  the  printing  of  contour  lines  only, 
alternate  contour  bands,  or  all  print  positions  within  the  map.  P$|TCQN2 
produces  maps  with  contour  intervals  delineated  by  combinations  of 
overprinted  characters  to  give  a  graded  density  from  light  (low  areas)  to 
dark  (high  areas).  Since  the  overprint  carriage  control  is  only  available 
on  the  CDC  6500,  this  program  is  restricted  to  that  machine. 

The  CALCOMP  plotting  routines  are  much  mere  effi@ient  on  the  IBM  709U; 
plotting  speeds  are  at  least  twice  as  great ?  and  plot  tapes  are  written  at 
a  higher  density.  Under  these  circumstances  program  CONTUR  has  not  been 
converted  to  run  on  the  CDC  6500,  although  no  unusual  difficulties  are 
anticipated  for  such  a  conversion. 

SUMMARY 

A  series  of  integrated  routines  are  available  to  perform  machine  con- 
touring of  three-dimensional  data  on  either  the  printer  or  the  CALCOMP 
plotter. 
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PART  II 

PROCEDURES  FOR  MAP 
CONTOURING  UTILISING 
REGRESSION  EQUATIONS 


NOTE  •  Program  POLYFIT  (and  associated  subroutines)  should,  be  used 

if  polynomial  regression  surfaces  are  to  be  fitted  to  regularly 
or  irregularly  spaced  data. 

Program  RESMAP  (and  associated  subroutines)  should  be  used  in 
conjunction  with  program  POLYFIT  to  analyze  residual  and 
original  data. 


CONTENTS  OF  THIS  SECTION 

1)  Program  POLYFIT  -  including  subroutines 
READ1,  DIVER,  DJVCK,  POLY,  EVALU, 
and  PLOT 


2)  Program  RESMAP 


PURPOSE  — 

TO  PERFORM  TREND  SURFACE  ANALYSIS  USING  POWER  SERIES  EXPANSION,, 
POLYNOMIALS  UP  TO  FIFTH  DEGREE  CAN  BE  REQUESTED.  THIS  PROGRAM  HAS 
EVOLVED  FROM  AN  EARLIER  PROGRAM  OBTAINED  FROM  PROFESSOR  R.B.JOHNSON* 
FORMERLY  HEAD  OF  THE  DEPARTMENT  OF  GEOSC IENCES. PURDUE  UNIVERSITY. 

PURDUE  UNIVERSITY  VERSION  4.2       A  KEITH  TURNER* 

CIVIL  ENGINEERING, 
FEBRUARY*  1968. 


THIS  VERSION  CAN  FIT  A  SURFACE  TO  1000  DATA  POINTS.  UP  TO  FOUR 
DIFFERENT  VARIABLES  CAN  3E  ANALYZE:-  AT  EACH  POINT.  DATA  CARDS  (ONE 
PER  POINT)  SHOULD  INCLUDE  IDENTIFICATION, X-CO0RDINATEtY-COORDINATE» 
ANH  VARIABLES  1-4  ( Z 1 ?Z2 »Z3»Z4 ) .  AN  UNLIMITED  MUMPER  OF  DATA  SETS, 
EACH  CONTAINING  UP  TO  1000  POINTS*  CAN  3E  PROCESSED. 

THE  PROGRAM  WILL  (1)  FIT  THOSE  POLYNOMIAL  SURFACES  REQUESTED* 
(2)  LIST  THE  MATRICES  USED  TO  DETERMINE  THE    COEF- ICI ENTS. < 3)  LIST 

COEFFICIENTS*  14)  ESTIMATE  THE  ERROR  J.N  EACH  COEFF ICIENT*  ( 5 )  PER- 
FORM ANALYSIS  OF  VARIANCE*AND(6)  LIST  THE  ORIGINAL  DATA, COMPUTED 
VALUES.AND  RESIDUALS  FOR  ALL  POINTS, 

SUBROUTINES  EVALU  AND  PLOT  GENERATE  AND  MACHINE-PUNCH  6RIDDED 
VALUES  FOR  SELECTED  TRENC  SURFACES  AND  PRODUCE  PRINTER-CONTOUR  MAPS 
OF  T;       IRFACES.A  VARIETY  OF  CONTOUR  MAPS  CAN  BE  SPECIFI ED.MULT- 
PLE  'APS  OF  ANY  SURFACE  MAY  BE  PROuM'       HE  GRIDDED  VALUES  MAY  PRO- 
DUCE CALCOMP  CONTOUR  MAPS  USING  PROGRAM  CONTUR. 

ALL  ORIGINAL  AND  COMPUTED  VALUES*       LSIDUALS*  FOR  THE  DATA  CAN 
BE  STORED  ON  A  USER-DESIGNATED  TAPE  (TAPES)  FOR  LATER  ANALYSIS  • 


TAPE  REQUIREMENTS  — 

TAPE1  IS  A  SCRATCH  TAPE  USED  TO  STORr.  AND  REGENERATE  DATA  DEST- 
ROYED DURING  THE  OPERATION  OF  SUBROUTINES  EVALU  AND  PLOT.  NO  REQUEST 
CARD  IS  NECESSARY  SIMCE  A  DISK  FILE  WILL  AUTOMATICALLY  BE  GENERATED- 

TAPES  IS  USED  TO  STORE  THE  ORIGINAL  INPUT  VALUES .COMPUTED  VALUES 
AND.  RESIDUALS  FOR  ALL  POINTS. THF  TAPE  SHOULD  nF  '•/RITTEN  IN  BINARY 
AT  556  BPI.  THE  REQUEST  CARD  FORM  Is  — 

REQUEST  (TAPE8,HI)  PLEASE  MOUNT  USE RTAPENAME .FILE  XXX 

THE  FIRST  TAPE  RECORD  WILL  CONTAIN  THE  NUMBER  OF  DATA  POINTS  (N)  BE- 
ING ANALYZED»THERE  THEN  FOLLOW  N  RECORDS  CONTAINING  THE  DATA  ID.X.Y* 
ZtlST  DEGREE  ESTIMATEtlST  DEGREE  RESIDUAL* »5TH  DEGREE  RESIDUAL. 

THIS  SEQUENCE  WILL  Br  REPEATED  FOR  EACH  VARIABLE  (ONE  UP  TO  FOUR 
TIMES)  FOR  EACH  DATA  SET  OF  1000  POINTS  OR  LESS.  AN  END-OF-FILE  IS 
WRITTEN  AFTER  THE  LAST  DATA  SET. 

IF  IT  IS  DESIRED  TO  IMMEDIATELY  ANALYZE  THE  RESIDUALS*  PROGRAM 
RESMAP  MAY  BF  INCLUDED  UNDER  THF  SAME  JOB  CARD.  IN  THIS  CASE  NO  TAPE 
REQUEST  CARD  IS  REQUIRED. 


C    RQUTIN""*!  ~r~UT~r'"  — 

C     A)  MAIM  PROGRAM   (POLYFIT) 

C 

C     ^)  SUbROUTI 

C         1?  READ1  -  READS  THE  DA 

2}  IMVER  -  DETERMINES  INVERi       'X. 

C         ?■)     INVCK  -  PRINTS  OUT  >■        X  MATRIX* 

C         ')  POLY  -  DETERMINES  ESTIMATED  VALUES  AND  RESIDUALS. 

b>  EVALU  -  EVALUATES  UD    POINTS. 

C         6*  PLOT  -•  PRODUCES  PRINTER  CONTOUR  .•.APS. 

C         7)  STATS  -  CALCULATES  STATISTICAL  MEASURES. 

C    DESCRIP1        CONTROL  CARDS  —  ALL  CARDS  FOR  EACH  DATA  SET; 

C   1}  LARLE  CARD  - 

'FT  (COLS  1-6  5  COD        INDICATING        OF  DATA 


- 


5  MASTER  CONTROL  CARD  - 
POLYFT  (COLS  1-6  1 

C         LIMI T ( 1-5 )(COLSl 1-15)  IF  LIMIT (K)=l .SURFACE  OF  DEGREE  rf.  FITTED. 
ISURFU-5HCOLS16-20)  (F  ISURF(K)=1 .GRID  VALUES  FOR  SURFACE  K 

ARE  CALCULATED  (THESE  WILL  BE  PUNCHED  IF  SPECIFIED  ON 
CARD  9).  HOWE        WAP  WILL  3E  PRINTED. 

C  I  SURF  {<.)  =2  .GRID  VALUES  WILL  BE  DE7ER- 

C  MINED  AND  MAP(S)  OF  SUFACE  K  WILL  BE  PRINTED. 

C         INC       25)  IF  INCK»1»X-PRIME-X  MATRICES  PRINTED. 

C         NLI        iL30)NUMBER  OF  GENERAL  TITLE  CARDS  TO  FOLLOW  tMXIKUM"5 

C         NZ.U.-4)  (COL  537-40)  IF  HZ  CO  =1 » VARIABLE  K  WILL  BE  ANALYZED* 

C         NRES  (COL  50)  IF  SluUAL  VALUES  LISTED* 

C         TOUT 1 (COL  60)  IF  [  TEN  ON  TAPE8  . 

C 

C   3)  GENERAL  TITLE  CARD(S)  - 

C  CARDStNUMSER  EQUALS  NLINES).EACH  CARD  CONTAINS  ANY  ALPHA- 

C  NUMERIC  TITLE  IN  COLS  1-72 

C 

VARIABLE  IDENTIFICATION  CARD  - 

C         FOUR  SETS  OF  ALPHANUMERIC  TITLES  (COLS  1-20.21-40.41-60*61-80) 

C  WHICH  WILL  IDENTIFY  VARIABLES  Z1-Z4. 

C  NOTE-IF  NO  FURTHER  IDENTIFICATION  REQUIRED. INSERT  bLANK  CARD 

C   5)  FORMAT  CARD  -  ,  , 

C  AT  OF  DATA  (COLS  7-78)* SHOULD  oE  IN  FORM  ( A 10 » &FLOATI NG 

C  POINT  FIELDS{X»v3.Zl9Z29Z39Z4)  ) 

C 

C   b)  DAT/  CARDS  - 

C         (  ME  CARD  FOR  EACH  POi 

C 

C   7)  BLANK  CARD  - 

C  ID  TO  DISTINGUISH  END  OF  DATA  SET. 

C    ~ — ■ ■ — ' — — — 


IF  THI  GRlDDING  OR  GRIDOING  AND  "f\Prj:\r    3*>Tir.'!S  ARF  ""'E^UESTCr 

(IF  IbURF(K)=l  OR  2)  •  THE  FOLLOWING  CARDS  MUST  BE  3EP*"ATT"»  e"3!!  f-'C 

REOUFST  — 

6)  ARFA  DEFINITION  CARD  - 

MAXIMUM  AND  MINIMUM  VALUES  OF  X  AND  Y  COORDINATES  (XMAX»XMIN» 
YMAX»YMIN)  WITH  DECIMALS  PUNCHEDCCOLS  l-10»ll-20»21-30»31-40) 

9)  GRID  SPECIFICATION  CARD  - 

ROWSICOLS  1-3  J  INTEGER  NUMBER  OF  ROWS  IN  GRID i MAXIMUM=100 ) e 
COLS(COLS  4-6)  INTEGER  NUMBER  OF  COLUMNS  IN  GRID (MAXIMUM=1C0 ) • 
NMAP(COLS  7-9)  NUMBER  OF  DIFFERENT  MAPS  OF  THIS  SURFACE  TO  BE 

PRINTED, 
IPUNlCOL  10)  IF  IPUN=1»GRID  VALUES  PUNCHED  (EACH  SURFACE  IDENT- 
IFIED BY  PUNCHED  TITLE  CARDS) 

NOTE-A)  FOR  1  INCH  GRID  COLS-MAP  WIDTH   (INCHES)  +  1 

ROWS»MAP  LENGTH  (INCHES)  +  1 
B)  IF  COLS«GT»42, PRODUCTION  OF  PRINTER  MAPS  IS  IMPOSSIBLE 
DUE  TO  LIMITED  WIDTH  OF  PAPER  (   HOWEVER 
CALCOMP  PLOTS  ALLOW  UP  TO  lOQXlOO  ARRAYS) • 


REPEAT  THE  FOLLOWING  CARDIS)  FOR  EACH  MAP  REQUESTED  BY  NMAP  - 

10)  PLOT  CONTROL  CARD  - 

CON  (COLS  1-2)  NUMBER  OF  CONTOUR  I  NTERVALS»MAxXMlJM  =  19 

IF  CON  =  0*  CON  SET  TO  6* 
TOUR (COL  3)  IF  TQUR-OsMAX*  ANu  MlN.  DATA  ELEVATIONS  USED  AS 

CONTOUR  LIMITS, 

IF  TOUR=l»DESlRED  TOP  AND  BOTTOM  CONTOURS  READ  IN 

FRO^  SPECIAL  CARD  (SFF  BELOW) 

IF  TOUR=2»VARIABLE  CONTOURS  "FAD  IN  (LOW  TO  HIGH). 

IN  THIS  CASE  SET  CON  EQUAL  TO  THE  NUMBER  OF  VALUES 

TO  BE  READ  IN. 

IF  T0UR*3»C0NT0UR  INTERVAL  BECOMES  1/2  STANDARD 

DEVIATION  INCREMENTS* I  MAX=+3  SIGMA,  MIN=-3  SIGMA  ) 
!  INFS(COL4)  IF  LlNPS*0t  CONTOUR  LIMES  ARE  PRINTED* 

IF  LlNES=ltALTERNATE(EVEN-VALUED>  BANDS  ARE  PRINTED 

IF  LlNfrS  =  2»  COMPLETE  BANDS  ARE  PRINTED,. 
IMCHESCCOLS5-9)  WIDTH  OF  MAP  IN  I NCHES»MAXlMUM=12*7 

11)  SPECIAL  CARD(S)   (  REQUIRED  ONLY  IF  TOUR  =  1  OR  2  ) 

IF  TOUR  =  l»  THIS  CARC  SPECIFIES  MINIMUM  AND  MAXIMUM  ELEVATION 
VALUES*  ZMIN(COLS  1-10 ) »ZMAX(COLS  11-20)  DECIMALS  PUNCHEDe 
IF  TOUR=2,  THIS  CARD  SPECIFIFS  A  FORMAT  (INCON)  USED  IN  READING 
A  SERIES  OF  IRREGULARLY  SPACED  CONTOUR  VALUES*  THE  VALUES  ARE 
ARRANGED  FROM  LOWEST  TC  HIGHEST  ON  SUBSEQUENT  CARDS*  SET  CON= 
NUMBER  OF  LINES* 


C  PROGRAM  PCLYFIT  MAIN  P;*06RA^ 

PROGRAM  POLYFIT( INPUT»CUTPUT»PUNCH»TAPE1 ,TAPE8»  APE5  = I NPUT,TAPE6= C 
1UTPUT) 
C 
C 

DIMENSION  XX(21»215»  P(21),  A(21»5),  XINV( 21 ,21 » 5 ) ,  PAR(21),  ZR(10 
1),  SSR(5),  BB<21»4),  SSQZ(4)»  X(lOOO),  Y(lOQn),  Z(1000»4),  IDClOOO 
2),  LIMITC5)*  NZ(4),  FRR(21),  SSF<5),  DFR£E(5) 

DIMENSION  ISURFJ5),  TITL<895),  ORDIN(4,2) 

DIMENSION  DUM1C3000) 

COMMON  X  ,  Y  ,  Z  ,  I D  »DUM 1  ,X  I NV * A ,KSOR 

DATA  ND0G/6HP0LYFT/ 

XD»O«0 

YD»0,u 

ZD=0.u 
C 

C      SEARCH  THE  INPUT  TAPE  FOR  THE  PARAMETER   CARD  WITH  THE  CODE  POLyFT 
C 
5      READ  (5*370)  LA3LE 

IF  (E0F.5)  310,10 
10     IF  (LAPLE-NDOG)  5,15,5 

15     READ  (5,375)  LAbLE,LIM I T» I  SURF, INCK,NLINES»NZ »NRES, IOUTl 
C 

C    WRITE  OUT  HEADING 
C 

WRITE  (6,455) 

DO  20  I=1,NLINES 

RFAD  (5,32))  ( TI TL (  I  I  ,  I  )  *  I  1  =  1 »8  ) 
20     WRITE  (6,460)  (TITL { I  I  ,  I  )  .  I I«l»8 ) 

WRITE  (6,380) 

READ  (5,465)  ( ORDIN ( I , 1 ) ,ORDl N ( I ,2 ) , I =1 ,4 ) 

SS0Z(D  =  o.0 

SSQZ(2)*0,0 

SSQZ(3)=U„0 

SSQZ(4)=u„0 

DO  30  1=1,5 
J=6-T 

IF  (i.IMITC  J)  )  25,30,25 
25        N=( (J+l)*< J+2) >/2 
JUMP*  I 
GO  TO  35 
C 
C 
30     CONTINUE 

GO  TO  305 
C 
C 

C       READ  THE  DATA  POINTS 
C 

35     CALL  RFAD1  (NCARDS) 
C 


C      COMPUTE  THE  COEFFICIENT   MATRIX 
C 

PAR(1  >  =  '•.. U 
DO  45  1^1  »N 

DO  4^  J=I »N 
40        XX(I»J)=0„0 

BB(I,1>=0.0 

BBC  I »2) =0.0 

BB(I»3)=0*0 
45     BB(I,4)=o.O 

DO  95  K=1»NCARDS 

XD=X(K) 

YD=Y(K) 

PAR(2)=XD 

PAR{3)=YD 

GO  TO  (  50  *  50  »  50*50*7.0)  ,  JUMP 
5u        PAR(4)=XD*XD 

PAR{5)=XD*YD 

PAR(6)=YD*-YD 

GO  TO  (55? 55,5  5,70*70)*  JUMP 
55        PAR(7)=PAR{4)*XD 

PAR<8)=PAR(4)*YD 

PAR(9)=PAR(6)*XD 

PAR(l0)=PAR(6}*YD 

GO  TO  (60,60*70,70,70) ,  JUMP 
60        PAR{11)=PAR{7)*XD 

PAR(12)<=PAR(7)*YD 

PAR{13)=PAR(41»PAR(6> 

PAR«14)=PAR{10)*XD 

PAR(15)=PAR{ 10)*YD 

GO  TO  {65,70,70,70,70),  JUMP 
65        PAR«16)=PAR( 11)*XD 

PAR(17)=PAR(11)*YD 

PAR < 1 8  >  =p  AR ( 7 ) #PAR { 6 ) 

PAR«19)sPAR(10)*PAR(4) 

PAR(2U)=PAR(15)*XD 

PAR(21)=PAR{15)*YD 
70        DO  75  1=1, N 

DO  75  J=I,N 
75        XX(I»J)=XX( I»J)+PAR(I)*PAR(J) 
C 

c    compute.;  the  p  vfctors  for  each  z  to  re  fittfd. 
c 

DO  90  JIM=1„4 

IF  (NZ(JIM))  80,90,80 
80  ZD=Z(K»JIM) 

DO  35  I  =  1.N 
8  5  RBCIt JIM)sPB(I»JIM)+ZD*PAR( I ) 

SSQZ ( J I M ) =SSQZ ( J  i  M ) +ZD*ZD 
90        CONTINUE 
95     CONTINUE 


SO 


WRITF  (6*325) 
DO  100  I-1»N 

WRITE  (6.330)  I , (XX ( I » J ) » J-l.N ) 
100    CONTINUE 

WRITE  (6,335) 
DO  105  JIM=1,4 

IF  (NZ(JIM).EO.O)  r-0    TO  105 
DO  105  1*1, N 

WRIT!!  (6*340)  I,BB(  I  .JIM) 
105    CONTINUL 
NM1=N-1 
DO  110  1=1, NM1 

TP1=I+1 
DO  110  J=IP1,N 
110    XX( J,I)=XX( I.J) 
WRITE  (6.345) 
DO  115  1=1, N 

WRITE  (6,330)  I » ( XX ( I . J ) , J=1»N) 
115    CONTINUE 

DO  14u  K=l,5 

TP  (LIMIT(O)  120,140,120 
120       NTERMSS{ ( K+l )*<K+2 ) )/2 
KSOR=0 

CALL  INVER  { XX»XINV9MTERMS,K) 
IF  (KSOB)  125,130,125 
125       WRITE  (6,385)  < 

LIMITlKJsO 
130       CONTINUE 

IP  (INCK)  135,140,135 
135       CALL  INVCK  (XX,XlNV*lO 
140    CONTINUE 
REWIND  1 
C 

C    THE  JIM  DO  LOOP  RUNNING  TO  STATEMENT  300  CONTROLS  WHICH  Z  IS  BEING 
C      FITTED. 
C 

DO  300  JIM=1,4 

IF  (NZ(JIM))  145,300,145 
145       DO  150  I=1,NCARDS 

150       RFAD  (1)  ID( I),X(I) ,Y( i ) ,2( !,1),Z{ 1,2) .Z( I,3>,Z(I»4) 
SSZ=S5QZ« JIM) 
DO  155  1=1, N 
155       BiI)=b6(I,JIM) 
C 

C       COMPUTE  COEFFICIENTS 
C 

DO  1"?0  K  =  l,5 

I!*  (LIMIT(K))  16C*i?0,l&0 
160  NYERMS=(  (K+l)*(K+2U/2 


IL 


DO  165  I«1.NTERMS 

A(I,K)=0.0 
DO  165  JaltNTERMS 
165  A(I»K)=A< I,K)+XINV(  i »J»K)*b£ J  J 

170       CONTINUE 

WRITE  (6e455) 
DO  175  I=*1»NLINES 
175       WRITE  (6,460)  ( T F TL ( I  I • I) » I  I =1 ,3 J 

WRITE  (6*395)  ORDIN{ JIM, 1 > tORDIN< JIM.2) 
WRITE  (6,390)  NCARDS 
00  l;»5  K*l»5 

II'  (LIMIT(K))  180*185,180 
180  NTERMS*((K+l)*tK+2))/2 

WRITE  (6,400)  K,(A< l,K),I*l»NTERMS> 
185       CONTINUE 
C 

C      ANALYSIS  OF  VARIANCE 
C      ERROR  ESTIMATES  FOR  COEFFICIENTS 
C 

ZSUM=u,Q 

DO  190  K=1»NCARDS 
190       ZSUM=ZSUM+Z(K,JIM) 

ZBAR=Z SUM/FLOAT (NCARDS) 

SSZBAR=ZSUM*Z3AR 

SSZ=SSZ-SSZ3AR 

WRITE  (6,405) 

DO  210  K  =  l»5 

TF  (LIMIT(O)  195»210»195 
195  NTERM5"((K+l)*(K+2))/2 

SSR(K)=0.0 
DO  200  J»1,NTER?.*S 
200  SSR(K)=SSR(K)+A( J9K)*B(J) 

SSR ( fc ) =SSR ( K ) -SSZBAR 
SS£(K)*SSZ-SSR(K) 
DFREF(K)=NCARDS~NTFRMS-1 
Si GMAE«SGRT ( SSE { K ) /DFREE ( K ) ) 
DO  205  JaltNTERMS 
205  PRR(  J)«STGV»AE*SrtpT(XI\'V(  J»J»K)  ) 

WRITE  (6,4105  K»(ERR(J)  »J-1»NTER.^S) 
210       CONTINUE 
C 

C      ANALYSIS  OF  VARIANCE 
C 

WRITE  (6*455) 
DO  215  I=1,NLINES 
215       WRITE  (6*460)  ( TITL( 1 1 » I ) * T 1  =  1,8 ) 

WRITE  (6,470)  ORDIN ( JIM,1 ) ,0RUI N ( JIri,2 ) 

WRITE  (6»415) 

NDF=NCARDS-1 

WRlTi!    (6»350)    NDF.SSZtJAR 


00  225  K  =  l,5 

IF  CLIMITIK))  22C»225,22C 
22u  NTERMS=( (K+l)*(K+2> )/2 

XSE=SS£(K)/DFREE(K) 
XSR=SSR(K) /FLOAT (NTERMS) 
F=XSR/XSE 
NDF=DFREE(<) 

WRITE  (6,420)  K ,NTERMS»NDF»SSR (K ) »XSR ,SSE( K) ,XSF,F 
22  5       CONTINUE 

JACK=6-JUMP 
WRITE  (6»355)  SS2 
WRITE  (6»360> 
00  2'5  K=l,5 

II  (LIMIT(K))  23~»235*230 
230  RR=SSR(K)/SSZ 

PCENT=PR*l00o0 
R=SQRT(RR) 

WRITE  (6,365)  K,PCENT9RR»R 
235       CONTINUF 

WRITE  (6»425)  JACK 

KSWT»1 

00  255  K=l,5 

IF  (LIMIT(K)J  24C*255,240 
240  NTFRMS=( (K+l)#(K+2) )/2 

GO  TO  (245,250) ,  KSWT 
245  KSWT*2 

NTERMS=MTERMS-1 
XMSQ-SSR ( K ) /FLOAT ( NTERMS ) 
F=XMSQ/ (SSE( JACK )/DFREE( JACK) I 
WRITE  (6,430)  K, NTERMS, SSR(K),XMSQ,F 
NTE=NTERMS+1 
J=K 

GO  TO  255 
C 
C 
250  DF=NTERMS-NTE 

SSQ=SSR(K)-SSR(J) 
XMSQ=SSQ/DF 

F=XM5Q/ <  SSF  ( JACK ) /DFREE ( JACK ) 3 
NDF=OF 

WRITE  (6,430)  K,NDF,SSOtXMSQ»F 
J=K 

NTE=NTERMS 
255       CONTINUF 

IF  (NRES)  260,300,260 
260       WRITE  (6,455) 

DO  265  I=1,NLINES 
265       WRITE  (6,460)  ( Ti TL( I  I , 1 ) » 11  =  1 *8 ) 

WRITE  (6,435)  ORDI N ( JIM • 1 ) sORO IN ( JIM»2 ) 
WRITE  (6,440) 


nn  ?70   1=1 ,ic 
270  ZRU>30.0 

IF  (IOUT1.NE.1)  GO  TO  275 
WRITE  (8)  NCARDS 
275       DO  285  1=1, NCARDS 
XD=X(I) 
YD  =  Y(  I) 
ZD=Z( I»JIM) 

CALL  POLY  {ZR*LIMIT,XD,YD,ZD) 

WRITE  (6,445)  ID(I)«X(I)»Y(I)*Z(I»JlM)*(ZA(K)»K«lf 10) 
IF  (ICUT1.NE.1)  GO  TO  280 

WRITE  (8)  ID( I) »X(I)»Y( I)tZ(!»JIM) 9(ZR(K)»K*1»10) 
280  O'NTINUE 

285       CONTINUE 

DO  295  1  =  1,5 

IF  (ISURF(I).NE«1.AND.ISURF(I J.NE.2)  GO  TO  290 
CALL  FVALU  ( A  *  I , I  SURF, Tl TL *NLI NES .ORDIM » Jl M) 
290  CONTINUE 

REWIND  1 
295       CONTINUE 
300    CONTINUE 

GO  TO  5 
C 
C 
305    V/RITF  (6,450) 

GO  TO  5 
C 

c 

310    CONTINUE 

IF  (IOUTl»NE.l)  GO  TO  315 

FND  FILE  8 

REWIND  8 
315    CONTINUE 

STOP 
C 

c 
c 

320    FORMAT  (8A10) 

325    FORMAT  (1H1,13HTHE  XX  VATR IX » /// »4H  ROW) 

33u    FORMAT  (  1H0  *  I  3  ,2X,  1 1E1C  .3  »/6X,  10E10.3  ) 

335    FORMAT  (1H1,13HTHE  6B  MATRIX »/// »4H  ROW; 

340    FORMAT  ( 1HC *  I  3 ,2X ,E16.4> 

345    FORMAT  (1H1,21HTHE  REVISED  XX  MATRIX »/// ,4H  ROW) 

350    FORMAT  (1H0,12H   MEAN     1  * , I4,F17»8 ) 

355    FORMAT  ( ///30X ,47HC0RRECTED  TOTAL  SS=  TOTAL  SS  -  SS  DUE  TO  MEAN  =  * 

1F20*6) 

360    FORMAT  ( ////1X,81HP£RCENT  SS  EXPLAINED  BY  REGRESSlON=MULTI PLE  CORR 

1ELATION  COEFFICIENT  (  R-SQUARED  )»///*54H  DEGREE    PERCENT  VARIATI 

20N  FXPLAINFD   R-SOUARFD      R) 

365    FORMAT  < 1H0 s 14, 15X ,F10 .3 , 10X ,F7«5 , 3X *F7.5 ) 


370    FORMAT  <A6) 

375    FORMAT  ( A6 *4X»5 1 1 ,5  1 1 ,4X » I  1 » 3X ♦ I2«6X »4I 1 ,9X » II ,9X 9 1 1 5 

380    FORMAT  (  lHr. ,///  ,5X  »5?HTREND  SURFACE  ANALYSIS  USING  POWER  SFRIFS  EX 

1PANSI0N. »/,5X»48HP0LYN0MIALS  UP  TO  FIFTH  DEGREE  CAN  BE  REQUESTED./ 

2) 
385    FORMAT  (46H  THE  COEFFICIENT  MATRIX  FOR  THE  FIT  OF  DEGREE  ,IU70H  I 

IS  MACHINE  SINGULAR.   THE  PROGRAM  WILL  TRY  THE  OTHER  FITS  ASKED  FOR 

2.  ) 
390    FORMAT  (1H0»15X»39H   THE  NUMBER  OF  POINTS  IN  THIS  GROUP  IS»I5,//// 

1) 
395    FORMAT  •  lHC»2e>X»?8HTRe,ND  SURFACE  EQUATIONS  FOR  9    2A10» 

181HTHE  COEFFICIENTS*  WiTH  THE  CONSTANT  TERM  LISTED  FIRST.   THE  FOR 

2M  OF  THE  POLY  IS  //12ZH         Z*Al+A2X+A3Y+A4X2+A5XY+A6Y2+A7X3+A8X2Y+ 

3A9XY2+A10Y3+A11X4+A12X3Y+A13X2Y2+A14+XY3+AI5Y4-J-A16X5+A17X4Y+A18X3Y 

42+  ...  ETC.  //1X»7H  DEGREE) 
400    FORMAT  i 1H0 t I4»3X »7E16 .7//E24, 7?6E16 .7//E24.7 »6E16.7 ) 
405    FORMAT  (///51H  ESTIMATE  OF  ERROR  IN  THE  COEFFICIENTS  OF  THE  POLYS/ 

1/10H     DEGREE) 
410    FORMAT  ( 1H0 » 16  *  IX *7Fl£ .2 s//F24„? »6F16.2 »//F24.2 »6E16.2 ) 
415    FORMAT  { 56H0  ANALYSIS  OF  VARIANCES  FOR  ALL  REGRESSION 

?>S////lX»99H   SOURCE   D.F.      SS  DUE  REG.       MS  DUE  REG* 

2  SS  FRROR  MS  FPROP        F-RATIO  ) 

420    FORMAT  < lH0»7HDEGREE  ,  I  1 » I  3  »  1H , , 14 »E17.8 ,EI6.8,E19.8 ,F18.5»F14.3 ) 
425    FORMAT  <////34H  ANALYSIS  OF  VARIANCE  FOR  ,U»12H  DEGREE  PO 

1LY////75H  DEGREE  OF  SOURCE   D*F.    SUM  OF  SQUARES         MEAN  SOUA 

2RE         F-RATIO  ) 
430    FORMAT  ( 7X» ! 1» 12X» I2»F20.8 fE20.8»F17.5 J 
435    FORMAT  ( 1H0 » 18X »82HCOORDINATES»OR IGI NAL  AND  COMPUTED  VALUES»AND  RE 

1SIDUALS  FOR  ALL  SURFACES  REQUESTED. »// »40X • 12HDATA  IS  FOR  ,2Al0t 

2//) 
440    FORMAT  { lhO »8X»1 15HID      X       Y         Z       Z-LINEAR   RES-L 

1  Z-QUAD*   RES-QD   Z-CUEIC   RES-C   Z-QUAR.    RES-GR   Z-QUINT    RES- 

2QT/1X) 
445    FORMAT  ( 1H  , A10 *2F8 .2 *F9.3 *  5 t F9.3 *F9.3 5 ) 
45G    FORMAT  ( 104HTHE  LIMIT  PARAMETER  IS  ZERO  FOR  THIS  GROUP.  THIS  MEANS 

1  NO  FIT  WAS  ASKED  FOR.  GOiNG  ON  TO  NEXT  GROUP.     ) 
455    FORMAT  {  1H1  ) 
460    FORMAT  { 1H  »20X»8A10) 
465    FORMAT  (8A10) 

470    FORMAT  UH0»30X»  9HDATA  FOR  »2A10»/) 
END 


£  5y£R2yTiN§_sEA2i .«„«„— 

C  PURPOSE  - 

C         TO  READ  ID,X»Y»ZltZ2»Z3,Z4  DATA  FROM  THE  DATA  CARDS* 

C 

C         THE  DATA  CARDS  SHOULD  BE  PRECEDED  BY  A  FORMAT  CARD  (COLS  7-78) 

C  AMD  BE  FOLLOWED  BY  A  BLANK  CARD.  THE  SUBROUTINE  DETERMINES  THE  NUM- 

C  BER  OF  POINTS  IN  THE  DATA  SET  BY  SEARCHING  FOR  A  CARD  CONTAINING 

C  BLANKS  IN  COLUMNS  1-10  „ 

C    USAGE  - 

C         CALL  READKNCARDS) 

C 

C       NCARDS  -  THE  NUMBER  OF  DATA  POINTS  BEING  ANALYZED*  THIS  IS  DETF"< 

C  MINED  BY  READ1» 

C 

SUBROUTINE  READ1  (NCARDS) 

DIMENSION  X(l'CO),  YUJOO),  Z<1000»4>»  10(1000) 

DIMENSION  FMT(12) 

COMMON  X»Y,Z,!0 

DATA  NRLANK/10H  / 

1=1 

READ  (5*20)  FMT 
5      READ  (5»FMT)  ID( I  )  *X(  I ) »Y( I } »Z ( I «1 ) »Z (  I  9? ) «Z ( I t3) tZ ( I »4) 

WRITE  (1)  ID(I)»X(I)»Y(I)»Z{I»1)»Z{I»2)»Z(I»3)»Z<I.4) 

IF  (ID(I J-NRLANK)  10,  15 1 10 
10     1=1+1 

GO  TO  5 
C 

c 

15     NCARDS=I-1 

RETURN 
C 

c 
c 

2  0     FORMAT  (6X,12A6> 
END 


C  ROUT  I  ME  INVtR 

C  PURPOSE  - 

C  TO  INVERT  THE  XX  MATRICES  FOR  ALL  SURFACES  BEING  FITTED, 

C 

C  USAGE  - 

C         CALL  INVER(XX,XINV»K,NDEGRF) 

C 

C       XX  -  MATRIX  FOR  WHICH  I NVFRSF  IS  REQUIRED, 

C       XINV  -  THE  CALCULATED  INV'PRSF. 

C       K.  -  NUMBER  OF  TERMS  IN  EGUATION  BEING  DETERMINED=SIZE  OF  MATRIX* 

C       NDEGRE  -  DEGREE  OF  EQUATION  BEING  DETERMINED. 

C 


SUBROUTINE  INVER  ( XX»X I NV»K »NDEGRE ) 

DIMENSION  XX(21*21)»  X  I NV ( 21 .21 ♦ 5 > ,  AM(21»22> 

DIMENSION  N(21),  DUM'^Yf  12310  > 

COMMON  DUMMY »<SOB 

DO  5  I  =  1»K 

DO  5  J=1»K 
5      AM(I,J)*XX(ItJ) 

K.K.=K+1 
C 

C      INITIALIZE  THE  ROW  PERMUTATION  CELLS  TO  ZERO. 
C 

DO  10  1*1 »K 
10     NCI)«0 

DO  85  J-ltK 

I  =  1 
JJ=J-1 

15        DO  2o  L=1»JJ 

IF  (NCL)-i)  20,25»2C 
20        CONTINUE 

IF  (AM(I»1)»  35,25»35 
25        1=1+1 

IF  (K-I)  30,15,15 
30        KS0B=1 

RETURN 
3  5        AMA«ABS( AM(  I  »i  )  ) 

II  =  1 

40         11*11+1 

IF  (K-II1  60,45,45 
45        DO  50     IZ*1,J 

IF  (NUZ)-II)  50,40,50 
5  0        CONTINUE 

AMB=A6S(AM( 11,1)) 
C 
C      TEST  FOR  MAXIMAL  SIZF  PIVCTAL  ELEMENT. 

r 

IF  (AMA-AMB)  35,40,40 
5  5        AMA=AMB 
1  =  11 

GO  TO  40 


c 
c 

60        N(J)=I 

C 

C      ADJOIN  THE  CORRECT  IDENTITY  COLUMN  VECTOR. 

C 

DO  65  IZ=1,K 
65        AM(IZ.K.K.)=T90 

AMU»KK)el.O 

DA=AM(I t 1) 
C 

C      NORMALIZE  THE  PIVOTAL  ROW. 
C 

DO  7u  IZ=2,KK 
70        AM(I»IZ-1)=AM{Z.IZ)/DA 
DO  85  IZ=1SK 

IF  (IZ-I)  75,85,75 
75        DR=AM(IZ»1) 
C 

C      ZERO  OUT  MON 'PIVOTAL  COLUMN  ELEMENTS. 
C 

DO  80  JZ=2,KK 
9  0        AM(IZtJZ-l)*AM( lZ,JZ)-AM(I,JZ-l)#Db 
8  5     CONTINUE 

DO  100  1=1, K 

DO  00  J=1,K 

IF  (N(J)-l)  90,95,90 
90        CONTINUE 
95     DO  100  J.  =  1»K 

LL=N>  L) 
ICO    XINV(L,* ,NDEGRE)=AM(LL,J) 
RETURN 
END 


C  SUBROUTINE  INVCK. 

C  PURPOSE  - 

C         TO  DETERMINE  AND  PRINT  OUT  THE  PRODUCT  OF  THE  XX  AND  XX  INVERSE 

C  MATRICES  FOR  ALL  SURFACES  BEING  FITTED.  THIS  SUBROUTINE  CALLED  ONLY 

C  WHEN  INCK  (COL  25>MASTER  CONTROL  CARD)  IS  SET  EQUAL  TO  1  (SEE  DES- 

C  CRIPTION  OF  CONTROL  CARDS  PAGE     I  . 

C 

C         THESE  LISTINGS  CAN  r?E  ExAMINFD  TO  DETERMINE  IF  STRONG  DISCREP- 

C  ANCIES  FROM  AN  IDENTITY  MATRIX  OCCUR  FOR  ANY  SURFACE,  IF  DISCREP- 

C  ANCIES  OCCUR?  THE  COEFFICIENTS  SHOULD  3E  REGARDED  WITH  SUSPICION  » 

C  USAGE  - 

C         CALL  INVCK{A»B„K) 

C 

C  A  -  ORIGINAL  XX  MATRIX* 

C  B  -  INVERSE  MATRICES  FOR  EACH  DEGREE  SURFACE. 

C       K  -  NUMBER  OF  TERMS  IN  EQUATION  OF  EACH  SURFACE=SIZE  OF  IDENTITY 

C  AND  OTHER  MATRICFS. 

C 

SUBROUTINE  INVCK  (A9B»K) 
DIMENSION  A(21t21)»  P{21»2U5)»  C(21) 
N=( «K+l)*(K+2) )/2 
WRITE  <6»15J  K 
DO  10  1=1 »N 
DO  5  J  =  ls>N 

C« J)=Ge0 

DO  5  L=1*N 
5         C(J)=A{I$.L}*B(L*J»K)+C(J) 
10     WRITE  J(S»20)  I  S(C(J),J=1»N) 

RETURN 
C 
C 

c 

15     FORMAT  (1H1»77H  THE  PRODUCT  OF  THE  COEFFICIENT  MATRIX  AND  ITS  INVE 

1RSE  FOR  THE  FIT  OF  DFGREF  6I2»//MH  ROW) 
20     FORMAT  ( 1H0 » I3»2X»1 1FK„7»/6X» 10F10.7 ) 
END 


C  SUBROUTINE  POLY 

C  PURPOSF  - 

C         TO  EVALUATE  THE  EQUATIONS  OF  THE  REQUESTED  SURFACES  AT  EACH 

C  OBSERVED  .'ATA  LOCATION  IN  ORDER  TO  DETERMINE  THE  ESTIMATED  VALUES 

C  AND  THF  RLSIDUALS  FOR  EACH  DATA  POINT  . 

C    USAGF  - 

C  CALL  PCLY(ZR»LIMIT,XDD,YDD,ZDD} 

C 

C  ZR  -  UP  TO  TEN  VALUES  WHICH  ARE  THE  CALCULATED  AND  RESIDUAL  VALUES 

C  FOR  THIS  DATA  POINT  FOR  ALL  SURFACES  REQUESTED. 

C  LIMIT  -  PARAMETER  SPECIFYING  WHICH  SURFACES  ARE  BEING  FITTED* 

C  XDD  -  X  COORDINATE  OF  THIS  DATA  POINT* 

C  YDD  -  Y  COORDINATE  OF  THIS  DATA  PO: 

C  ZDD  -  Z  COORDINATE  OF  THIS  DATA  POINT* 

C 

SUBROUTINE  POLY  (ZR »LIMI T* XDD,YDD ,ZDD J 

DIMENSION  A(21,5J»  ZR(10),  XD<1)»  YD!D,  ZD(1)»  PTU) 

DIMENSION  XDOtl),  YDD(l)t  ZDDC1) 

DIMENSION  LIMIT(5)»  DUMMY* 12205) 

COMMON  DUMMY t A 

XD-XDD 

YD=YDD 

ZD=ZDD 

no  40  K=lt^ 

IF  (LIKITCK) )  5*40,5 
5         GO  TO  (10.15,20,25,30) i  K 
10        PT  =  A(  1»K)+XD*A(2»K)+YD*A(  3»K) 

GO  TO  2  5 
C 
C 

15        PT=A( l»2i+XD*A(2»2)+YD*A(3,2>+XD*XD*A(4»2)+XD*YD*A(5»2)+YD*YD*A 
1    (6,2) 

GO  TO  3  5 
C 
C 
20        XD2»XD*XD 

YD2=YD*YD 

PT=A(  l»3S<-XO*A(2»3)+YD*A<3,3)+XD2*A(4»3)+XO*YD*A(5»3)+YD2*A(6»3 

1  ) +XD2*XD*A (7,3 )  +  XD2* YD* A  < 8,3) +XD*YD2*A (9,3) +YD2*YD*A (10,3) 
GO  TO  3  5 

C 
C 
25        XD2=XD*XD 

XD3=XD2*XD 

YD2*YD*YD 

YD3sYD2»YD 

PT=A(  1*4)+XD*A(2»4)+YD*A(3»4)+XD2*A(4»4)+XD*YD*A<5,4)+YD2*A<6»4 

)+XD3*A(7.4)+XD2*YD*A{E»4)+XD*YD2*A<?»4)+YD3*A(l0»4)+XD2*XD2*A( 

2  H»4)+XD3*YD*A(12,4)+XD2*YD2*A(13»4)+XD*YD3*A(l4t4)+YD2*YD2*A{l 

3  5,4) 

GO  TO  3  5 


c 

30        XD2«XD*XD 
XD3=,".D2*XD 
XD4=^D3*XD 
YD2*YD*YD 
YD3=YD2*YD 
YD4=YD3*YD 

PT=A ( 1 9 5 ) +XD*A ( 2  9  5 ) +YD*A (3,5) +XD2*A (4 » 5  ) +XD*YD* A { 5 , 5 ) +YD2*A ( 6  » 
1    )  +XD3*A  {  7  „5  }  +XD2*YD*A  f  8  9  5  )  +XD*YD2*A  {9,5)  +YD3*A  (  10»  5  ) 

PT=PT+XD4*A{H»5)+XD3*YD*A(  12  »  3  J+XD2*YD2*A  (  13  ,5  )  +XD*YD3*A  ( 14,5 

1  +YD4*A<  I5,5)-i-xn4#Xr,>*Ai  16 »5 )+XD4*YD*A( 1 7 ,5 ) *XD3*YD2*A<  18,5)+XD2- 

2  YD3* A { 1 9  9  5 ) +XD*YD4*A ( 2C  »  5 ) + YD*YD4*A ( 2 1 » 5 ) 
35        ZR(?*K-1)*PT 

ZR(2*K>=ZD-PT 
40     CONTINUE 
RETURN 
END 


C  S'JtiROUTlNh    EVALU 

C  PURPOSF  - 

C         TO  FVALUATE  THE  EQUATIONS  OF  THE  SURFACES  AT  A  SERIES  OF  EQUAL- 

C  LY  SPACED  (GRIDDED)  LOCATIONS  SO  THAT  CONTOUR  MAPS  OF  THE  SURFACES 

C  CAN  BE  PRLPARED, 

C  USAGE  - 

C  CALL  EVALU(A,K»ISURF,TITL»NLINES»0RDIN»JIM) 

C       A  -  MATRIX  CONTAINING  POLYNOMIAL  COEFFICIENTS 

C       K  -  COUNTER  USED  TO  IDENTIFY  SURFACES  IN  MAP  TITLES* 

C       ISURF  -  PARAMETER  USED  TO  SPECIFY  MAPPING  OR  GRIDDING  OPTIONS. 

C       TITL  -  GENERAL  ALPHANUMERIC  TITLE(S) .NUMBER  SPECIFIED  BY  NLINES 

C       NLINES  -  NUMBER  OF  TITLE  LINES  «  NUMBER  OF  CARDS) »  MAXIMUM"  5  « 

C       ORDIN  -  SPECIFIC  ALPHANUMERIC  TITLES  USED  TO  IDENTIFY  VARIABLES- 

C       JIM  -  COUNTER  USED  TO  IDENTIFY  WHICH  VARIABLE  Cl-4)BEING  ANALYZED 


SUBROUTINE  EVALU  ( A >K, I S URF ,TI TL.NLI NES .ORDIN - 

INTFGER  COLS, ROWS 

DIMENSION  ISURF(5i 

DIMENSION  A(21r5)9  C{lCO»lOO)»  DUMMY (2205) 

DIMENSION  TITL (8,5 },  ORDIN 14*2) 

COMMON  C .DUMMY , A 

READ  (5.9.3  XMAX.XMIN  .YMAX.YMIN 

READ  (5.95)  ROWS.COLS.NMAP* I  PUN 

XlNC=(XMAX~XMIN)/FLOAT{COLS~i) 

Y I  NC=  (  YMAX-YMI  N  )  /FL.OA T  (  ROWS-I  > 


c 

c 

ZERO  C  ARRAY 

c 

DO  5  KK=1»95 

DO  5  JJ=1»95 

5 

C( JJ»KK)=0.0 

c 

c 

DEFINE  STARTING 

c 

X  AND  Y 

XD=XMIN 
YD=YMAX 
DO  45  JJ=l»ROWS 

DO  4J  KK=l»COLS 

GO  TO  (10,15,20,25,30),  K 
10  PT=A(l,l)+XD*A(2tl)+YD*A(3»l) 

GO  TO  35 
C 
C 

15  PT«A(1!2)+XD*A(2»2)+YD*A(3»2)+XD*XD*A(4»2)+XD*YD*A(5»2)+YD*Y 

1       D*A(6,2) 
GO  TO  35 


c 

?0         xn?=xo*xo 

YD2=YD*YD 

PT=A ( 1 , 3 ) +XD*A ( 2 ,3 ) +YD*A (3*3) +XD2*A (4,3) +XD* YD* A (5,3) +YD2*A ( 

1  6»3)+XD2*XD*A(7,3)+XD2*YO*A(8»3)+XD*YD2#A(9,3)+YD2*YD*A(lO*^ 

2  ) 

GO  TO  33 
C 
C 
25  XD2=XD*XD 

XD3-XD2#XD 

YD2=YD*YD 

YD3=YD2*YD 

PT=A ( 1 ,4 )  +XD*A( 2 ,4  >+YD*A  (  3  »4  }  +XD2*A  (4,4)  *XD*YD#A  (  5  ,4  J  +YD2*A  ( 

1  6»4)+XD3*A(7,4)+XD2*YD*A(8>4)+XD*,rD2*A(9.>4)-»-YD3*A(  10,4  J+XD2* 

2  XD2*A{  11,4)+XD3*YD*A(12,4)+XD2*YD2*  M  13  «4  J+XD*YD3*A(  14,4)+YD 

3  2*YD2*A(15,4) 
GO  TO  35 

C 
C 
30  X02=XD*XD 

XD3=XD2*XD 

XD4=XD3*XD 

YD2=YD*YD 

YD3=YD2*YD 

YD4=YD3*YD 

PT  =  A{1  ,5)+XD*A(?,5)+Yni*A(3»5  )+X02*A  {  4,5  ) -i-X'">*YL*M  5  ,5  )+YD2*A  ( 
1       6»5)+XD3*A(7,5)+XD2*YD*AC8t5)+X0*YD2*A(9,5)+YDi^AU0i5) 

PT=PT+XD4*A{  11,5  KXD3*YD*A(  12 , 5 ) +XD?#YD2*A  (  13,5  >+XD*YD3*A ( 14 

1  , 5  )  +YD4*A  ( 1 5  ,  5  )  +XD4*XD*A  (16,5)  +XD4*  rD*A  (  17,5}  +XDi*:-YD2*A  (13,5 

2  >+XD2*YD3*A( 19,5)+XD*YD4*A(20»5)+YD*YD4*A(21,5) 
35           C(JJ,KK)=PT 

40        XD=XD+XINC 

XD«XMIN 
45     YD=YD-YINC 

WRITE  (6,130) 

DO  50  I=1,NLINES 
30     WRITE  (6,135)  ( TI TL ( 1 1 , I ) , I  1  =  1  ,8  ) 

WRITE  (6,145)  ORDIN(JIf'.,l)  ,ORDlN  ( JIM,2  ) 

WRITE  (6*100)  K»COLS,ROWS 

WRITE  (6»110)  XMAX»XMIN»YMAXtYMIN 

DO  55  JJ=l,ROWS 

V/RITE    (6»115)     (C(JJ,KK>  SKK=1,C0LS) 
55  WRITE    (6»120) 

IF  (IPUN)  60,75*60 
50     DO  65  I=1»NLINES 
65     PUNCH  l-.O,  CTITL(IItI).II*l»6) 

PUNCH  1L,U,  ORDIN(JIM,l) ,0RDIN(JIM,2) 

PUNCH  105 »  K, COLS, ROWS 

DO  70  J=l,ROWS 

PUNCH  125,  (C( J,KK) ,KK>l,COL$) 
70     CONTINUE 


75     CONTINUE 

IF  (ISURF(K)«Nc.2.0R.COLS«GT.42)  GO  TO  85 
DO  80  I*1»NMAP 
80     CALL  PLOT  (Cf ROWS»COLS»K »TITL,NLlNES.ORDI N» JIM) 
85     CONTINUE 

RETURN 
C 
C 
C 

90  FORMAT  «4F10»3) 
95  FORMAT  t  3 I 3  *  I  1  1 
100    FORMAT  (1H0»36HGRIDDED  VALUES  OF  SURFAl       [GREE  »I1»9H«  GRID 

UI3.12H  COLUf  ROWS.) 

105    FORMAT  ( 3X»36HGRIDDED  VALUES  OF  SURFACE  OF  DEGREE  »Il,9H.  GRID 

1I3,12H  COLUMNS  8Y  SI3»6H  ROWS, ) 
110    FORMAT  (1HC55HMA)  COORDIN/       ...JFS  READ  U 

1MAX*»F10«3»6H  XMIN«»F1C.3»6H  V'MAX~»F10.3 »6H  YMIN»»F10.3// ) 
UK  »12F10«3) 
(1H0) 
(8F10.3) 
(  1H1  ■• 

(  1H  ,2r;X»8A10) 
(8A10) 

(1HC»1--X»13HANALYSIS  FOR  »2Al0  «•///// > 
(2vJX*2Al0) 


115 

FORMAT 

120 

FORMAT 

12  3 

FORMAT 

130 

FORMAT 

135 

FORMAT 

140 

FOPMAT 

145 

FORMAT 

15  0 

POPMAT 

„_„ -.- SUEEQUIJ       r..„ .——.-. -™— _-« 

PURPOSE  - 

TO  DEVELOP  PRINTER  CONTOUR  '-"APS  OF  5TED  SURFACES. 

SUBROUTINE  IS  <K    SPECIAL  VERSION  OF  SUBROUTINE  PLTCON1  .  THUS  THE 
SAME  BASIC  DESIGN  IS  RETAINED  AND  ACCORDINGLY  THE  INSTRUCTIONS  FOR 
PLTCON1  •  FOUND  ON  PAGE     »  SHOULD  BE  FOLLO  . 

JI RED  SUBROUTINES  — 

STATS-  TC  CALCULATE  CERTAIN  STATISTICAL  MEASURES. 

CONTROL  CARD  ORDER  — 

TH*S  >  ITRCL  CARDS  . 

\STER  CONTROL  C 
2) SPECIAL  CARD(S)  -  REQUIRED  ONLY  IF  TOUR  *  1  OR  2  • 

r  FUNCTIONS  OF  IRIBED  ON  PAGE 


ROUTINE    PLC 


>WS»C0LS»K0LD9TJTL 


INES»ORDlNv 


D I MENS I  -  131)  ,                              .31) , 

1    SYM  i  1  3 1 )  •  L  (  8  ,  5  )  »  ORD I N  ( 4  •  2 ) 

INTEGER  COLMAXs  .       M,PS»GRIDdSYM»bLAN  , 

IMINUSsSTAF!  »VSTOP 
REAL  MINZ»! 

DATA  (SYMT  ), lHl »1H2 » 1H3 s 1 H4 • 1H5  » 1H6  »1H7 » 1H8  » 1H9 , 1H 
lA»lHB»lHCtlHDPlHP,lHFjlHG,lHH,lHl»lHJ»lHK/»BLANK/lH  /.PLUS/1H+ 

S/lH-/i  1*/ 

.j  15.170)  CON»TOUR8LINES»lN< 
JJ=0 
-Q 

L  S 

=  NUMBER  OF  DATA  POINTS 
■is-************ 


1C 


WRITE  (6*175) 
10  I=1»NLIN 

i  TE  ( 6  » 1 80  )  ( T I TL (I  I » I  )  » II  =  1 ,8  ) 
WRITE  (6»185       N(JIf'»l)»0RDIN(JIM»2) 

fE  [6s>190)  ows 

CALL  STATS  t  ROWSs>COLS»10u  »100,Z,ZBAR  »XI 


c  NE    P'LLV 

C  TOUR    VALUE    DETERMINES    ME  5ED    IN    DEF  D.NTOURS 

C      ********************  #*-i?******.s##*.5j 

c 

IF  tTOUR.GT.3„OR.TOUR.LF.O)  GO  TO  30 

GO  TO  (15*20,25) ■  TOUR 
lb  WRITE  (6  „  195} 

READ  (5»2  MAX 

GO  TO  35 
C 
20     READ  (5*20  5?  INCON 

WRITF  (6,210} 

READ  (5, INCON)  (C        » 1=1, CON) 

ZMIN»DRANG{1) 

ZiMAX=DRANG(CC 

GO    TO    35 
C 

25  ZMIN«ZR*,R-3.0*ZS1 

ZMAX*ZBr,R  •'  IG 

COM =1? 

WRITE  (6,215  5 
TO  35 
C 
30     WRITE  (6,220) 

ZMIN=MINZ 

ZMAX=MAXZ 
35        TE  (6,225)  R0WS,C0LS?ZBAR,2SIG,MAXZ,MINZ 


r 

********* 

C    CALCUL  (COL) 

****** 
C 

IF  UNCHESaLE.0.0)  COL- 

IF  (INCHES0GT»i29        127 

IF  (INCHES«»GT.O«0»AN       [  :S»LE»12«7 3       ^CH^S^lO 

PS~ ( COL-COLS  J / ( COL  S 

c 

C      **********  *******i 

TEST  GRID  SIZE— IF  GRID  SMALLER  THAN  2  CHARACTERS  BETWEEN  COLUMNS 
C    AN  ERROR  STATEMENT  IS  PRWil'ED  fv.MD  PLC        ED 

*************  »******#^*J!-*#4i 

c 


IF  (PS.GT.l)  GO  TO  40 
WRITF  (6,2: 
GO  TO  165 


******** 

C    CALCULATE  RE        CONST.-. 

****************«.#**##4J.. 

40     COLMAX=PS* I  COLS-  I  • +COLS 


YlNC=ie^/( ( FLOAT (PS+1) )*0.6) 

GRID=PS+1 

XINC*1.0/FLOAT«6RIO) 

RANG2=ZMAX-ZMIN 
RHIGH=ZMAX+RANG2 

RLOW=ZwlM-RANGZ 

IF  (C0N.EQ.0.OR.C0N.GT.19)  CON  =  5 

FCON=CON 

CI=RANGZ/FCON 

IF  (TOURe£Q»2)  GO  TO  50 

JJJ=C0N+1 

no  45  T=1»JJJ 
45     DRANG(I)=ZMlN-5-FLOAT(I~l)*CI 
50     CONTINUE 

WRITE  (6.235) 

WRITE    (6.240)     ZMIN,ZMAX,CI  *  GRID, XI  NC  ,  YlNC  ,CON  t>PS.C0LMAX 
C 
C  #**»###*#*#####*#**•*♦##*###***#****■#*#* 

C         WHITE    OUT    CONTOUR    SYMBOL    TABLE 

c 

WRITF    (os>2A5 ) 

IF  (LlNESoFQ.l)  WRITF  (6.250) 

WRITF  (6*255) 

IF  (LlNES.NE.n)  GO  TO  55 

WRITE  (6s260> 

WRITE  (6,265)  (SYMTABi I )tDRANG(I)f I»19C0N) 

WRITE  {6 ,270) 

GO  TO  60 
C 
5  5     WRITE  (6,275) 

WRITE  (6.280}  (SYMTA9( I ) ♦DRANG f I ) .DRANG { 1+1) ♦I=19C0N) 

WRITE  (6.285)  DRANG  (CON-*  1 ) 
60     WRITE  (6.175) 

DO  65  I=1»NLINES 
65     WRITE  (6.180)  < TITL ( 1 1 » I ) *  1 1  =  1 .8 ) 

WRITF    (6.185)    ORDINUIN'sl)  *0RDIN(JIM,2) 

WRITE    (69190)     K OLD, N, COLS, ROWS 

XSTOP=GRID 

GRID=C0LS-1 

C0L=3 

DO    75    J=1»GRID 

DO    7o    I=29XSTOP 
SYM(COL)=ELANK 
70  C0L=C0L-M 

SYM(COL)=PLUS 
75  COL=COL+1 

M=C0LMAX+2 

SY.'M2)=PLUS 

SYM(1)=BLANK 

SYM(N)=dLANK 

WRITE    (6.290)     ( SYM( I ) » I =1 ,N> 


C    ESTIMATE  VALUE  OF  ALL  PRINTING  POSITIONS  ON  MAP 

C         ##***#**•***#•!!■*#*■#*#*#*»•»#**  ***■**•»#*«#** 

C      LOOP  THROUGH  MAP  LENGTH ,  ROW  BY  ROW 

ITOP*l 

IRO^sROWS+1 

DO  160  MM«2tIROW 

YSTOP-L/YINC+«,5 

IF  (MM.EQ.IROW)  YSTOP-1 

M*MM-1 

YP=0*j 

PSYMC 1)=PLUS 

PSYM(N)=PLUS 
no  160  KJ*1,Y^T0P 

C0L=1 

Z 1  * YP*  (  Z  (  MM  »  1 )  -Z  <  !*  »  1  )  >  +Z  {  M  , 1  ) 
C      #*#»######*#**####*##*###■****»#*•###*#*# 
C      LOOP  THROUGH  MAP  WIDTH .COLUMN  by  COLU 

C         *#***«#*#■*##*#*#«#########*##**#*:#»#*•*# 

II=C0LS-1 

JLFT=2 

DO  1*5  J*1*II 
J.*»J+1 

Z--YP*(Z(Mvt,JJ)-Z<  M,JJ)  )+Z(M»JJ) 
ZDELT=Z?-Z1 
XP*0.0 

DO  130  KK«1»XST0P 
JLFT=JLFT-1 
C0L=C0L+1 
ZEST=XP*ZDELT+Z1 

IF  (ZEST.GE8ZMlN.ANDaZFST«LF«>ZMAX)  GO  TO  95 
IF  {ZEST.LT.ZMINoAND. { I TOP.FG. I  *  OR. MM* EG. IROV/.ORe JLFT»EQ« 
1  1)}  GO  TO  80 

IF  (ZEST»GT.ZMAX»AND. { iTOP.EU.l*OR«Mi«i.EQ. IROW.OR» JLFT.EQ- 
1  1)!  GO  TO  85 

PSYM(COL)=BLANK 
ZPREy(C0L)=ZE5 T 
ZLAST»ZEST 
GO  TO  130 
C 

80  IF  (ZEST.LE.RLOW)  GO  TO  90 

PSYM(C0L)=MINUS 
ZPREV(COL)»ZEST 
ZLAST*ZEST 
GO  TO  130 


85  IF     (Z£SToGE„RHlGr>)    60    TO    90 

PSYM(C0L)=PLU5 
ZPREVlCOL)=ZEST 
ZLAST=ZEST 
GO    TO    130 
C 

90  P5Y*1(C0L)=STAP 

ZPREV(COU~ZEST 
ZL  \ST  =  ZEST 
GO    TO    130 
C 

95  I=<<ZEST-ZMIN)/CI >+lo 

LVAL=LINES+1 

IF  { IT0P.FQ.1.0R.MNUFG.IR0W.0R.JLFT-E0.1)  LVAL=3 
GO  TO  (100,120, 125) t  LVAL 
100  DO  105  1*1, CON 

ZSTEP=DRAMG(I } 
SW1=1» 
SW2=1* 
SW3*1. 
SW4*1« 

IF  (ZSTEP.GE.ZLAST.OR.ZSTEP.GE.ZPRFVCCOL))  Swl  =  0. 
IF  (ZSTEP.LE.7E5T)  SW2*0. 

IF  (ZSTEP.LT«ZLAST.OR.ZSTtP.LT»ZPREV(COD)  SW3*0» 
IF  (ZSTEPoGE.ZEST)  SK4=C* 

IF  (  (SWl»EO»C««AND»SW2*ECJ«Q.  )  »GR®  {  SW3eEu«0«  «  AND.  SW4»EG 
1  .0.) )  GO  TO  110 

105  CONTINUE 

PSYM(COL)=BLANK 
GO  TO  115 
C 

110  PSYM(COL)=SYMTAB(I ) 

115  ZLAST=ZEST 

ZPR£V(COL)=ZEST 
GO  TO  130 
C 
120  K-I/? 

IF  <(K*2).NE.I>  PSYM(C0L)»SY^1A8(I> 
IF  {  «K«-2JoE06!  3  PSYM(COL)«BLANK 
GO  TO  130 
C 
12  5  PSYM(COL)*SYMTABU) 

ZPREV(COL)=ZF"ST 
130  XP=XP+X1NC 

7l=Z2 
135       JLFT=0 

C0L=C0L+1 

IF  (Z{yn.COLS)«GE«ZMIN.AND-Z(M»COLS)«LE.ZMAX>  GO  TO  150 

IF  (Z(Mt.COLS)*GT*ZMAX)  GO  TO  140 

IF  (Z<rt»COLS).LE.RLOW)  GO  TO  145 

PSYM( COL) 'MINUS 

GO  TO  155 


c 

140 


c 

145 

C 
150 

155 


160 
165 

C 

170 
175 
180 

185 
190 


195 
200 
205 
210 
215 
220 
225 


230 
235 
240 


245 
250 
255 
260 
265 
270 
275 
280 

28  5 

29  0 


IF  (Z{M,COLS).GE.RHIGH)  GO  TO  145 

PSYM<COL)=PLUS 

GO  TO  155 

PSYM(COL)=STAR 
GO  TO  155 

I  =  « ( Z ( H ,COLS ) -ZM I N ) /C I ) +1 • 
PSYM(COL)=SYMTAB(I ) 
WRITE  (6>290  3  (PSYM(L) »L*1«N) 
YP=yP-J-YINC 
PSYM( 1)=BLANK 
PSYM(N)=BLANK 
IT0P=U 

WRITE  (6*290)  (SYM(L) »L=1»N) 

CONTINUE 
RETURN 

FORMAT  (I2,2I1,F4.0) 

FORMAT  flHl) 

FORMAT  ' 1H  »2ox,8A10) 

FORMAT  i  lHr»l''X»21HANALYSIS  OF  DATA  FOR 

FORMAT  ( 1HO»33HCONTOUR  MAP  OF  SURFACE  0 

HON  HAS  BEEN  EVALUATED  AT  *I4»2lH  LOCAT 

2MNS  BY  tI396H  ROWS.//) 
FORMAT  (25H  CONTOUR  LEVELS  SPECIFIED) 
FORMAT  (2F10.0) 
FORMAT  C10A6J 

FORMAT  (37H  VARIABLE  CCNTCuR  INTERVALS 
FORMAT  (3uH  CONTOUR  INTERVAL  =  Ul/2  SIG 
FORMAT  {26H  CONTOUR  LEVELS  CALCULATED) 
FORMAT  U4H  GRIDDED  DATA/*42H  LINEAR  I  ,\ 

1ARES  ///»17H  OBSERVED  VALUES  //»8H  R0.4S 

?M    ZBAR  =  sF10.3s.10H    ZSlG  ■  tFlu.3»luF 

3  =  9F10.3) 
FORMAT  (25H  TOO  MANY  COLUMNS  IN  GRID) 
FORMAT  (1H  /,16H  VALUES  EMPLOYED) 
FORMAT  <8H  ZMIN  =  .F10.3.I0H    Z«MAX  ~    * 

1GRID  =  »I5,10H    XINC  =  ,F10«3»lV/H    Yl 

2H    PS  -  »I5»12H    COLMAX  =  »I5) 
FORMAT  (35H4CONTOUR  SYMBOL  TABLE  IS  AS 
FORMAT  (43HCX0NLY  EVEN  SYMBOLS  PRINTED 
FORMAT  (7H0  SYMBOL  »HX  «  5HVALUE) 
FORMAT  (5HC    *»6X»7HM0  DATA,/,«5H 
FORMAT  <4X»A1»2X*F12,4) 
FORMAT  (4X»1H+»4X.10HAE0VE  HIGH) 
FORMAT  (5Hu    *.25X»7HMO  DATA,/,5H 
FORMAT  (4X9A1»2X,F12.496H   TO   9F12,4) 
FORMAT  (4X,1H+,2X,F12.491CH   OR   OVER) 
FORMAT  ( Ih9»129Al) 
END 


»2A10»///) 

F  DEGREE  ,I1,37H.  THE  EQUAT 
IONS  ON  A  GRID  ,  I3,12H  COLU 


SPECIFIED) 

M  4  ) 

TERPOLATION  W I THIN  GRID  SOU 

=  .I3.10H    COLS  =  »I3,/i8 

MAXZ  ~    »F10.3»10H    MINZ 


F1G„3,8H    CI  =  »Fl0e3/»8H 
NC  =  »F10.3,/,7H  CON  =  15,8 

FOLLOWS) 

IN  BODY  OF  MAP) ) 

.4X.10HBELOW   LOW) 


■»23x*lOHdELOW   LOW) 


C  JOB  CONTROL  CARDS  NECESSARY  TO  PUN  POLYFIT  AND  RFSMAP  AS  ONE  JOB 

C  COL  1 

C  ACCT»NAME*T800,CM70:00»PRIORITY 

C  RUN{ $»♦»»* ?77777} 

C  REQUEST (TAPE8, HI;    USER  TAPE  IDENTIFICATION   (THIS  CARD  OPTIONAL 

C  L60e 

C  RFL»7u-„ou. 

C  RFWINDILGO) 

C  RUN{Ss>  s>  »*>* 77777  ? 

C  LGO» 

C  7-8-9  CARD 

C  PROGRAM  PCLYFIT  SOURCE  DECK 

C  7-8-9  CARD 

C  POLYFIT  CONTROL  AMD  DATA  CARDS  CS^F  RFLOWJ 

C  7-8-9  CARD 

C  PROGRAM  RESMAP  SOURCE  DFCK 

C  7-8-9  CARD 

C  RESMAP  CONTROL  CARDS 

C  6-7-8-9  CARD 

C 


c 

c    - — — — 
C    MOTE  •- 

PROGRAM  PCLYFI7  -  SAMPLE 

DATA  Cifi 

:RDS 

'HE 

SE  CARDS 

WILL  FIT 

TRENDS  1- 

■4 » 

PRODUCE 

PRINTER  MAPSiP'JNCH 

c   rRinnrp 

TREND  VALUES  FOR  A 

1/2  INCH  GRI 

DsAND  SI 

ORE  DATA  FOR 

RESMAP 

C         FOI 

C 

C         THi 

JR 

DIFFERENT 

VARIABLES  WILL  BE 

ANALYZED. 

•  DATA  REPRt 

SENTS  WATER  WELL  OSS 

ERVATIOf 

IS  DURING  THE 

LAST 

C         FOUR 
C 

MONTHS  OF 

1967  IN 

BARTHOLOiV 

IEM 

COUNTY 

> INDIANA. 

c 

POLYFT 

POLYFT     1111 

22?? 

I             1 

1111 

1 

1 

BARTHOLOMEW 

COUNTY* 

INDIANA 

ANALYSI 

S  OF 

■ER-WELL  LEVELS 

SEPTEMBER 

1967     OCTOBFR   1967 

NOVEMBER 

19  67        DECEMBER 

{A109i 

2F1 

u.2»4F10. 

3) 

01 

00485 

00115 

597SOQ 

596600 

596600 

596700 

02 

00719 

CO  134 

590200 

589300 

591400 

592900 

03 

O0828 

v.  0  231 

60901*0 

607700 

608700 

61C9C0 

4 

00674 

00  253 

5930^0 

592400 

592100 

595500 

07 

00414 

C0369 

602600 

601700 

601100 

601100 

77 

00573 

00347 

592700 

592700 

592700 

592700 

08 

00684 

00331 

5964O0 

596000 

395900 

597900 

10 

00587 

00485 

600000 

599500 

5992C0 

600600 

11 

00727 

00570 

6099^0 

609300 

608800 

609000 

14 

0ul66 

00739 

61?2'.  0 

613200 

613200 

613200 

16 

00414 

00"?36 

611600 

611300 

611100 

612900 

18 

00762 

00666 

614500 

6139Q0 

613500 

614600 

20 

00467 

00845 

614500 

614300 

614400 

615300 

24 

00456 

01049 

623300 

623200 

626500 

626600 

26 

00759 

01116 

646C!'-'0 

644900 

644600 

646800 

27 

00990 

01229 

652200 

651800 

651700 

6519C0 

28 

00478 

01225 

634100 

633800 

633400 

635  500 

30 

U0221 

01163 

6296  00 

628800 

629100 

631700 

36 

00410 

0 1510 

6468  0  0 

645400 

644800 

644400 

37 

00292 

01^:2  0 

6437^0 

642  800 

642100 

641900 

38 

00414 

01417 

6455 0  0 

644400 

644100 

643800 

48 

00444 

00656 

6057o0 

604700 

604800 

604400 

50 

00706 

01360 

645300 

644200 

643900 

643800 

51 

00705 

01513 

651000 

650700 

650700 

652100 

S2 

00634 

00736 

617600 

617300 

617000 

617300 

93 

00630 

00761 

6128w0 

612900 

612700 

613700 

94 

00639 

01167 

635400 

63*4900 

634400 

634000 

93 

>17 

00740 

615400 

614900 

614400 

6144C0 

96 

00574 

00934 

627900 

.627400 

627000 

627000 

li«0 

0. 

0       16 

,0 

CO 

33  23   11 

loom. 2 

11.0 

0. 

0       16 

.0 

C*.0 

33  23   11 

•100111,2 

1967 


11.0  0.0  16.0  CO 

33    23       11 
100111,2 

11.0  0.0  16,0  CO 

33    23       11 
100111.2 

11.0  OoU  16.0  CO 

33    23      11 
100111.2 

11.0  0»0  16.0  CO 

33    23       11 
100111.2 

11.0  0*0  16. C 

33    23       11 
100111.2 

11.0  0.0  16.0  CC 

33    23       11 
100111.2 

11.0      O.u       16.0       0.0 
33  23   11 
100111.2 

li.o    co     16.0     r.o 

33  23   11 
100111.2 

11.0      CO       16.0       0*0 
33  23   11 
100111.2 

11.0      0.0       16.0       0a0 
33  23   11 
100111.2 

11.0      0.0       16.0       0.0 
33  23   11 
100111.2 

11.0      O.u       16.C       0,0 
33  23   11 
100111.2 

11.0  O.u  16.0  CO 

33    23       11 
10nlH.2 

11.0  CO  16.0  0«0 

33    23       11 
100111.2 
C 

C  PROGRAM  RESMAP  SHOULD  FOLLOW  IMMEDIATELY 


SAMPLE    OUTPUT   FROM    PROGRAM       POLYFIT  32A 


pnara>  m  w.h  >*  f-ti 


TREND  r"^E==~^£  TRENDS 


::TRI;ND:::^       '=     I  ,,: TREND  4 


C  F-'iC-**.***     RESMAp  -  GENERAL  INSTRUCTIONS 

C  PURPOSE  - 

C  TO  BE  USED  IN  CONJUNCTION  WITH  PROGRAM  POLYFIT  •  THIS  PROGRAM 

C  ANALYZES  THE  ORIGINAL  AND  RESIDUAL  VALUES  STORED  ON  TAPE8  6Y  PROGRAM 

C  POLYFIT.  GRIDDED  VALUES  AND  PRINTER  CONTOUR  MAPS  ARE  PRODUCED  BY 

C  SUBROUTINES  GRID  AND  PLTCON1  OR  PLTCON2  •  A  VARIETY  OF  CONTOUR  MAPS 

C  MAY  BE  SPECIFIED*  MULTIPLE  MAPS  OF  ANY  DATA  MAY  BE  PRODUCED  „  THE 

C  GRIDDED  VALUES  MAY  PRODUCE  CALCOMP  CONTOUR  MAPS  USING  PROGRAM  CONTUR 

C 

C  PURDUE  UNIVERSITY  VERSION  2»0      A,  KEITH  TURNER, 

C  CIVIL  ENGINEERING, 

C  FEBRUARY*  1968  • 

C  TAPE  REQUIREMENTS  - 

C         TAPE8   IS  USED  TO  STORE  DATA  DEVELOPED  BY  PROGRAM  °OLYFIT» 

C  SUBROUTINES  REQUIRED  — 

C         DGRID  -  TO  GRID  THE  IRREGULARLY  SPACED  DATA 

C         2)PRTCONl-  TO  PRODUCE  PRINTER  CONTOUR  MAPS 

C  OR 

C  PRTCON2-  TO  PRODUCE  GRADED  DENSITY  PRINTER  CONTOUR  MAPS 

C         3)STAT5-  TO  CALCULATE  CERTAIN  STATISTICAL  VALUES 

C 

C  NOTE  1  —  GRID  CONTROL  CARDS  ARE  REQUIRED  FOR  EACH  ANALYSIS  REQUEST- 
C  ED  •  SEE  SECTION  III. 

C  NOTE  2  —  PRTCON1  OR  FRTCON2  CONTROL  CARDS  REQUIRED  FOR  EACH  MAP 
C  REQUESTED  (NUMBER  OF  MAPS  =  NUMBER  OF  MAPS  PER  REQUEST 

C  TIMES  NUMBER  OF  ANALYSES  REQUESTED).,  SEE  SECTION  IV  . 


C  PROGRAM  RE5MAP  -  MAIN  PROGRAM 

PROGRAM  RES  tAP< INPUT. OUTPUT 9TAPE8,PUNCH5TAPE5=INPUT9TAPE6™OUTPUT > 
DIMENSION  X»100C),Y(10CO)*ZUOOO>,IVAL(6)  »TITL(8) 

5  READ(5,65)  TITL 
IF(EOF,5J  55,6 

6  READ(5,85)  { I VAL ( I ) • 1=2 ,6 ) 
DO  50  1=1,6 

IF(IVAL( I)*NE.l)  GO  TO  50 
READS8)  NUM 
IF(EOF,3)  55*10 

10  WRITE(6>60!  NUM 
WRITE(6,70)  TITL 

GO  TO  {2o921j22?23»24s25}  ,  I 

20  DO  11  N=1»MUM 

11  READ  (8) lD*XCN)»Y«N)»Z(N)»ZCl»ZRl»ZC2»ZR2fZC3tZR3»ZC4.ZR4tZC5»ZR5 
GO  TO  45 

21  DO  12  N=lfNUM 

12  READ  <8)IDtX(N)*Y{N)9ZD»ZCl»Z(N)  ,ZC2 5ZR2 ,ZC3 »ZR3*ZC49ZR4»ZC5,ZR5 
GO  TO  45 

22  DO  13  N=1»NUM 

13  READ  (8>ID»Xt?4>*Y<N),ZD»ZCl»ZRl»ZC2»Z<N)  »ZC3tZR3»ZC4»ZR4»ZC5tZR5 
GO  TO  45 

23  DO  14  N=1»NUM 

14  READ  (8)ID»X(N}»YCN)«ZDtZCl»ZRltZC2»ZR2»ZC3»Z(N)  »ZC4.ZR4»ZC5,ZR5 
GO  TO  45 

24  DO  15  N" Is NUM 

15  READ  C8>ID.X{N)»Y(N),ZD.ZCl,ZRl»ZC2tZR2»ZC3»ZR3»ZC49Z(N)»  ZC5eZR5 
GO  TO  45 

2  5     DO  16  N* It NUM 

16  READ  (8)ID»X(H)fY(N) jZD»ZCl ,ZR1„ZC2»ZR2»ZC3  9ZR3»ZC4,ZR4.ZC5»Z (N) 
45     REWIND  8 

CALL  GRID  (NUMfXtY»Z) 
50     CONTINUE 

GO  TO  5 
55     STOP 

60     FORMAT(1H1»10X.32H  NUMBER  OF  POINTS  TO  BE    READ  IN-»I4,////) 
65     FORMAT(8A10) 
70     FORMATC1H  ,5X»8A10J 
8  5     F0RMAT(6I1) 

END 


C  PROGRAM  RESMAP  -  SAMPLE  DATA  CARDS 

^  — _-^— >.„_„«^..« 

C  MOTE  -THFSF  CARDS  6RF  D^SIfWFD  TO  ANALYZE  DATA  PRGDUCEH  BY  PROGRAM 
C         POLYFIT  FROM  THE  SAMPLE  nATA  SHOWN  THFRF.  THESF  CARDS  RFQUFST 
C         THE  PRODUCTION  OF  ORIGINAL  AND  FOUR  DIFFERENT  RESIDUAL  MAPS  FOR 
C         THE  SAME  FOUR  MONTHS*  PUNCHED  CARD  OUTPUT  FOR  THE  GRIDDED 
C         VALUES  WILL  ALSO  BE  PRODUCED. 
C 

^  -„_.«..._„ 

BARTHOLOMEW  COUNTY  WATER  SURFACE  ANALYSIS  FOR  SEPTEMBER  1967 

33  23     101  1 

BARTHOLOMEW  COUNTY  GROUND  WATER-OR IG INAL  WELL  ELEV.  DATA .SEPTEMBER  1967 

11«0     0.0        16*0        0.0 
100111.2 

33  23     lui  1 

BARTHOLOMEW  COUNTY  GROUND  WATER-  1ST  DEGREE  RESIDUALS  *  SEPTEMBER  1967 

H»0     090  16. C        0*0 

100111,2 

33  23     101  1 
BARTHOLOMEW  COUNTY  GROUND  WATER-  2ND  DEGREE  RESIDUALS  *  SEPTEMBER  1967 

H»0     0.0        16.0        O.u 
100111,2 

33  23     101  1 
BARTHOLOMEW  COUNTY  GROUND  WATER-  3RD  DEGREE  RESIDUALS  ,  SEPTEMBER  1967 

11*0     0.0        16,0        0.0 
100111.2 

33  23     101  1 
BARTHOLOMEW  COUNTY  GROUND  wATcR-  4TH  DEGREE  RESIDUALS  •  SEPTEMBER  1967 
H»0     0.0        16.0        OoO 

loom. 2 

BARTHOLOMEW  COUNTY  WATER  SURFACE  ANALYSIS  FOR  OCTOBER  1967 
11111 

33  23     101  1 
BARTHOLOMEW  COUNTY  GROUND  WATER-ORIGINAL  WELL  ELEV.  DATA, OCTOBER  1967 

11*0     0.0        16.0        0eC 
100111.2 

33  23     101  1 
BARTHOLOMEW  COUNTY  GROUND  WATER-  1ST  DEGREE  RESIDUALS  ,  OCTOBER  1967 

11-0     0.0        16»0        0.0 
100.111.2 

33  23     lui  1 
BARTHOLOMEW  COUNTY  GROUND  WATER-  2ND  DEGREE  RESIDUALS  ,  OCTOBER  1967 

H»0  0.0  16.0  0.0 

100111.2 

33    23  1>11 

RARTHOLOMFW  COUNTY  GROUND  WATER-  3RD  DEGREE  RESIDUALS  ,  OCTOBER  1967 

H»0  0.0  16„0  0.0 

1001U.2 
33    23  101    1 


BARTHOLOMEW  COUNTY  GROUND  WATER-  4TH  DEGREE  RESIDUALS  *  OCTOBER  1967 

11.0     0.0        16,0        0.0 
100111,2 

BARTHOLOMEW  COUNTY  WATER  SURFACE  ANALYSIS  FOR  NOVFMBER  1967 
11111 

33  23     101  1 
BARTHOLOMEW  COUNTY  GROUND  WATER-ORIGINAL  WELL  ELEVe  DATA , NOVEMBER  1967 

11.0     o.J        16.0        0.0 
100111,2 

33  23     101  1 
BARTHOLOMEW  COUNTY  GROUND  WATER-  1ST  DEGREE  RESIDUALS  »  NOVEMBER  1967 

11.0     0.0        16*0        0.0 
100111,2 

33  23     101  1 
BARTHOLOMEW  COUNTY  GROUND  WATER-  2ND  DEGREE  RESIDUALS  *  NOVEMBER  1967 

11.0     0.0        16.0        0,0 
100111,2 

33  23     101  1 
BARTHOLOMEW  COUNTY  GROUND  WATER-  3RD  DEGREF  RESIDUALS  >  NOVEMBER  1967 

11.0     0,0        16.0        0.0 
100111,2 

33  23     101  1 
BARTHOLOMEW  COUNTY  GROUND  WATER-  4TH  DEGREF  RESIDUALS  *  NOVEMBER  1967 

11.0     0.0        16.0        0,0 
100111,2 

BARTHOLOMEW  COUNTY  WATER  SURFACE  ANALYSIS  FOR  DECEMBER  1967 
11111 

33  23     101  1 
BARTHOLOMEW  COUNTY  GROUND  WATER-OR IGTNAL  WELL  ELEV«  DATA .DECEMBER  1967 

11.0     Oe0        16.0        0.0 
100111.2 

33  23     101  1 
BARTHOLOMEW  COUNTY  GROUND  WATER-  1ST  DcGREE  RESIDUALS  »  DECEMBER  1967 

11.0     O.v,         l&e0        0.0 
100111.2 

33  23     101  1 
BARTHOLOMEW  COUNTY  GROUND  WATER-  2ND  DEGREE  RESIDUALS  »  DECEMBER  1967 

11*0     0*0        16,0        0,0 
100111.2 

33  23     101  1 
BARTHOLOMEW  COUNTY  GROUND  WATER-  3RD  DEGREF  RESIDUALS  »  DECFM3ER  1967 

11»0     0.0        16^0        0.0 
100111.2 

33  23     101  1 
BARTHOLOMEW  COUNTY  GROUND  WATER-  4TH  DEGREE  RESIDUALS  »  DECEMBER  1967 

11.0     ue0        16*0        0,0 
100111.2 
C 


SAMPLE   OUTPUT   FROM    PROGRAM      RESMAP 
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PART  III 

PROCEDURES  FOR  ASSIGNING 

IRREGULARLY  SPACED  (RANDOM)  DATA 

TO  A  SQUARE  GRID 


NOTE  -  This  procedure  most  precede  any  attempt  to  machine -contour 
irregularly  spaced  data.  Since  potential  users  may  wish  to  clearly 
understand  the  underlying  principles  ased  in  subroutine  GRID,  a 
description  of  the  algorithm  is  included. 


CONTESTS  OF  THIS  SECTION 

1)  Description  of  Interpolation 
Algorithm 

2)  SUBROUTINE  GRID 

3)  SUBROUTINE  STATS 
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AN  INTEHPQLATiasr  ALGORITHM  FOR  GEOGRAPHICAL  DATA 

The  documentation,  for  this  algorithm  and  a  MAD  language  program  were 
obtained  from  Professor  W»  R.  Tobler,  Geography  Department,  University  of 
Michigan.  The  MAD  language  program  formed  the  basis  of  the  Subroutine 
GRID.  The  algorithm  involves  a  weighted  moving  average  procedure,  and 
appears  to  give  results  which  are  comparable  to  those  obtained  by  the 
fitting  of  local  polynomials,  but  is  computationally  much  simpler.  Follow** 
ing  considerable  (empirical)  testing,  the  main  disadvantage  appears  to  be 
a  slight  dampening  (or  smoothing)  effect  on  the  data*  This  approaches 
10$  of  the  data  range  for  extremely  erratic  data,  but  for  commonly  en- 
countered data  distributions  is  much  less.  Competing  schemes  are  described 
in  the  references  cited. 

The  following  discussion  was  obtained  from  Tobler:  - 
Definitions: 

An  observation  is  recorded  as  X. ,  Y. ,  Z. ;  i  =»  1,  ...  ,  N,  where 
X,  Y  are  locational  coordinates  and  Z  is  a  value  at  S,  Y, 

Lattice  points  are  labelled  jk  in  standard  isatrix  notation  (j  =  row- 
subscript,  and  K  »  column  subscript).  The  problem  is  to  estimate  Z.,  . 

Procedure : 

(1)  Calculate  A,  the  area  of  the  rectangle  bounded  by  the  maxi- 
mum and  minimum  of  the  coordinate  values  X  and  Y. 

(2)  Let  the  grid  size,  G,  equal  the  square  root  of  the  area 
divided  by  the  number  of  observations, 

G  -  (A/K)l/2 

(3)  The  number  of  rows  in  the  matrix  is  computed  as  the  (truncated) 
integer  part  of 

R  *  (  A  Y/G)  +  1.5 


-  2  - 
and  the  number  of  columns  as  the  integer  part  of 

Cs(A  X/G)  +  1.5 
This  yields  a  matrix  which  has  approximately  N  entries . 
(k)       In  general,  because  of  truncation, 

(C  -  1)  •  G  ^  A  X 
(R  -  1)  -G/iY 
Split  the  difference,  letting 

A"  X  ■  (  A  X  «  (C   1)  •  G>/2 
A"  Y  «  {.  A Y  -  (R   1)  •  G)/2 
and  assign  X  and  Y  coordinates  to  positions  in  the  matrix „  For 

"bill  "fc5l. 

example  the  j   row  and  the  k   column  will  have  the  coordinates 


K    "  Xnrfr,  *  ^"  X  *  (k  -  1) 

k    mxn 


G 


Y .  «  Y    -  A!'  Y  -  U   -  1)  •  G 
3        max        *°    * 

(5)   For  each  position  in  the  mtrix  calculate  the  square  of  the 
distance  to  each  observation  using 

D^  -  (Xk  -  X±f  +  (Y.   -  Y..)2 

and  rank,  the  points  in  order  of  increasing  distance  from  the 

th 
jk   lattice  point  as  they  are  computed,  saving  the  results  for 

2 
only  the  nearest  six  points,  except  if  D.  <   e  in  which  case 

set  Z .     »  z.  and  proceed  to  the  next  position  in  the  matrix 

since  the  assignment  has  been  completed.  (We  have  taken  €  ■  G/25, 

but  this  is  arbitrary). 


-  3  - 
(6)  Assign  a  value  to  the  jk  l  position  by  using  a  weighted  mean 
of  the  value  at  the  nearest  point  with  the  weighted  mean  of 
the  nearest  six  points.  Symbolically,  using  the  notation 
IX  a  to  denote  the  rank  a  of  the  distance  from  the  jk 
position  to  the  point  i  and  2   for  the  observed  value  at  the 


nearest  point,  this  can  be  written  as 


z      +    a^i 

a 


6 

V     Jk, 
Dia 


z., 


L         ha 

asl 
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C  SUBROUTINE  GRID 

C  PURPOSE  - 

C  TO  ASSIGN  SCATTERED  OBSERVATIONS  TO  A  SCUARE  GRID,  THE  PROGRAM 

C  INVOLVES  A  WEIGHTED  MOVING  AVERAGE  TECHNIQUE  AND  LINEAR  INTERPOL- 

C  ATION  WITHIN  JRID  SQUARES.. 

C 

C  PURDUE  UNIVERSITY  VERSION  2oQ       A*  KEITH  TURNER* 

C  CIVIL  ENGINEERING. 

C  FEBRUARY*  1968* 

C.  USAGE  —  CALL  GRID  {N,X»Y»W3 

C 

C  DESCRIPTION  OF  PARAMETERS  - 

C  N  -  NUMBER  OF  SCATTERED  OBSERVATIONS 

C  X  ~  X  CO-ORDINATES  OF  SCATTERED  OBSERVATIONS 

C  Y  -  Y  CO-ORDINATES  OF  SCATTERED  OBSERVATIONS 

C  w  -  ELEVATION  (Z  CO-ORDINATE)  vALyES  OF  SCATTERED  OBSERVATIONS 

C  DESCRIPTION  OF  CONTROL  CARDS  — 

C 

C  1}  GRID  SPECIFICATION  CARD  - 

C  ROWS  (COLS  1-31   NUMBER  OF  ROWS  IN  GRID  • 

C  COLS  J COLS  4-6)   NUMBER  OF  COLUMNS  IN  GRID  • 

C  GSlZE(COLS  7-10 >  SIZE  OF  GRID  IN  X-Y  INPUT  DATA  UNITS  » 

C  BND  (COLll}   IF  BND=0»GRID  BOUNDARIES  CALCULATED  FROM  DATA. 

C  IF  BND=1*GRID  BOUNDARIES  SPECIFIED  ON  CARD  3  • 

C  LISTK'OLIZi       IF    Ll5T»i»INPUT   IMTA    IS    LISTED. 

C  IF  LIST*2* INTERPOLATED  VALUES  ARE  LISTED. 

C  IF  LIST*3*B0TH  INPUT  AND  INTERPOLATED  VALUES  ARE 

C  LISTED. 

C  IPUNJC0L13!   IF  IPUN*1»GRID  VALUES  ARE  PUNCHED  • 

C  NN'AP(C0L14-15>  NMAP  -  NUMBER  OF  DIFFERENT  MAPS  DESIRED  • 

C  IPL0HC0L16)     IF    IPL01=:UGRID   CALLS    SUBROUTINE    PRTC0L1    • 

C  IF  IPL0T«2»6RID  CALLS  SUBROUTINE  PRTC0L2  • 

C  NOTE—  IF  ANY  ONE  OF    ROWSfCOLS»OR  GSIZE  SPECIFIED tPROGRAM  WILL  COM- 

C  PUTE  THE  OTHERS*  IF  NONE  ARE  SPECIFI ED, PROGRAM  SELECTS  A  GRID  SIZE. 

C 

C  2  5  TITLE  CARD  - 

C  ANY  ALPHANUMERIC  TITLE  IN  COLS  1-72. 

C 

C  3  J  BOUNDARY  SPECIFICATION  CARD  dOPTlOMAL-REQUlRED  ONLY  IF  BND*1  ) 

C  MAXIMUM  tND  MINIMUM  X-Y  COORDINATE  VALUES  i XMAX »XMIN» YMAX *YMIN } 

C  IN  COLS  *-lG,ll~20e21~30»31-40  WITH  DECIMALS  PUNCHED. 

C    REQUIRED  SUBROUTINES  

C 

C  1}  SUBROUTINE  STATS 

C 

C  2)    MAPPING  SURROUTINFS  (  IF  NMAP  GREATER  THAN  Z^RO  ) 

C  IF  IPL0T«1*SUBR0UTINE  PRTC0N1  (PRODUCES  BANDS  OR  LINES) 

C  IF  IPL0T«2»SUBR0UTINE  PRTC0N2  (PRODUCES  SHADED  DENSITY  MAPS) 

C 


r 


SUBROUTINE  GRID  (N,X,Y»W) 


DIMENSION  X(IOOO),  Y(1000}»  W(1000)»  ZUOu»lOO)»  D(20)»  NK20)»  FO 
lRM(lG),  TITL(12)»  A(6)»  B<6)»  C{6) 

REAL  MAXX»MINX»MAXYtMINY»MAXWtMINW»MAX2»y.INZ 

INTEGER  TAPE»YES»9ND»ROWS»COLS 

READ  (5*110)  ROWStCOLS»GSl2E»BND»LIST»IPUN»NMAP»IPLOT 

READ  (5*115)  (TITL(I) ,|«1.12) 

IF  (LISToLT.1qOR*LIST,»GT83}  GO  TO  5 

IF(LIST.NE„2)    WRITE(6*120)     CT!TL< I )» 1=1*12) 
5  IF    (IPUN.EQol)    PUNCH    125*    CTlTMl  )  tl*ltl2 ) 

DO  10  I«ltN 
10     IF  (LlSTeEQel<,ORaLIST*EQ»3)  WRITE(6*130)  I  *X(  I )  »Y(  I  )  *W(  I  ) 

CALL  STATS  ( N  » 1  *  10 JO  » 1 , X  » XBAR ,DUM * XS I G  »MAXX * MINX ) 

CALL  STATS  ( N  ,  1 »1030 ,  1 ,  v » YBAR  »D«JM»  YS  IG»MAXY  ,MINY ) 

CALL  STATS  ( N » 1 • 10UO » 1 * W »WBAR > DUM • WS I G  sMAXW  »MI NW  3 

WRITE  (6,135)  (T1TL(I)*!:=1»12)  ,N 

IF  (BNDoNE®  '.)  GO  TO  15 

WRITE  (6»14J) 

READ  (5*145)  XMAX»XMIN»YMAX.YMIN 

GO  TO  2U 


15     WRITE  (6*150) 

XMAX=MAXX 

XMIN=MINX 

YMAX«MAXY 

YMIN*MINY 
20     FN*MAXW-MINW 

WMAXeFN+MAXW 

WMIN=MINW~FN 

MAXZ«*MAXW 

MINZ"MINW 

WRITE  (6*155)  NMAPtLISTtBND*TAPE 

M=0 

FN=N+0 

DX=XMAX~XMIN 

DY=YMAx-YMIN 

AREA=DX*DY 

BEST=SORT«AREA/FN) 

AVEDISs6EST*l»07346 

WRITE  (6»l60>  XBAR,XSlG*MAXX»MINX»YBAR»YSlG»MAXY»MINY»WBARtWSIG»MA 
lXW*MINW»AREA»AVEDlS»BEST»GS.t2E»R0WS*C0LS 

IF  (GSIZEoLEaOS  CO  TO  25 

COLS*Cle5+DX/GSlZE) 

R0WS=(1«5+DY/GSIZE) 

M=l 

GO  TO  ^5 
C 

25     IF  (COLS.LE»l)  GC  TO  35 
30     GSIZE«DX/FLOATCCCLS-l) 


R0WS=<1»5+DY/GSIZE) 
M»l 

GO  TO  45 
C 

35  IF    (ROWSoLE.l)    GO    TO    4C 

GSIZE=DY/FL0AT(ROWS~I) 
COLS=ao5+DX/GSlZE) 
M=l 

GO    TO    45 
C 
4C  GSIZES8£ST 

COLS~(lo5+DX/GSlZE) 
ROVi/S=U.5+DY/GSlZE) 
45  CONTINUE 

IF    {COLS. LE* 100)    GO    TO    50 
C0LS=100 
GO    TO    30 
C 
50  WRITE    f6»l65) 

IF    (M«EQ»1)    WRITE    <6»170) 

IF    (M.NE.lJ    WRITE    { 6*175 i 

WRITE    (6,180)    XMAX*XMIM     ■:  MAX  »YMIN,GS  IZE»R0v,'S  »COLS 

IF    (LXST.EQ.2«OR.LISToEG.3?    WRITE(6»185)     { TI TL( I) » I =1 , 12 ) 

DY= ( DY-( ROWJ-1 . ) *6S IZF) /2.0 

0X=(DX-(C0LS-1. >*GSlZE)/2„0 

YMAX=YMAX-DY 

XM!N=XMIN+DX 

HGRID=GSIZE/25.0 

F  N= { WMAX-MAXW ) / 2  e  C 

MAXX-MAXW+FN 

MlNX*MINW-FN 

00  o~  MsltROWS 

YPSYMA.X~{FL0AT(M-1)  )*GSIZE 
DO  90  J«l»COLS 

XP=XMIN+(FLOAT(J-l))*GSIZE 
DO  55  1=1,9 
55  D(I)=1.E10 

Nlil)*l 
Nl(2)=l 
00  75  I=1,N 

DIST=< (XP-XC I) )**2)+< (YP-YU 5 5**2) 
IF  (DISToGE.HGRID)  GO  TO  60 
Z(M,J)=W(I) 
GO  TO  85 
C 

60  YES»MINO(I»e) 

DO  70  K*1»YES 

IF  (ftIST«GF»D(IC))  GO  TO  70 
DO  65  JJ»KtYES 
KK=YES+K-JJ 
MM=KK+1 
D(MM)=D<K.K) 
65  N1(M.M)*N1{KK> 


D(K)*DIST 

GO  TO  75 
C 

70  CONTINUE 

7  5  CONTINUE 

DUMlaO.O 
PSUM=0.0 
DO  80  1=1.6 
YES=N1(I) 

MAXW«1«0/(SQRT«D(I)) S 
PSUM=PSUM+MAXW 
80  DUM1=DUM1+W(YES)*MAXW 

DUM6=DUM1/PSUM 
Nll=Nl(l) 
DUM2=W(N11} 
DUM*  ( DUM2+DUM6  )/2a0 

IF  (DUM.GT«MAXX»CR.DuM.LToMINX)  Z(M»U)«DUM2 
IF  (DUM.LE.MAXX.AND.DUM.GE-MINX)  ZtM»J}=DUM 
85  CONTINUE 

90        CONTINUE 

IF  CLIST«EOe2«OR*LIST-EQ»3)  WRITE(6»190)  C Z (M» I } » 1*3 »C0LS ) 
IF  (IPUN.EQ.l)  PUNCH  195*  ( Z«M* I ) » 1*1 »COLS) 
95     CONTINUE 

NEWT0T»R0WS*C0LS 

CALL  STATS  <ROWStCOLS»100»100, ZtZ9AR»XP»ZSIGoZMAX«ZMIN) 
YMIN~YP 

WRITE  { 6  s20 J )  NEWTOT  » XM1 N *  YM I N  »ZMAX »ZMI N » ZBAR ,ZS I G 
WRITE  (6,205) 
IF  CNMAP«LT«1)  GO  TO  1C5 
DO  100  JNUM«1,NMAP 

IF  JIPLCToEQ.l)  CALL  PRTCON1 (Z»ROWS»COLS »TITL1 
IF  (IPLOT.EQ.2J  CALL  PRTCCN2(Z»R0WS»C0LS»TITL> 
100    CONTINUE 
105    CONTINUE 

RETURN 
C 

110    FORMAT  (213, F4o0, 311,12*11) 
115    FORMAT  (12A6) 

120    FORMAT  ( 1H1  .  12A6,//,3Xt24H  COORDINATE  OBSERVATIONS) 
125    FORMAT  (6X»12A6) 
130    FORMAT  (5X*4H  I  =»I3»7H  X(I>  =»Fl0.3»7H  Yd)  «»F10»3»7H  w(I)  =  »Fl0 

1*3) 
135    FORMAT  { 1H1 » 12A6*// »3X *  I  5 » 13H  OBSERVATIONS,// ) 
140    FORMA'!  (3X,25H  EOUNOARIFS  PREDETERMINED) 
145    FORMAT  (4F10.0) 

150  FORMAT  (3X,32H  BOUNDARIES  CALCULATED  FROM  DATA) 
155  FORMAT  (3X»6HNMAP  -,I3»7H  LIST  -»I3*6H  BND  =*I2,7H  TAPE  =  »!'3) 
160  FORMAT  { 3X» 15H05SERVED  VALUES, //»3X»6HX3AR  = ,F10«3»6HXSIG  =sFl0°3, 
16HMAXX  a»F10,3»6HMINX  «»F1Q.3»//3X»6KY3AR  - ,FlO*>3 *6HYS 13  =,F10»3»6 
2HMAXY  a»Fl0.3»6HMINY  = »F10.3 »//3X»6HWBAR  = ,F10. 3 »6HWSIG  =,F10.3,6H 
3MAXW  -,Flu,3,6HMINW  »,F1 Je3»//3X»6HAR£A  ■ ,F1C 3»8HAVEDI S  =»F10.3»6 
4H8EST  =»F1C.3,//3X,7HGSIZE  «»F10«3t6HROWS  =»I4»6HC0LS  --*I4,/J 


175         FCT>"A7    (*X??1J    CiLCULAjFTi    PRQ.V    r>ATAe/) 

180    FORMAT  <3Xt,16H  VAL"Ee  7?-  yw*y  =«Fln.3 ,3**7^  v**I»;  r.:.p 

U0e3j.3X»7H    YMAX    =»F1^»3  ,»?vr>8"    Gc.!.r-~    -••*c  1 

-*I4*8X*7H   COLS    = >I<+8 
185         FORMAT    < 1H1»1?A6»//:  "  *f>) 

190         FORMAT    ( 1H0 * C 10F1O.3)  ) 
195         FORMAT    58F10»3) 

FORMAT    (1H0»3X»3HN    =  9I3*6HX.M]  6HZ 

1MA>     =>FlU»3*3X>6H  =  »F10»3s3X*6HZStG    '■ 

205         FORMAT    < 13H1    £ND    OF    DA] 


C  PURPOSE  - 

C  TO  CALCULATE  STATISTICAL  MEASURES  OF  ONE  OR  1  IO.NAL  DATA 

C  ARRAYS*  ME  .      DETERMINED  ,N,2 JVARI ANCE»3 )STANDARD  D 

C  .  LUF»5)M|N  ;'■'!!!-  VALUF0 

C         U 

C  L    STAT'  COLSvIUJJ,  R»XSIG»>  MIN) 

C 

C  /MBcR    OF 

C  -    MUW3ER    C'r    COL 

C  II 

c  J j  •  y/  nu^pp 

C  X  -  DATA  ARRAY  BEING 

C  XBAR  -  MEAN  OF  AR? 

C  XVAR  -  VARIANCE  OF  A9RAY. 

C  XSIG  -  STANDARD  DEV 

C  X.MAX  -  MAXIMUM  VALUE  OF 

C  XMIM  - 

SURI 

C 

C      PRC        3  FIND  .  FS 

C 

COLS 
.  J ) 
fS*COLS 
XSl 

1) 
! { 1  *  1 J 
DO    5     I«1»RC 
5    J=l»COLS 

[»Ji.nT,:  =x{ i* J) 

»J).LT.XMI  ,  J) 

S>  J) 

9  SQR+XU 

iL 
XV'  *-2«.0*XI 

■QRKXV 


PART  IV 

PRODUCTION  OF  CONTOUR  MMPS 
ON  1SS  PRBKSSR 


contents  of  ass  section 

l)  Subroutine  FKECOB  1 

2}  Subroutine  FRTCON  2 

3)  Subroutine  EXPSO 

•    h)  Programs  MAPPER 

5)  SAMPH5  am 


ff 


^SUBROUTINES  PRTCONX  AND  PRTCON2  -  GENERAL  INSTRUCTIONS 

THESE  PROGRAMS  ARF  BASED  ON  A  MAD  LANGUAGF  PROGRAM  OBTATNFn 
FROM  PROFESSOR  W.  R.  TOBLER,  DEPT.  OF  GEOGRAPHY,  UNIVERSITY  OF 
MICHIGAN,  ANN  ARBOR.  BOTH  SUBROUTINES  PRODUCE  CONTOUR  MAPS  OF  GRIDD- 

cDuaSJc  °VHE  PRIwter  usinc  Unear   interpolation  WITHIN  THE  gr?d 

tSpp5f;  H0«EyER  PRTCON1  ASSIGNS  A  SINGLE  CHARACTER  TO  EACH  CONTOUR 
luIER/Ah  < /LLOrflNG  A  MAXIMUM  OF  19  INTERVALS)  WHILE  PRTCON2  USES 
THE  OVERPRINT  OPTION  TO  PRODUCED  PATTERNS  FOR  EACH  INTERVAL  THAT 
ARE  LIGHT  FOR  LOW  AREAS  AND  DARK  FOR  HIGH  AREAS  {  A  MAXIMUM  OF  10 
CONTOUR  INTERVALS  IS  POSSIBLE).  PRTC0N1  HAS  THREE  OPTIONS  -  CONTOUR 
LINES  MAY  BE  PRINTED,  ALTERNATE  CONTOUR  BANDS  MAY  BE  PRINTED,  OR Ml 
?2JN^d2S-I,II°NS  °M  THE  MAP  MAY  BE  PANTED.  PRTCSN2  HAS  T  0  OPUO^ 
JSnTO^R  INTERVALS? fNTE°  WITH  °R  WITH°UT  WHITE  B°RDERS  SEPARATING  THE 

rprr  dSImt^R?UTINES  C*H    CCNTOuR  IRREGULARLY  SHAPED  AREAS  IF  THE 
?2I?t£?I?.T&  FALLING  OUTSIDE  THE  AREA  OF  INTEREST  ARE  GIVEN 

R     !  VERY  HIGH  OR  LOW  VALUES  (  VALUES  MUST  BE  ONE  DATA  ELEV- 
ATION RANGt  BELOW  LOW  OR  ABOVE  HIGH)  <s 

THE  SUBROUTINES  PREPARE  A  FREQUENCY  TABLE  SHOWING  THE  PERCFNT 

/mr!^,^"1  AREA  °F  INTERE5T  F0R  IRREGULAR  AREAS,"  N  MER 
ICAL  INTEGRATION.  THE  VALUES  MAY  BE  FURTHER  ANALYZED  BY  SUBROUTTNF 
HYPSO  WHICH  PERFORMS  HYPSOMETRIC  (  AREA-ALTITUDE)  ANALYSIs! 


SPECIAL  NOTE  ~ 

THE  MAXIMUM  MAP  WIDTH  IS  12.7  INCHES. THE  WIDTH  OF  THE  PAPER   rHF  mavtmiim  m..u 

IF  PS-(C0LUc2f,?^^rf  °ATA  IS  A  ^lON^THE  IpECIFIED^M.p'wiD  hT  *"*" 

A   ^     in  ;1!   (  SEE  D^I«TIONS  BELOW  >  IS  LESS  THAN  2,  FEWER 

MMCA-?c£?IrIn'g   P0SITI°NS  WOULD  BE  ASSIGNED  TO  EACH  COLUMN  SPACE  AND  SO  AN 
UNSATISFACTORY  MAP  WOULD  RESULT. THUS  THE  PROGRAM  WILL  NOT  PLOT  A  MAP   FP   i« 
LESS  THAN  2    ,  AND  ACCORDINGLY  THE  MAXIMUM  NUMBER  OF  COLUMNS  FOR A  MAP  OP 
MAXIMUM  WIDTH  IS  42.  IF  MORE  COLUMNS  ARE  REQUIRED  TH*  MAP  MUST    BE  MADE  IN 
t^UL?  LBth"S'v;cCT-Y-PARTd°F  THE  °ATA  IN  EACH  PASS.  THE  NUMBER  COLUMNS 
!LTBE  TH~  SAME  IN  EACH  PASS  FOR  THE  MAPS  TO  BE  AT  THE  SAME  SCALE.  TO  DO 
THIS  IT  MAY  BE  NECESSARY.  AND  IS  ADVISABLE,  TO  REPEAT  SOME  COLUMNS. 


C  SUBROUTINE  PPTCONl 

C  PURPOSE  - 

C  TO  PRODUCE  CONTOUR  MAPS  ON  THE  PRINTER. 

C 

c    — „__. . _ __ —_-__- .„—._„ 

C  USAGE  - 

C  CALL  PRTC0N1(Z.R0'/S?C0LS,TITL) 

C  Z       ARRAY  OF  GRIDOEO  ELEVATION  VALUES  (MAXIMUM  ~  100  X  100  ). 

C  ROWS  -  NUMBER  OF  ROWS  IN  Z  ARRAY 

C  COLS  -  NUMBER  OF  COLUMNS  IN  Z  ARRAY. 

C  TITL  -  ALPHANUMERIC  TITLE  1-72  CHARACTERS. 

C 

C  PUOUE  UNIVERSITY  VERSION  2.0         A.  KEITH  TURNER, 

C  CIVIL  ENGINEERING, 

C  FEBRUARY  1968. 

C  CONTROL  CARD  ORDER  - 
C 

C  1)  MAST  :R  CONTROL  CARD  - 

C  C0,4  (COLS  1-2)  THE  NUM6ER  OF  CONTOUR  INTERVALS  (MAXIMUM"  19 )  • 

C  (  CON  =  DATA  RANGE  /  DESIRED  CONTOUR  INTERVAL). 

C  TOUR  (COL  3)    IF  TOUR=0, PROGRAM  USES  MAXIMUM  AND  MINIMUM  FROM 

C  DATA  FOR  UPPER  AND  LOWER  CONTOUR  LIMITS. 

C  IF  T0UR*1»  DESIRED  TOP  AMD  BOTTOM  CONTOURS  READ 

C  IN  FROM  SPECIAL  PROGRAM  CARD  (SEE  INFLOW). 

C  IF  T0UR«2.>  VARIABLE  CONTOURS  READ  IN  (LOW  TO 

C  HIGH).  IN  THIS  CASE  SET  CON  EQUAL  TO  THE  NUMBER 

C  OF  VALUES  TO  BE  READ  IN  . 

C  IF  T0UR»3»  CONTOUR  INTERVAL  BECOMES  1/2  STANDARD 

C  DEVIATION  ?  MAX=+3  SIGMA,  Ml N=— 3  SIGMA), 

C  LINES  (COL  4)   IF  LtNES=0»  CONTOUR  LINES  ARE  PRINTED. 

C  IF  LINES*!  '-ALTERNATE  (EVEN  VALUED)  CONTOUR  BANDS 

C  ARE  PRINTED,, 

C  IF  LINES=2->ALL  BANDS  PRlNTEDtALL  PRINT  PCSNSL 
INCHESCCOLS5-8) INCHES  IS  MAP  WIDTH  IN  INCHES.  IF  GREATER  THAN 

C  12*7  OR  BLANK, SET  TO  10.0 

C  (COL  =  MAP  WIDTH  In  CHARACTERS* I NCHES*1( 

C  IHYPSO  (COL10)  IF  IHYPS0*lSU8R0UTINE  HYPSO  CALLED. 
C 

C  2)  SPECIAL  CARD'S)  - 

C  REQUIRED  ONLY       rOUR  ~  1  OR  2 

C  IF  TOUR* 1,  THIS  CARD  SPECIFIES  MINIMUM  AND  MAXIMUM  ELEVATION 

C  VALUES.  —  ZMIN  (COLS  1-10}  AND  ZMAX (COLS1 1~20 ) . 

C  IF  T0UR=2,THIS  CARD  SPECIFIES  A  FORMAT  (INCON)  TO  BE  USED  IN 

C  READING  IN  A  SERIFS  OF  IRREGULARLY  SPACED  CONTOUR  LINES.  THE 

C  VALUES  ARE  ARRANGED  FROM  LOWEST  TO  HIGHEST  ON  SUBSEQUENT  CARDS* 
C 


SUBROUT I NE  PR  TCON1 ( Z , ROWS • COLS » T I  T  L  ) 

DIMFNSION  FORM (10) »TITL(12) ,SYMTAB(21) »PSYM( 131 ) »Z( 100*100 J » 
1  iMCON(lu) ,DRANG<25> ,TA3E<25) »ZPREV(131) »SYM(13l) 

DIMENSION  CUMH(27) 

INTEGER    ROWS, COLS, CON, TOUR »C0LMAX, COL, PSYM* PS, GRID, TOUT, TLOW, 
1    SYM,  1LANK,PLUS»MINUS»STAR,SYMTAB 

INTEGcR    XSTOP,YST0P 

REAL    MINZ»MAXZ, INCHES, INCON 

DATA    CSYMTABC  I)  »  1  =  1  *21  )/lHO,lHl  ,1H2*  1H3  ,  IWi  ■•  1HC,  ,  1H6  »1H7»  1H8  »  1H9  , 

1  lHA»lH3»lHC*lHD»lHF*lHF»lHGtlHH»lHI »1HJ»1HK/»FLANK/1H  /9PLUS/IH+/ 

2  *MINUS/1H~/,STAR/1H*/ 

I    READ  (5»loOO)  CON » TOUR »L I NES» INCHES  9  I HYP SO 
lOOu   FORMAT ( I  2  * 2  1 1 , F40  C , 2X ♦ i 1 > 

C       **********************#***************** 

C    ZERO  COUNTERS  AND  ARRAYS  USED  IN  FREQUENCY  TABLE 

C       *****************  ***#*****#***-!(•********* 

JJ»0 

MM«0 

SUM»0«0 

TLOW=0 

TOUT=0 

DO  5  1=1,25 
5      TA3ECI)=0.0 

N=R0WS*C0LS 

C      N.  *  NUMBER  OF  DATA  POINTS 

C     *******#********************»*******#*** 

WR I  TEC 6 » 1030)  TITL,N,CCLS,R0WS 

1030   F0RMAT(1H1?12A6*//,22H  MAP  IS  DEVELOPED  FROM, I = *30H  GRIDDED  VALUES 
I, CONSISTING  OF  ,I3»12H  COLUMNS  BY  *I3,6H  ROWS*//} 
CALL  STATS { ROWS tCOLS* V  0 , 100 ,Z ,Z3AR ,XP,ZS IG*MAXZ *MI H?. ) 
C  **************************************** 

C    DETERMINE  CONTOUR  LINE  ELEVATIONS 

C    TOUR  VALUE  DETERMINES  METHOD  USED  IN  DEFINING  CONTOURS 
C     **  *******  *********  ****  *****  ****  ********* 

IF  CTOUR  .GT.  3  .OR*  TOUR  ,L£-  05  GO  TO  60 

GO  TO  (30, 40, 50), TOUR 
30     WRITE  (6*1060) 
1060   FORMAT  C25H  CONTOUR  LEVELS  SPECIFIED} 

READ  (5*1070)  ZMIN»ZMAX 
1070   FORMAT  (2F10.0) 

GO  TO  70 
40     READ  (5*1020)  INCON 
1020   FORMAT  C10A6) 

WRITE  (6,1080) 
1080   FORMAT  ( 37H  VARIABLE  CONTOUR  INTERVALS  SPECIFIED) 

READ  C5» INCON)  i DRANG ( I ) , I =1 »CON> 

ZMIN=DRANGC1) 

ZMAX=DRANG(COM) 

GC  TO  70 
50     ZMJN»ZRAR-3.0*ZSIG 

ZMAX*ZBAR+3.0*ZSIG 

C0»M«13 


WRITE"  {6,1090) 
1090   FORMAT  { 30H  CONTOUR  INTERVAL  ~  Ul/2  SIGMA) 

GO  TO  7^ 
60     WRITE  (6.1100) 
1100   FORMAT  (26H  CONTOUR  LEVELS  CALCULATED) 

ZMIN*MINZ 

ZMAX=MAXZ 

7C     WRITE(6»1110)  R0WS*C0LS,Z6AR,ZSIG»MAXZ*MIN? 

Ul-   FORMAT  (13H  GRIDDED  DATA/»42H  LINEAR  INTERPOLATION  WITHIN  GRID  SOU 
1ARES  /,17H  OBSERVED  VALUES  /»   8H  ROWS  -  *I3*I0H    COLS  =  , 
213, /,8H  ZbAR  =  ,F10.3,10H    ZSIG  -  »Fl0o3»l0H    MAXZ  -  ,FlC«3»loH 
3   MINZ  n  ,F10.3) 

C    CALCULATE  MAP  WIDTH  IN  CHARACTERS  (COL) 
C     ♦****•■*****#*#«•##•**•*■&#■*#*  *#**«•***** "»#«*** 

IF  ( INCHES»LE.0»0)  COL«100 

IF  (INCHES  .GTe  12.7)  COL*127 

IF  (INCHES.GT.O.O.AND.  I NCHFS-LE.12.7 )  COLs T NCHES*l0«0 

PS= ( COL-COLS ) / ( COLS-1 ) 
C     #***•****#*#•*##*#  *##■»•#■*##  »•#«#«■#*#«**#*#*♦ 

C    TEST  GRID  SIZE— IF  GRID  SMALLER  THAN  2  CHARACTERS  BETWEEN  COLUMNS, 
C    AN  ERROR  STATEMENT  IS  PRINTED  AND  PLOT  ABORTED 
C    **#*•»#«•***#####♦#-*•«■*■»*##■»*•»##*«•*#*#*«■**«• 

IF  (PS*GT«.l)  GO  TO  80 

WRITE  (6,1120) 
1120   FORMAT  ( 25H  TOO  MANY  COLUMNS  IN  GRID) 

GO  TO  999 
C 

C    CALCULATE  REQUIRED  CONSTANTS 
C 
80     COLMA/sP S*{ COLS-1 J+COLS 

YINC=U-'/C  {  FLOAT  (PS+1)  )*0#6) 

GRID*PS+1 

XINC»1„0/FLOA7(GRIO) 

RANGZ**ZMAX-ZMIN 

RHIGH=ZMAX+RANGZ 

RLOW*ZMlN-RANGZ 

IF  (C0N-EQ.0.0R.C0N.GT.19)  C0N--5 

FC0N=C0N 

CI-RANGZ/FCON 

IF  (T0uR«EQo2)  GO  TO  ICO 

JJJaCOM+1 

DO  90  T=1»JJJ 
90     DRANG(I)=ZMIN+FL0AT(I-}.  )*CI 
100    CONTINUE 

WRITE  (6»1130) 
1130   FORMAT  (1H  /,16H  VALUES  EMPLOYED) 

WRITE  (6,1140)  ZM I N s»ZMAX »CI ,GR ID»X I NC , Y I NC » COM  » PS »  COL  MAX 
1140   FORMAT  (  8H  ZMIN  ■  »F10.3»10H    ZMAX  =  »F"       I    CI  =  »Fl0.3/t8H 
1GRID  ■  »I5,10H    XINC  ■  •F10.3»10H    YINC  -  ,Fl0.3,/f7H  CON  =  15, 
2  8H    PS  =  »I5»12H    CCLMAX  *  »?5) 


C    WRITE  OUT  CONTOUR  SYMBOL  TABLE 

C     ***»«■**«***###•»■»■*#*#•»»*****»*#**##***#«■* 

WRITF  (6.1150) 
1150   FORMAT  (35H4CONTOUR  SYMBOL  TABLE  IS  AS  FO-LOWS) 

IF  (LINES. EO.l)  WRITF  (6,1160) 
1160   FORMAT  (43HC(0NLY  EVEN'  SYMBOLS  PRINTED  IN  BODY  OF  MAP)  ) 

WRITE  (6.1170) 
1170   FORMAT  <7H0SYMF}0L»11X»5HVALUE) 

IF  (LINES.NE.O)  GO  TO  110 

WRITE  (6.1180) 
1180   FORMAT  (5 HO    *»6X»7HN0  DATA»/»5H     -,4       JELOW   LOW I 

WRITE  (6.1190)  <SYMTAB(I>»DRANG(I )»I*1»CQN) 
1190   FORMAT  ( AX  ,Al »2X ,FI2e4 ) 

WRITE  (6,1200) 
120u   FORMA  ('  <4XtlH+.4X»10HA30VE.  HIGH) 

GO  TO  I2u 
110    WRITE  (6,1210) 
121o   F0RMAT(5H0    **25X*7HN0  DATA./.5H     -.23.X.10HBELOW   LOW) 

WRITE  (6.1220)  (SYMTABd  ).DRANG(I  )  »DRANG(  1+1)  .I~lt.COM) 
1220   FORMAT  ( 4X • Al » 2X ,F12»4 ,6H   TO   .F12e4? 

WRITE  (6*1230)  DRANG(C0N+1) 
1230   FORMAT  (4X»1H+»2X»F12.4»10H   OR   OVER: 
120    WRITE(6.1030)  TITL.N.CCLS.ROWS 

XSTOP  «  GRID 

GRID-COLS-1 

C0L  =  3 

DO  140  J=ltGRIO 

DO  130  I  =  2 .XSTOP 

SYMI COL) "BLANK 
130    COLsCOL+1 

SYM( COL) =PLUS 
140    COL-COL+1 

K-COLMAX+2 

SYM(2)=PLUS 

SYMID-BLANK 

SYM(N) -BLANK 

WRITE  (6,1240!  (SYM( I) ,  Is! »N) 
1240   FORMAT  (1H9.129A1) 
C 

C    ESTIMATE  VALUE  OF  ALL  PRINTING  POSITIONS  ON  MAP 
C 

C      LOOP  THROUGH  MAP  LENGTH*  ROW  BY  ROW 
C 

ITOP-1 

I  ROW  =  ROWS  +  1 

DO  340  MM  =  2, 1  ROW 

YSTOP  =  l./YINC  +  .5 

IF  CMM.EO.IROW)  YSTOP-1 

M-MM-1 


'.-; ...  -  ■ 


145         ypssO.J 

PSYM(l)    ■    PLUS 

PSYM(N)    ■    PLUS 

DO    34U    KJ=l»YSTOP 
150         CQL*1 

Zl    =    YP*(Z{M,v;,l)-z(r/Ul))-5-ZCM,l) 
C 

C      LOOP  THKOUGH  MAP  WIDTH »C0LUMN  BY  COLUMN 
C 
190    II*C0LS-1 

JLFT»2 

DO  330  J=13II 

JJ=J+1 

Z2  =  YP* «  Z  « MM  t J J ) -Z ( M » J J ) ) +Z ( M t J J 3 

ZDELT=Z2-Z1 

XP~0o0 

DO  480  KK  =  1»XST0P 

JLFT*JLFT-1 
195    C0L=CQL+1 

200  ZEST=XP*ZDELT+Z2 
IF{ZEST»GE.ZMIN«AND.ZEST«LE*ZMAX)  GO  TO  325 
IFCZEST«GT<,ZMAX)  GO  TO  203 

IF  (ZEST.LE»RLOW)  GO  TO  206  . 

IF(  ITOP«EO»l»ORe?/M.EQ8  I  ROWoORi  JLFT«EO*l )  GO  TO201 

PSYM(COL)*BLANK 

GO  TO  2o2 

201  PSYMC COD "MINUS 

202  TLOW=TLOW+l 
ZPREV(COL)=ZEST 
ZLAST'ZEST 

GO  TO  124 

203  IF(ZEST.GE.RHIGH)  GO  TO  2?6 
IF(ITOP.EQ,1.0R.MM,EO#IROW,OR«JLFTaEQ.l)  GO  T02C4 
PSYM(COL)=BLANK 

GO  TO  2U5 

204  PSYM(COL»=PLUS 

205  TABE(C0N+U«TABE(C0N+1)+1.0 
ZPREV(COL)=ZEST 

ZLAST  'ZEST 
GO    TO    124 

206  IF{ ITOP.EO.l.OR.MM.EO. I ROW.OR. JLFT.EQ.l )    GO    TO207 
PSYM«COL)*BLANK 

GO    TO    20  8 

207  PSYM( COL) "STAR 

208  T0UT=T0UT+1 
ZPREV(COL)=ZEST 
ZLAST=ZEST 

GO    TO    480 
325         IM(ZEST-ZMIN)/CI)  +  1. 
TABECI)  =  TABE(  D+1.0 
LVAL«LIN£S+1 


I F  (  I  TOP-  EQ.  1  .OR.  KM.  EOo  I  ROW  .  OR.  Ji_FT  e  EQ  «,  1  ) LVAL= 
GO    TO    (326,327*328  5     »LVAL 

326  DO    309    1=1, CON 
ZSTEP«DRANG( I) 
SW1=1. 
SW2=1. 
SW3=1. 
SW4=1» 

IF     (ZSTEP.GF.ZLAST.OR.ZSTEP.GF.ZPRFVtCOL) )     SWl*0. 

IF    (ZSTEP.LE.ZESTJ    SW2»0c 

IF    (ZSTEP.LT.ZLAST.OR.ZSTEP.LT.ZPREV(COL))    SW3*0. 

IF    (ZSTEPoGEoZEST)    SW4=0» 

IF(  (SWloEU»0o.AND.SW2.EG.C*)«0R.  {  SW3.E0.CU  .AND.SW4.EG.0o  )  )G0    T0310 

309  CONTINUE 

PSYM( COL) "BLANK 
GO    TO    320 

310  PS  YM  { COL)=S  Y.^  TABU) 
320        ZLAST=ZEST 

ZPREV(COD=ZEST 
GO    TO    124 

327  K*I/2 

IF((K*2).NE.I)  PSYM(COL)~SYMTAB(I ) 
IF{  (K*2>  »EQ.I)  PSYM< COL)  "BLANK 
GO  TC  124 

328  PSY^(COL)»SYMTAB(I) 
ZPREV(COD=ZEST 

124    SUM=SUM+1.0 
480  XP=XP+XlNC 

Zl  =  Z2 
330    JLFT=0 

COL  =  COL  +  1 

IF{Z(M»COLS).GE«ZMIN.AND»Z(M»COLS)«LEoZMAX)  GO  TC  333 

IF(Z(M.COLSUGT«ZMAX)  60  TO  335 

IF(Z(M„CQL5).LE.RL0W)  GO  TO  336 

PSYM(C0L)=MINUS 

TL0W*TL0W+1 

GO  TO  334 

335  IF(Z(M.COLS)«GE.RHIGH)    GO    TO    336 
PSYM(COD'PLUS 
TAB£(CON+1)=TABE(CON+1)+I.O 

GO    TO   334 

336  PSYM(-COL)*STAR 
T0UT«T0UT+1 

GO  TO  332 

333  I  =  ((Z{M*C0LS)~ZMIN)/CI)+1« 
TABEM)  =  TA3E(  Ij+1.0 

PSYM  ( COD  =SYM  TABS  I  ) 

334  SUM-SUM+1.0 

332    WRITE  (6,1240)  (PSYM(L) tL*l»N) 
YP=YP+YlNC 

PSYM(l)  *  BLANK 

PSYM(N)=BLANK 
340    ITOP=0 

WRITE  (6,1240)  CSYM ( L) »L=1 PN) 


C    CALCULATE  AND  PRINT  FREQUENCY  TABLE 

C     *#****##-*#*###*##**-»**#*#**##*###■*:- •»**•#** 

WRITE  (6»1250>  TOUT 
1250   FORMAT  ( 16H1FREQUENCY  TABLE»/«12X»8HINTERVAL  »16X»5HUNITS  *6X*l0HP 
1ERCENTAGE  ,4X. 10HCUMULATIVE  »/»16X»7HNQ  DATA  ,UX»I8) 

PTAB= ( FLOAT (TLOW) /SUM )*1 00,0 

CTAB=*TA6 

WRITE  (6*1260)  TLOW»PTAS»CTAB 
1260   FORMAT  (14X»9H3EL0W  LOV*»  11X*  I8»2  (3X»F12«4  >  ) 

DO  350  I«lfCON 

VALUE=DRANG« I) 

VAL1=DRANG(I+1> 

VALl«=DRANG(I  +  l) 

PTAB=(TABE( I  J /SUM ) *100.0 

CTAB=CTAB+PTAB 
350    WRITE  (6V1270)  VALUE9VAL1 »TABE( I ) •PTABtC 
1270   FORMAT  <1X*F12.4»6H   TO   ,F12.493X»F8.0,>2  <  3XtF12*4)J 

I=C0N+1 

PTA3*(TABE< I )/SUM)*100*0 

CTAB=PTAb+CTAB 

WRITE  (6»1280)  VAL1»TABE( I )»PTA6»CTAB»SUMfCTAB»CTAB 
1280   FORMAT  ( 1X»F12.4»10H   OR   OVER  »11X*F8.0 »2( 3XtF12.4 ) »/»18X»5HTOTAL 

1*11X»F8<»0S2(3X*F12.4)  ) 
C     *#*H$  *#*#***•*###***#*  »#*#*****■»***** 
C    IF  IHYPSO  =  1  »  PERFORM  HYPSOMETRIC  ANALYSIS 
C     *********  #########*-s-##tf#-**#  #*■»***#* 

IF  (IHYPSO)  450*999*450 
450    CALL  HYPSO(0RAN6,2.MAXtZ.viIN,CONiCUMHtTITL) 
999    CONTINUE 

RETURN 

END 


C  SUBROUTINE  PRTCON2 

C  PURPOSF  - 

C  TO  PRODUCE  GRADED  DENSITY  CONTOUR  MAPS  ON  THE  PRINTER 

C 

C  USAGE  - 

C  CALL  PRTCON2<Z<:ROVS9COLS,TITL> 

C  Z       ARRAY  OF  GRIDDED  ELEVATION  VALUES  J MAXIMUM  =  100  X  100  J* 

C  ROWS  -  NUMBER  OF  ROWS  IN  Z  ARR  8 

C  COLS  -  NUMBER  OF  COLUMNS  IN  Z  ARRAY* 

C  TITL  -  ALPHANUMERIC  TITLE  1-72  CHARACTERS. 

C  PUDUE  UNIVERSITY  VERSION  1.0                  TURNER* 

C  CIVIL  ENGINEERS 

C  -UARY  1968. 

C  CONTROL  CARD  ORDER  - 

C 

C  1)  MASTER  CONTROL  CARD  ~ 

C  CON  (COLS  1-23  THE  NUMBER  OF  CONTOUR  INTERVALS  {MAXlMUM-lO  >  • 

C  {  CON  =  DATA  RANGE  /  DESIRED  CONTOUR  INTERVAL) • 

C  TOUR  (COL  3)    IF  T0UR=0*PR0GRAM  USES  MAXIMUM  AND  MINIMUM  FROM 

C  DATA  FOR  UPPER  AND  LOWER  CONTOUR  LIMITS* 

IF  TOUR-U  DESIRED  TOP  AND  BOTTOM  CONTOURS  READ 

C  IN  FROM  SPECIAL  PROGRAM  CARD  {SEE  BELOW). 

C  IF  TOUR-2.  VARIABLE  CONTOURS  READ  IN  (LOW  TO 

C  HIGH).  IN  THIS  CASE  SET  CON  EQUAL  TO  THE  NUMBER 

C  OF  VALUES  TO  BE  READ  IN  • 

C  LINES  (COL  4)   IF  L \ NES*0»BLANK  BORDERS  ARE  PLACED  BETWEEN 

C  CONTOlR  INTERVALS* 

C  "  LINES»1»BLANK  BORDERS  NOT  DEVELOPED. 

INCHESCCOLS5-8) INCHES  IS    MAP  WIDTH  IN  INCHES.  IF  GREATER  THAN 

C  12*7  OR  BLANK* SET  TO  10 oO 

C  (COL  *MAP  WIDTH  IN  CHARACTERS*INCHES*10»0» 

C  IHYPSO  CCOL10)  IF  IHYPSO=lSUBROUTTNE  HYPSO  CALLED. 
C 

C  2)  SPECIAL  CARD(S)  - 

C  REQUIRED  ONLY  WHEN  TOUR  -    1  OR  2 

C  IF  TOUR=l9  THIS  CARD  SPECIFIES  MINIMUM  AND  MAXIMUM  ELEVATION 

C  VALUES.  --  ZMIN  (COLS  1-10)  AND  ZMAX (C0LS11-20) • 

C  IF  T0UR=2»THIS  CARD  SPECIFIES  A  FORMAT  fINCON)  TO  BE  USED  IN 

C  READING  IN  A  SERIFS  OF  IRREGULARLY  SPACED  CONTOUR  LINES.  THE 

C  VA.UES  ARE  ARRANGED  FROM  LOWEST  TO  HIGI       M  SUBSEQUENT  CARDS. 
C 


SUBR0UTlNE(PRTC0N2<Z*R0WS,C0LS,TlTL) 

DIMENSION    FORMUOsTITU  12)  ,Z{  100,100  )  .DRANG <  25  )  »SYMl  131 ) 

DIMENSION    SYMTAB^10,2),PSYM«131,2) 

DIMENSION  ZPREV{131>,INCON(10),TABE<25)»CUMH<27) 

INTEGER  ROWS. COLS, CON, COLMAX »COL*PSYM,PS »GR I D,SYM,3LANK, PLUS t 
1SYMTAB,XST0P,YST0P 

INTPPFP  TOUR t STAR 

INTEGER  MI NUS»D0T50WF , EYE, LPAR,RPAR, SLASH, EGGS, EQUAL 

REAL  INCHES»MINZ,MAXZ,INC0N 

DATA  BLANK/ 1H  / ,PLUS/1H+/ .MINUS/1H-/ ,D0T/1H*/,0WF/1H0/,EYE/1HT/ * 
1EGGS/1HX/»^0UAL/1H=/, SLASH/1 H//,LPAR/lH(/«RPAR/lHV/*STAR/lH*/ 

SYMTAB(ls>l)  =  DOT 

SYMTAB(1»2)=BLANK 

SYMTAB(2»1)»LPAR 

SYMTAB(2,2)*PLANK 

SYMTAB(3»l)=OWE 

SYMTAB{3»2)»OWE 

SYMTAB(4,1)=EGGS 

SYMTAB(4»2)«EGG5 

SYMTABC5 »1)=LPAR 

SYMTAB{5*2>=RPAR 

SYMTAB(6»1)=0WE 

SYMTABJ6»2)=MINUS 

SYMTAB(7,1)*EYE 

SYMTAB(7t2)=EQUAL 

SYMTAB«3,1)=EGGS 

SYMTAB{8,2  5sMINUS 

SYMTAB(9»1)=0WE 

SYMTA3(9,2)=SLASH 

SYMTA J(lu,l)*OWE 

SYMTAt3(10,2)*STAR 
C 

1    READ  (5.1000)  COM,TOUR*LINES»INCHES  .IHYPSO 
1000   FORMAT? 12.2  1 1 »F460,2X . 1 1 ) 

C    ZERO  COUNTERS  AND  ARRAYS  USED  IN  FREQUENCY  TABLE 

JJ=0 
MM=0 

SUM»0aO 
TLOW=0 
T0UT=u 
DO  5    2-1 3.25 
5      TABE<I)=0»0 
N=ROWS*COLS 

C      ■»■»#*#*•»«*•»**•«•■!?•*•*'«■***  **«#*'»*•»*»#•&-#*««*•»** 

C      N  *  NUMBER  OF  DATA  POINTS 

WRI TE ( 6  » 1030 )  T I TL ,N  f  COLS, ROWS 
1030   F0RMAT(1H1,12A6»//,22H  MAP  IS  DEVELOPED  FRCM»I5t30K  GRIDDED  VALUES 
1, CONSISTING  OF  .I3.12H  COLUMNS  BY  ,I3,6H  ROWS«//J 
CALL  STATSCROwStCOLS9lv  0»lOO»Z»ZBAR,XPtZSlG,MAxZ»MINZ) 


C    DETERMINE  CONTOUR  LINE  ELEVATIONS 

C    TOUR  VALUE  DETERMINES  METHOD  USED  IN  DEFINING  CONTOURS 

C    **#***##***•*#*#*#**#***#************•**** 

IF  {TOUR.GT.2.0R.TOURoLE»0)  GO  TO  60 

GO  TO  (30,40) ,TOUR 
30     WRITE  (6,1060) 
1060   FORMAT  <25H  CONTOUR  LEVELS  SPECIFIED) 

READ  < 5,1070)  ZMIN.ZMAX 
1070   FORMAT  (2F10..0) 

GO  TO  70 
40     READ  (5,1020?  I NCON 
1020   FORMAT  (10A6! 

WRITE  (6,1080) 
1080   FORMA*  (37H  VARIABLE  CONTOUR  INTERVALS  SPECIFIED} 

READ  »5»INC0N)  < DRANG f ! ) » 1 =1 oCON) 

ZMIN=DRANGU: 

ZMAX=DRANG(CONJ 

GO  TO  70 
60     WRITE  (6,1100) 
1100   FORMAT  (26H  CONTOUR  LEVELS  CALCULATED) 

ZMIN=M1NZ 

ZMAX=MAXZ 
70     WRITE(6»1110)  ROWS»COLS»ZBAR»ZSIG»MAXZ,MlNZ 

lllu   FORMAT  (13H  GRIDDED  DA?A/,42H  LINEAR  INTERPOLATION  WITHIN  GRID  S 
1ARES  /,17K  OBSERVED  VALUES  /,   8H  ROWS  =  <,  I3»1QH    COLS  ■  • 
213,/, 8H  ZBAR  =  »F10.3»30H    ZSIG  =  »F10.3»10H    MAXZ  =  ,Fl0«3,l0H 
3   MINZ  =  »F1D83) 

C    CALCULATE  MAP  >/IDTH  IN  CHARACTERS  (COL) 
C     *#*»#**#*■#*  **#*#**#  *#**x  #***#«*#•#  *##*•*## 

IF  ( INCHES. LE.O.O)  COL=100 

IF  (INCHES  .GT«>  12,7)  COL=127 

IF    (INCHES«GT«0.0«AND.lNCHES«LEo12,7)    COLssrNCHPS*l0»0 

PS= ( COL-COLS ) / ( COL  3-1 ) 
C  #*•*##****«####***###**#*##*****•*#***##** 

C    TEST  GRID  SIZE— IF  GRID  SMALLER  THAN  2  CHARACTERS  BETWEEN  COLUMNS, 
C    AN  ERROR  STATEMENT  IS  PRiMTED  AND  PLOT  ABORTED 
C     *#***###*#****##**#  **###*#*sHt#*-!Hc-#**«*** 

IF  (PSoGTol)  GO  TO  80 

WRITE  £6,1120) 
1120   FORMAT  (25H  TOO  MANY  COLUMNS  IN  GRID) 

GO  TO  999 
C 

C    CALCULATE  REQUIRED  CONSTANTS 
C 
60     C0LMAX*PS*{C0LS-1)+C0LS 

YlNC  =  1.0/(  (FLOAT  (fS+1)  )M).6> 

GRIDaPS+1 

XINC*1.0/FLOAT(GRID) 

RANGZ=ZMAX-ZMIN 


RHIGH-ZMAX+RANGZ 

RLOW*ZMlN-RANGZ 

IF  (CON.EQ«0«,OR.CONoGTelO)  C0N=5 

FCON*CON 

CURANGZ/FCON 

IF  (TOUR.EO«2l  GO  TO  ICO 

JJj=CON+l 

DO  90  1=1, JJJ 
90     DRANG(I)=ZMIN+FL0ATCI-1)*CI 
100    CONTINUE 

WRITE  (6,1130) 
1130   FORMAT  (1H  /,16H  VALUFS  EMPLOYED) 

WRITE  (6,1140)  ZMIN»2MAX»CI ,GRID»XINC tYlNC,CONtPS»COLMAX 
1140   FORMAT  { 8H  ZMIN  =  ,F10o3,10H    ZMAX  »  ,F10.3t8H    CI  *  »Fl0a3/,8H 
1GRID  =  »I5,10H    XINC  =  ,F10e3,10H    YlNC  ■  9F10e3s/»7H  CON  »  15, 
2  8H    PS  =  ,  I5»12H    COLMAX  *  tI5) 
C     ***#####**#^  #***■**#*****■#  *#*##***#****#* 
C    WRITE  OUT  CONTOUR  SYMBOL  TABLE 
C    **#*#*##**####*»#**#*#*#####*##*#*•**#*** 

WRITE  (6,1150) 
1150   FORMAT  (35H4C0NT0UR  SYMBOL  TABLE  IS  AS  FOLLOWS) 

WRITE  (6,1170) 
1170   FORMAT  ( 7H0SYMB0L, 11X»5HVALUE) 
110    WRITE  (6,1210) 
1210   F0RMAT(5H0    *»25X,7HN0  DATA,/,5H     -,23X»l0H8ELOW   LOW) 

DO  450  I*1«C0N 

WRITE  (6,1221)  SYMTABCI»2) 
1221   P0RMAT(1H+,4X»A1) 

WRITE (6, 1220)   SYMTAB( ! »1 ) , DRANG < I ) » DRANG ( 1+1 ) 
1220   F0RMATC1H  »4X » Al ,2X »F12e4,6H   TO   »F12»4) 
450    CONTINUE 

WRITE  (6,1230)  DRANGCCQN+1) 
123u   FORMAT  (4X«lH+,2XtF12.4*10H   OR   OVER) 
120    WRITE(6,1030)  TITL,N, COLS, ROWS 

XSTOP  =  GRID 

GRID=C0LS-1 

C0L=3 

DO  140  JsliGRID 

DO  13  )  I  «  2, XSTOP 

SYM( COL) "BLANK 
130    C0L=C0L+1 

SYMC COL) "PLUS 
140    COL-COL+1 

N=C0LMAX*2 

SYM(2)=PLUS 

SYM(1)=6LANK 

SYM(N)=BLANK 

WRITE  (6,1240)  CSYMC I > *I»1»N> 
1240   FORMAT  (1H9,129A1) 


c 

C    ESTIMATE  VALUE  OF  ALL  PRINTING  POSITIONS  ON  MAP 

C 

C      LOOP  THROUGH  MAP  LENGTH*  ROW  BY  ROW 

C 

ITOP=l 

I  ROW  =  ROWS  +  1 

DO  340  MM  =  2?IROW 

YSTOP  =  l./YINC  +  .•>& 

IF  (MM.EQ.IROW)  YST0P=1 

M*MM-1 
145    YP=0»0 

PSYM(1»1)=PLUS 

PSYMtl»2»»BLANK 

PSYMCN»ll=PLUS 

PSYM(N»2)=BLr. 


00  340  K  =  Ir/VSTOP 

150 

C0L  =  1 

C 
C 
C 
190 

21  =  YP*(2CMM91)-ZCM,1>J+Z(M»1) 

LOOP  THROUGH  MAP  WIDTH 9 COLUMN  3Y  COLUMN 

II  =  C0«.S-1 

JLFT=! 

DO  33J  J=1*I I 

JJ=J+1 

Z2= YP*( Z ( MM  »JJ) -Z (M  a J J) i +Z { M  » J J ) 

ZDELT=Z2-Z1 

xp=o.o 

DO  480  KLK  ■  UXSTOP 

JLFT=JLFT-l 

195 

COLsCOL+1 

200 

ZEST=XP*ZDELT+Z1 

IFCZEST.GE.ZMIN.AN9.ZEST.LE.ZMAX)  GO  TO 

IF(ZEST.GT.ZMAX)  GO  TO  : 

IF  (ZEST.LE.RLOW)  GO  TO 

I F { I  TOP • EQ * 1 .OR . MM *  EQ • 1  ROW * OR • JLFT • EG . 1 ) 

PSYM<C0L»1)*BLANK 

GO  TO  202 

201 

PSYM(C0L»1) "MINUS 

202 

PS YM I  COL u2)  =  BLANK 

TL0W=TL0W+1 

ZPREVCCOL)=ZEST 

ZLAST*ZEST 

GO  TO  124 

203 

IF(ZEST.GE.RHIGH)  GO  TO  206 

325 
GO  T0201 


IF( ITOP.EQ.l.OR.MV.EQ.IROW.OR.JLFT.EQ.l)  GO  TO204 

PSYM(COL»l)=BLANK 

GO  TO  205 

204  PSYM(COL,l)»PLUS 

205  PSYM(COL,2)=BLANK 
TABF(CON+l)-~TABF£C0N+l)  +  1.0 
ZPREV(COL)sZEST 
ZLAST-ZEST 

GO  TO  124 

206  IF(  ITOP.E(J.1.0R,,.';M.EQ.!kOW»OR.JLFT»EQ»l)  GO  TO207 
PSYM(COL»l)=BLANK 

GO  TO  20  8 

207  PSYM(COLtl)=STAR 

208  PSYM{C0L»2)=BLANK 
T0UT=T0UT+1 
ZPREV(COLJ=ZEST 
ZLAST=ZEST 

GO  TO  480 

325  CONTINUE 
LVAL=LINES+1 

IF{  ITOP.EQ.l.OR.MM.EQ.  1  ROW.OR.  JLFT«EQ«l  )L\/AL=2 
GO  TO  (326*327), LVAL 

326  DO  309  1*1, CON 
ZSTEP=DRANG(I! 
SW1=1. 
SW2=1. 
SW3=1» 
SW4=1. 

IF  (ZSTEP»G~sZLAST«.OR,ZSTEP«,GEoZPREV(COL)  )    SWlx0« 

IF  (ZSTEPoLE.ZFST)  SW2=0o 

IF  (ZSTEP»LT.ZLAST.OR,ZSTEPsLT«ZPREV{COL) >  SW3«0. 

IF  <2STEP.6E.ZEST)  SW4»0. 

IF(  <SWl.EQ.O..AND.SW2.EQ.0.).OR.(SW3.£Q.(H.AND.SW4.EG.0«>  )GO  T0310 

309  CONTINUE 

327  I«(I2EST-2MIN)/CI)+1. 
PSYM(COL»l)*SYMTAB{ 1,1) 
PSYM(C0L,2 )*SYMTAB( I»2) 
TABE(I)=TABE(I)+leC 

GO  TO  3  20 

310  PSYM|COL»l)*BLANK 
PSYM(COL»2)=BLAN!< 
TABE(I)  =  TABEU)+1.0 

320    ZLAST=ZEST 

ZPREV(COL)sZEST 
124    SUM=SUM+1#0 
48o  XP=XP+XlNC 

Z1  =  Z2 
330    JLFTa  ) 

COL  =  COL  +  1 

IF(2{M,COLS).GE.2MIN.AND.2(M»COLS).LE.2MAX)  GO  TO  333 

IF(2CM»COLS).GT.2MAX)  GO  TO  335 

IF(Z(M,COLS)«LE.RLOW)  GO  TO  336 

PSYM(C0L,1)=MINUS 

PSYM(COL,2)»3LANK 

TL0W=TL0W+1 

GO  TO  332 


335  IF<Z(M»C0LSUGE«RHIGH>    GO    TO    335 
PSYM(COL,l)>=PLUS 
PSYM(COL<.2)-3LANK 
TABE(CON+l}=TABE(CON+l)+loO 

GO    TO    332 

336  PSYM{COL*l)*STAR 
PSYM(COL»2)=BLANK 
TOUT=TOUT+l 

GO  TO  334 

333  I    =    ((Z<M»C0LS)-ZMIN)/CI>+1. 
PSYM ( COL  » 1 ) =SYMTAB ( I 9 1 ) 
PSYM(COU2)  =  SYMTAB(I,2) 

332         SUM=SUM+1.0 

334  WRITE(6»1241)     { PSYM ( L,2 ) »L=1 »N) 
1241       frOP,MAT(lh'+»12QAl) 

WRITE (6* 1240)  (PSYM< L,l ) »L=1 »N> 
YP=YP+YINC 

PSYMU»1)=8LANK 

PSYM(le2)=8LAN!< 

PSYM(N»1)=BLANK 

PSYM(N*2)=BLANK 
340    ITOP=0 

WRITE  (6»1240)  (SYM(L) »L=1»N> 
C     ###*#*#****#*#**###***###-*#***********## 

C    CALCULATE  AND  PRINT  FREQUENCY  TABLE 

C     ***##♦**##*###**##**#*#****###•*#*•***■»*#■* 

WRITE  (6*1250)  TOUT 
1250   FORMAT  ( 16H1FREQUENCY  TA3LE*/ » 12X •8HINTERVAL  ,16X»5HUNIVS  *6X»l0HP 
1ERCENTAGE  ,4X» lOHCUMULAT I VE  »/»16X»7HN0  DATA  ,UXtI8) 

PTAB=(FL0AT( TL0W)/5UM)*1 JC«0 

CTAB=PTA8 

WRITE    (6,1260)     TL0W9PTAB»CTAB 
126u      FORMAL    (14X»9H8EL0W    LOW* 11X *  18 »2 ( 3X »F12«4  )  ) 

DO    35J    I-UCON 

VALUF=F>RANG( I) 

VAL1=DRAMG( 1+1) 

VAL1=DRANGU  +  1) 

PTAB=(TABE( I ) /SUM )*100.0 

CTAB=CTAB+PTAB 
350    WRITE  (6*127,)  VALUE. VAL1 «TABF <I > »PTA8»CT£B 
127o   FORMAT  (!X9Fl2e4,6H   TO   »F12.4*3X >F8.0«2 ( 3X «F12.4) ) 

I=C0N+1 

PTAB=(TABF( I ) /SUM )*100 ,0 

CTAB=PTAB+CTAB 

WRITE  (6*1280)  VAL1 »TABEC I ) ,PTAB»CTAB ,SUM «CTAB,CTAB 
128u   FORMAT  ( 1X*F12«.4, 10H   CR   OVER  » 11X »F8.0 >2 ( 3X »F12.4 ) ,/ * 1 3X » 5HT0TAL 

l»llX»F8e0*2(3X»F12o4)) 
C     *###*#*♦«•#*******##*■**«•*#*«■****#**» 
C    IF  IHYPSO  ■  1  »  PERFORM  HYPSOMETRIC  ANALYSIS 
C     ##**»#*#■**###*##*#**•**•»#*  **■#**»#*♦# 

IF    (IHYPSO)    450»999945f 
450         CALL    HYPSO{DRANGtZV.Ax»Zi'iINtCON»CUMH.TITL) 
999         CONTINUE 

RETURN 

END 


C  SUBROUTINE  HYPSO 

c  „ „„ — . ~_„ — „ ... _„„..-„„.-_— 

C  PURPOSE  - 

C  TO  PERFORM  HYPSOMETRIC  (  AREA-ALT ITUDF  )  ANALYSIS. 

C 

C        REFERENCE- 

C  A.  N«  STRAHLER-1952-HYPSOMETRIC  i AREA-ALTITUDE  J  ANALYSIS 

C  OF  EROSIONAL  TCPOGRAPHY-GEOL.  SOC»  AMERICA  BULL*  VOL«63 

C  PAGES  1117-1142  * 

C  USAGE  - 

C  CALL  HYPSO(DRAN6»Z.MAX»2MIN»C0N»CU.,<H»TITL} 

C         DRANG-ARRAY  CONTAINING  CONTOUR  LEVELS* 

C         ZMAX  -MAXIMUM  ELEVATION  VALUE* 

C         ZMIN  -MINIMUM  FLFVftTION  VALUE* 

C         CON   -NUMBER  OF  CONTOUR  INTERVALS© 

C         CUMH  -ARRAY  CONTAINING  CUMULATIVE  PCT  AREAS  BELOW  EACH  CONTOUR* 

C         TITL  -  ALPHANUMERIC  TITLE  (  1-72  CHARACTERS) • 

C 


SU3ROUTtNE  HYPSOC  DRANGsZMAX  »ZMIN»CON  »CUMH  *TI  TL.5 

DIMENSION  IR0WC49.81) .SELHC25) »RELA{25) 

0 1  MENS I ON  T I TL ( 1 2 ) , CUMH (27). DRANG  I  27  J 

DATA  D0T»BLANK»STAR/1H.,»1H  »1H*/ 

REAL  I ROW 

INTEGER  CON 

C       ##•»»*#*##**##*#*#**#*■»**«•##*■«•«■•»■##**#«**•«••!*•*-}** 

INTEGER  CON 
C     *«•##*####  ###•*  #*####*#*  *********#*•*####***  •*■-*■** 

C    MAKE  PLOT  ALL  BLANKS »DRAy  BORDERS  WITH  DOTS»AND  PLOT  A  STAR  AT  TOP 

C    LFFT  AN >  BOTTOM  RIGHT  CORNERS 

C     #«•#***  >H:#  #*##*###******  «■*•»«■■»■#*«■  ♦#****•##*##*-**» 

DO  2  M=l*49 
DO  2  N=l»82 

2  IROW(M.N]«BLANK 
IROW(l»l)=STAR 
DO  3  N^2»81 

3  IROW(UN)  =  DOT 
IROW(49»81)=5TAR 
DO  4  N-l»80 

A      IROW(49»N)=OOT 
DO  5  M=2»Ae 
JR0W(M91)=D0T 

5      IROW(M»81)sDOT 

C     ***#*##*****#*####*#*♦****#**###*###•»•*##*#■*•** 

C    DETERMINF  POSITIONS  OF  STARS  WITHIN  THF  PLOT 

C       *#*#♦**#****#**#***#*##**#*##********#»***•*** 

K=CON+l 

DO  20  1=1 #K 


RELH< I >=f DRANG( I ) -ZMIN J / I ZMAX-ZMl N ) 

REL  A  d  )  » 1  •  <■-  U  CUMH  (  I  )  -CUKH  ( 1  )  )  /  (CUMH  ( C0N+ 1 )  -ClJMH  {  1  3  )  ) 

M=(49.0~(RELH{ I )*49.0) }+0,5 

N=<RELA< I )*81.0)+0.5 
20     IROW(M»NJ«STAR 

C     ##*****#**«■«•#){•****#*****■;■*»#«***#  it****^**-**** 
C    PRINT  OUT  ALL  THE  PLOT  SEQUENTIALLY 
C     *#*####-s.-K-#*##**##*';i#*#*H  «•#-*■*•*■«•■»#**■»•#*■#** * »**» 

WRITE(6durO)  TITL 

1000  FORMAT  (lHltZZHHYPSOMeTB.'C  CURVE  FOR  ,12A6l 
WRITE(6>1001) 

1001  F0RVAT(1H  ,8HRELe  HT«/) 
WRlTEC6»lw021 

1002  F0RMAT(1H  ,  8X*3H()  »0  »37X  »'-  H0s5  »37X  »3H1  ,0  ) 
WRITE  (6 » 1003) 

1003  F0RMAK1H  »9X»lH+t  39X»lH-i-»39X»lH+) 
WRITE  (6 *  1004  3  {  IR0W£  1»N)  *v>1<>81) 

1004  FORMA  T(1H  S4X»5H1»0  +  ,8)A!.,5H+  leO) 
00  10j  M=2,24 

100  WRITF(6»10f)5)  {  IROW.M»Ki  si  «1  »8l) 

1005  FORMAT  (1H  ,9X*81A].) 
WRITE(6»1006)  <IROW(25«N)  N=l,81) 

1006  F0RMAK1H  ,4X,5H0„5  +,81£.,5H+  0*5) 
DO  101  M=26*48 

101  WRITE(6»luC5)  ( IR0W«M,N)*f-l*8l) 
WR  I  TF  (  6  »  1  •-'  C  7  )  (I  ROW  (  4?  •  N  >  i  N*>  1  •  8 1  > 

1007  FORMATflH    »4X»5H0»0    +,81A1,5H+    0*0) 
WRITF(6»1003) 

WRITE16»1002> 
WRITE(6»1U08) 

1008  FORMAT  (1H  ?44X » 13HRELATP  .'  AREA  /) 

C      ***#**«-#«**-K->##**iJ#-***W**»r  #*•»■•«■  ##*•■«•  *#*•■&  ■»;■•*#*# 

C    CALCULATE  AREA  BELOW  THE  CUFv  E  USING  DECIMAL  VALUES  OF  PQk\TS  FOR 

C    THF  P^ST  OBTAINAPLF  ACCURACY,, 

C     **»♦#*#**#*##***##«##*#*#*«****###****#*«•##*# 

AREA=0.0 

DO  200  I «if CON 

A={  (RELHI  I  J+RELHJ  1-8-1)  )*(RE.A(  I)-RELAd  +  X)  J  )/2o0 
200    AREA=AREA+A 

WRITEC6tl010)  AREA 
lOlu   FORMAT (1H  ,77HTHE  HYPS0MET:IC  INTEGRAL  FOR  THE  AbOvE  CURVE  ,  THE  A 
IRE  A  BELOW  THE  CURVE  )  iS  »-".»4) 

C      *#**###***•»>##*#***■**■**  »•****  1  ft********** 

C    WRITE  OUT  COORDINATES  OF  HYP :-C  '.ETRlC  CURVE 

WRITE  {6»2001) 
2001   F0RMAT(1H1,50HCALCULATED  V<LL!S  OF  HYPSOMETRIC  CuRVE  C00RDIf\TES» 
1//»5X?12HP0INT  NuMBER»5Xtbt-A-\xlStlOX*6HH-AXlS  ) 

DO  21  I-lsCON 
21     WRITE  {6,2^00)  I«RELAd  >fR?.H(*) 
2000   FORMAT  ( lHu  ,  10X  »  I  2  , 1  Ox  * F6 e  4,10X  »»"S •  4  ) 

RETURN 

END 


C  PROGRAM  MAPPER 

C    PURPOSE  - 

C  A  SIMPLE  PROGRAM  TO  READ  IN  GRIDDED  DATA  FOR  CONTOURING* 

C 

C  CONTROL  CARD  ORDER  - 

C 

C  1)  SPECIFICATION  CARD  - 

C         ROWS4CGLS1-3)   NUMBER  OF  ROWS  IN  GRID*  (MAX"lOO) 

C         COL5«CQLS4-63   NU'IRER  OF  COLS  IN  GRID*  <MAX»100) 

C         LISTCCOL  75    IF  LI5T=1»  LIST  INPUT  DATA 

C         NMAPCCOLS3-9)  NUMBER  OF  MAPS  REQUSTED  FOR  THIS  DATA  SET, 

C 

C  2)  TITLE  CARD  - 

C         ANY  ALPHANUMERIC  TITLE  IN  COLUMNS  1-72- 

C 

C  3.  FORMAT  CARD  - 

C         SPECIF  I FS  FORMAT  OF  INPUT  DATA  (COLS  1-60) 

C 

C  4)  DATA  CARDS 

C 

PROGRAM  MAPPER? I NPUT, OUTPUT »TAPE5=I NPUT ,TAPE6^0UTPUT ) 

DIMENSION  FORM(10),TTTU 12 ) »Z( 100, 100 ) 

INTEGER  ROWS* COLS 
1      READ  (5>lU00)  ROWS»COLS»LIST»NMAP 
1000   FORMAT  (2I3*11»I2) 

IF  (c'0F»5)  999  5,2 
?      IF  (NMAP.LT.l)  NMAP^l 

READ  C5»10105  TlTL 
lOlu   FORMAT  (12A6) 

READ  C5»lO20)  FORM 
102'J      FORMAT  (10A6) 

DO  7  1  si, ROWS 
7      READ  C5»F0RMJ  {Zil » J) » J*l »COLS> 
C     *#*******«•*###****»■*##■»##**#**#»■»•###*#♦♦ 

C    IF  LIST  =1  LIST  INPUT  DATAsOTHERWlSE  GO  TO  20 
C    *#*#**##*******##»#■*####*##**#**»#**#### 

IF  (LIST»NE.1»  GO  TO  2C 

WRITF  (6tl030)  TITL 
1030   FORMA"  (1H1«12A6//) 

WRITE«6tlu25) 
1025   FORMAT(lH0»25Xf22KLlSTING  Or  INPUT  DATA,   ) 

DO  10  1=1 ♦  ROWS 
10     WRITE  (6»1C40)  (Z ( I *J) i J=l «COLS) 
1040   FORMAT(1HO,10F10.3»/(1X»10F10#3) J 
20     DO  25  I=1,NMAP 
25     CALL  PRTC0N1(Z»R0WS»C0LS»TITL) 

GO  TO  1 
999    STOP 

END 


C  SAMPLE  DATA  FOR  PROGRAM  MAPPER , SUBROUTINE  PRTCON1 

20  161  3 
DRAINAGE  5ASIN»  VERDUGO  HlLLSs  CALIFORNIA,, 
(16F4.0) 

99999 999999999 9999999999999999999999c9999999999 99999 188 099999999 
99999999999999 999 9999 99999999999 9999999999999999 182 5 186099999999 
999  99999999999999  99999<?99  99Q9  99999  9999999999 177  5 177  5 183099999999 
9999999999999999999999999999999999999999171517001775180099999999 
99999999999999  99999999999999161016001613162516751680174517609999 
999999999999999999999999157^157515  50152515751628163  8171217509999 
9999999999999999 153 31 550152 31 500 1488 14?5 1525 1525 16201 6S01700 2 712 
99999999999999991490 1483 14601455 1425 144014501 5 50  U>25 1600 16301675 
999999999999147514251450142  5140013751380150015751560155516189999 
99999999999914401400138  51375137013251400151214801475153299999999 
9999Q99999991 390 1375 132 5 1320 13 12 132 5 141 2 1425 1412 148 2 151 299999999 
9999999999991 3  80 13401 30012751 255 1325 1325136014401475999999999999 
99999999140'Jl  345 1280 125  5 1220 1225 12651325 139014509999999999999999 
999999991 3601 300 1250 1200 116 51 2 50 1325 13801425 99999999999999999999 
9999999913321275121211751208129013801395999999999999999999999999 
999 999991 312 1230 1155 11 66 1250 1325 9999  999  999999999999 9999999999999 
99999999 126i>  1170 112  5 1200 1275 132599999999999999999999999999999999 
99999999! 225 11 25 11 25 12 20 126299999999999999999999999^999999999999 
99999999 1175 11 20 UOOll 75999999999999999999999999999 9999999999999 
9999999999999999108099999999999999999999999999999999999999999999 
16IU10.0 


1080o 

*.880. 

161110.0 

1080. 

1880. 

161210.0 

1080. 

188ue 

SAMPLE  OUTPUT  FROM  SUBROUTINE  PRTCON 


65A 


CONTOUR  MAP  OF  DRAINAGE  BASIN  IN  VF.ROUGO  HILLS. STRAHLER  FIG  15  . 

320  OBSERVATIONS 

CONTOUR  LEVELS  SPECIFIED 

GRIDDEO  DATA 

LINEAR  INTERPOLATION  WITHIN  GRID  SQUARES 

OBSERVED  VALUES 

(16F4.0) 

ROWS  «   20    COLS  ■   lh 

ZBAR  a    6145.009    ZSIQ  »    4262.532    MAxZ  ■    9999,000    MlNZ 


1080.000 


VALUES  EMPLOYEO 

ZMIN  ■    1080.000    ZMAX  =    1880.000    CI  = 
GRID  ■      *    XINC  b       0,167    YINC  * 
CON  ■     16    PS  ■      5    COLMAX  ■     91 
CONTOUR  SYMBOL  TABLE  IS  AS  FOLLOWS 

(ONLY  EVEN  SYMBOLS  PRINTFD  IN  BODY  OF  MAP) 


50.000 
0,278 


SYMBOL 


VALUE 


• 

NO 

DATA 

m 

HFLOW   LOW 

0      10B0 

.0000 

TO      1130, 

,0000 

1      1130 

.0000 

TO      1180. 

,0000 

2      1130 

,0000 

TO      1230, 

,0000 

3      1230 

.0000 

TO      1280, 

,0000 

4      1280 

.0000 

TO      1330, 

,0000 

5      1330 

.0000 

TO      1380, 

,0000 

6      1380 

.0000 

TO      1430, 

,0000 

7      1430 

.0000 

TO      1480, 

.0000 

8      1480 

.0000 

TO      1530, 

,0000 

9      1530 

.0000 

TO      1580, 

,0000 

_  A              1580 

.0000 

TO      1630, 

,0000 

8      1610 

.0000 

TO      1680, 

,0000 

C      1680 

.0000 

TO      1730, 

,0000 

D      1730 

.0000 

TO      1780, 

,0000 

E      1780 

.0000 

TO      1B30, 

,0000 

F      1810 

.0000 

TO      1880, 

.0000 

4      1880 

.0000 

OR   OVER 

FREQUENCY  TABLE 

INTERVAL 

UNITS 

PERCENTAGE 

CUMULATIVE 

NO  DATA 

3081 

BELOW 

LOW 

0 

0. 

0. 

11)80. 0000 

TO 

1130.0000 

39. 

1.7751 

1.7751 

1130.0000 

TO 

1180,0000 

98. 

4.4606 

6.2358 

1180. 0000 

TO 

1230,0000 

171. 

7.7833 

14.0191 

1230.0000 

TO 

1280.0000 

184. 

8.3751 

22.3942 

1280.0000 

TO 

1330.0000 

223. 

10.1502 

32,5444 

1330.0000 

TO 

1380.0000 

230. 

10.4688 

43.0132 

1360.0000 

TO 

1430.0000 

230. 

10.4688 

53.4820 

1430.0000 

TO 

1480.0000 

213. 

9.6950 

63.1771 

1480,0000 

TO 

1530.0000 

182. 

8.2840 

71.4611 

1S3O-.O000 

TO 

1580.0000 

146. 

6.6454 

78.1065 

1580.0000 

TO 

1630.0000 

136. 

6.1903 

84,2968 

1630.0000 

TO 

1680.0000 

88. 

4.0055 

88,3022 

1840,0000 

TO 

1730.0000 

84. 

3.8234 

92.1256 

1730.0000 

To 

1780.0000 

67. 

3.0496 

95.1752 

1780.0000 

TO 

1830.0000 

34. 

1.5476 

96.7228 

1830,0000 

TO 

1880.0000 

17* 

0.7738 

97.4966 

1880.0000 

OH   OVFR 

55. 

2.503* 

100.0000 

■     '   :?. 

TOTAL 

2197. 

100.0000 

100,0000 

65B 
SAMPLE    OUTPUT    FROM    SUBROUTINES    PRTCON   I    AND  HYPSO 

CONTOU*    «A»    Of    0HAINA8E    8ASTN    IN    »f»OU«0    Hi  LIS . STMAHLEO    FI<1    If    . 


cccccce 
cce 


ui  .  _•_ 

ECCCC  • 

KEH  •• 

EEEEEE      t__ 

CCCCCC  • 

CCCCCCCCC     EEEEE  •• 

CCCCCCCCCC        t        f_ 

CCCCCCCC  • 

.......1.11.. ......        ccccc  •• 

AAAAAAAAAAAAAAAAAAAAA              CCCC  * 

AAAA                AAAAAAA                CCCC  • 

88         AAAAAAA            CCCCC  •♦ 

89688888        AAAAAAAAA      CCCCCC  _f 

88889399888888           AAAAAA     CCCCCCCC  • 

88688888888889888    8888888888(8    AAAA     CCCCCCCCCCCCC* 

888888  88888888888888  888888888     AAAA 

8888888888888888  888888    AAAAAAA 

8888888                     «66              888      AAAAAAAAAAAAAAA  •  • 

8686666668        888     AAAAAAAAAAAAAAAA  • 

6666666666666666      888       AAAAAA    AAAAAA  • 

66           666666666666        666    688                    AAAA  •• 

6666666  66*66666666666         6666   888(8  •_ 

666666666666666666666            6666   68888888688888888  ♦ 

66666666666666                4     666    66886688       88888  •• 

6666666666                      64     666     86             8688(6  * 

66666666                       444     6666                  898866  • 

66666        4444444444444444444    66666666666666*6     8868916  •• 

66  44444444444444444444      86666666666666     ( t_ 

6           444444444444444444444           66*******  • 

6        444444444            4444444444        6666  •• 

444*4444  444444444       6666  l_ 

44444444                   44444444       •***  * 

666      44444            72222222        44444      666*6  ». 

6    44444      2222222222    44444    **£  •_ 

444444       22222722222      4444*     »*«*  • 
44444      22222222222222    4444 
444*4       22222222222222     6444       6E 
44444     222222222222222    4444      46 
44444      2222222  222222?     444    6666666 
44444    22722222  222222    44C 


4444     22272 

222?? 

44* 

444     7772 

222 

4444 

4*   2222 

2222 

4444 

4   2227 

2222 

44444 

722      0 

2222 

444444 

232        0 

2222 

C 

222        00 

22222 

22?    0000000 

2222 

222    0000000 

224 

?2    oooooooo 

2? 

oooooonooo 

908 

*Y*SOMCT*ic  cuhve  roil     contouh  >»  or  mainasc  basin  in  ve»ou(0  hillIiStianlcr  fk  is 
ML.    HT. 

0.0  0.5 


0.0  0.5 

«LATI»£    AREA  _ 

M»P50NET»!C     INTEOMAL    FON    TNC    ABOVE    CUUVE     {     TNC    AREA    BEUO.    TMf    CUOVl    I     II    0.4J85 


PROECCTION  OF  CGETOSIR  MPS 
OS  SEE  CAICQMP  PLOTTER 

PROGRAM!  COOTJR 


C  PROGRAM  CCNTUR  -  GENERAL  INSTRUCTIONS 

C 

C 

C  PURPOSE  - 

C         TO  PLOT  CONTOUR  MAPS  OF  GRIDDED  DATA  ON  THE  CALCOMP  PLOTTER. 

C  THIS  PROGRAM  IS  RASED  ON  A  FORTRAN  II  PROGRAM  WRITTEN  BY  FRANK  Jo 

C  RENSt UNIVERSITY  OF  MICHIGAN  DEPARTMENT  OF  GEOGRAPHY. 

C 

C  PURDUE  UNIVERSITY  VERSION  1.2 

C  PROGRAMMED  BY  —     SANDRA  L.  TURNER®    A„  KEITH  TURNER, 

C  COMPUTER  SCIENCE.    CIVIL  ENGINEERING* 

C  FEBRUARY  1968, 

C 

C         THIS  PROGRAM  CAN  PRODUCE  MAPS  TO  ANY  SCALE  WITHIN  THE  LIMITS 

C  OF  THE  PLOTTER  PAPER  WI^TH  (MAXIMUM  ~    28  INCHES). 

C 

C         MANY  OPTIONS  ARE  AVAILABLE  ANU  CAN  bE  SPECIFIED  THROUGH  THE  USE 

C  OF  CONTROL  CARDS*  THE  USER  MAY  SPECIFY  CONTOURS  AT  REGULAR  OR  IRREG- 

C  ULAR  INTERVALS* AND  MAY  TRANSFORM  MAP  COORDINATES  BY  TRANSLATION  OR 

C  ROTATION  ABOUT  ANY  OF  THE  THREE  MAJOR  AXES  (  X?>Ys>GR  Z3„ 

C 

C         THE  PROGRAM  CAN  PRODUCE  ISOMETRIC  BLOCK  DIAGRAMS  OF  THE  CONTOUR 

C  MAPS  BY  ROTATING  THE  DATA  ABOUT  ALL  THREE  AXFS  BEFORE  PLOTTING.  THE 

C  USE  OF  A  PERSPECTIVE  OPTION  IN  COMBINATION  WITH  TRANSLATION  WILL 

C  ALLOW  THE  PLOTTING  OF  LEFT  AND  RIGHT  VIEWS  FOR  STEREOSCOPIC  {  3-D) 

C  EFFECTS. 

C 

C         THE  PROGRAM  IS  CAPABLE  OF  PLOTTING  A  NUMBER  OF  DIFFERENT  MAPS 

C  OF  A  SINGLE  SET  OF  DATA,  EACH  INVOLVING  DIFFERENT  MAP  SCALES. CONTOUR 

C  INTERvALSdROTATIONS  OR  TRANSLATIONS*  WITHOUT  REREADING  THE  DATA, 

C  VALUES  OF  CONTOUR  LINE  COORDINATES  ARE  CALCULATED  AND  STORED  ON  A 

C  SCRATCH  TAPE.  THESE  VALUES  ARE  MODIFIED  AND  PLOTTED  FOR  EACH  MAP 

C  REQUESTED  SO  THAT  THE  CALCULATIONS  ARE  MINIMIZED  AND  THE  PLOTTING 

C  OF  MULTIPLE  MAPS  IS  EFFICIENT* 

C  TAPE  REQUIREMENTS  - 

C  DA  SCRATCH  TAPE  (  TAPE  3)  IS  USED  TO  STORE  CONTOUR  COORDINATES 

C  PRIOR  TO  PLOTTING* 

C  2!  A  PLOT  TAPE  IS  REQUIRED. 

C         SINCP  TWO  TAPES  ARF  USFn.RFASST  '^NMENT  OF  TAPE  UNITS  IS  REQUIRED 

C         THIS  IS  ACCOMPLISHED  BY  A  MAP  ROUT INE • INCLUDED  IN  THE  LISTING. 

C 

C  ROUTINES  REQUIRED  - 

C  A)  MAIN  PROGRAM  (COMTUR)  PERFORMS  OPERATIONS  ACCORDING  TO  CONTROL 

C  CARD  SPECIFICATIONS. 

C  B)  SUBTOUTINES  - 

C         1J  SCAN   -SCANS  DATA  TO  LOCATE  CONTOURS. 

C         2)  TRACE  -TRACES  OUT  CONTOUR  LINE. 

C         3)  CALC   -PERFORMS  INTERPOLATION. 

C         4.)  DRAFT  -PERFORMS  ROTATIONS  AND  TRANSLATIONS  »WR  I TE5  PLOT  TAPE. 

C         5)  PKOJ   -DUMMY  SUBROUTINE  TO  AU-OW  ADDITION  OF  USER-DEVELOPED 

C  MODIFICATIONS. 

C 


c 
c 
c 
c 

c 
c 
c 
c 
c 
c 
c 
c 

c 
c 
c 
c 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


DESCRIPTION  OF  CONTROL  CARDS  AND  THFIR  OPERATIONS  - 


EACH 
SOME 


CONTROL  CARD  CONTAINS  A  CONTROL  INTEGER  « NCON )  IN  COLUMNS  1-2. 

OPERATIONS  REQUIRE  INPUT  DATA  IN  THE  FORM  OF  VARIABLES  AND/OR 
DATA  ON  ADDITIONAL  CARDS.  ANY  VARIABLES  REQUIRED  SHOULD  bE  PLACED 
IN  COLUMNS  3-12, 13-22 »23-32,33-42,  ETC  wlTH  DECIMAL  POINTS  PUNCHED. 
WHEN  THE  CONTROL  INTEGER  INDICATES  THAT  ADDITIONAL  CARDS  ARC 
ARY  FOR  A  COMPLETE  OPERATION,  THESE  CARDS  SHOULD  FOLLOW 


NECESS- 
IMMEDIATELY. 


NCON 
01 


REMARKS 


02 


03 


OA 


05 


08 


10 


OPERATION-  READ  AND  CONTOUR  DATA  MATRIX 
VARIABLES-1) NUMBER  OF  ROWS 

2) NUMBER  OF  COLUMNS 

3) WIDTH  OF  MAP  IN  INCHES 

MSC4LE  FACTOR  FOR  Z  VALUES  (Z=Z/ZSCALEl  (IF  BLANK=1> 
ADDL  CARDS-  THE  DATA  MATRIX  IS  READ  IN  ROW  BY  ROW  FROM  TOP  TO 

BOTTOM  {  EACH  ROW  IS  SCANNED  FROM  LEFT  TO  RIGHT). 

THE  MATRIX  IS  READ  ACCORDING  TO  A  PREVIOUSLY  DEFINED 

FORMAT  (SEE  NCCN*19)»  MAXIMUM  SIZE  IS  100X100  • 

OPERATION-ROTATE  THE  VIEWING  POSITION  OF  THE  SURFACE  THETA 

DEGREES  IN  THE  X  DIRECTION. 
VARIABLES-1)  THETA  IN  DEGREES- 
ADDL  CARDS-   NONE 

OPERATION-ROTATE  THE  V  EWlNG  POSITION  OF  THE  SURFACE  THETA 

DEGREES  IN  THE  Y  DIRECTION. 
VARIABLES-1)  THETA  IN  DEGREES* 
ADDL  CARDS-   NONE 

OPERATION-ROTATE  THE  VIEWING  POSITION  OF  TH^  SURFACE  THETA 

DEGREES  IN  THE  Z  DIRECTION* 
VARIABLES-1)  THETA  IN  DEGREES. 
ADDL  CARDS-   NONE 

OPcRATlON-CALL  SUBROUTINE  PROJ  FOR  USER-DEFINED  MODIFICATIONS. 
VARIABLES-  NONE 
ADDL  CARDS-   NONE 

OPERATION-PLOT  USERS  NAME  AND  ACCOUNT  NUMBER  ON  PLOT. 

VARIABLES-  NONE 

ADDL  CARDS-  NAME  AND  ACCOUNT  NUMBER 

OPERATION-PLOT  CONTOURS  WITH  PERSPECTIVE  TRANSFORMATION  TO  PROD 

-UCE  LEFT  OR  RIGHT  STEREO  MODELS. 
VARIABLES-1 JDISTANCE  TO  PERSPECTIVE  PLANE. 

2JDISTANCE  TO  OBJECT  PLANE. 
ADDL  CARDS-   NONE 

OPERATION-TRANSLATION  OF  CONTOUR  COORDINATES. 

VARIABLES-liTRAMSLATICN  OF  X  IN  INCHESINOT  SCALED, TRUE  MEASURE) 
2) TRANSLATION  OF  Y  IN  INCHES (NOT  SCALED, TRUE  MEASURE) 
31TRANSLATI0N  OF  Z  IN  INCHES(NOT  SCALED, TRUE  MEASURE) 

ADDL  CARDS-   NONE 


NCON  REMARKS 


11  OPERATION-  DEFINE  PLOTTER  TAPE  AS  OTHER  THAN  NORMAL  NUMBER. 
VARIAELES-1J  PLOTTER  TAPE  NUMBER 

ADDL  CARDS-   NOME 

12  OPERATION-  READ  TITLE  FOR  PLOT(PLOTTED  AT  TOP  OF  MAP). 
VARIABLES-  NONE 

ADDL  CARDS-  TITLE  PUNCHED  IN  COLUMNS  1-48 

13  OPERATION-  TERMINATE  PLOT  TAPE 
VARIABLES-  NONE 

ADDL  CARDS-   NOME 

14  OPERATION-  READ  LABELS  AND  COORDINATES {MAX*  PER  PLOT  =  10). 
VARIABLES-1)  X  COORDINATE  OF  L- 

21  Y  COORDINATE  OF    LABEL 

3)  Z  COORDINATE  OF  LABEL 

4)  HEIGHT  OF  LABEL  LETTERS. 

ADOL  CARDS-  A  LABEL  PUNCHED  IN  COLUMNS  1-3  6 

15  OPERATION-DEFINE  REGULARLY  SPACED  CONTOURS. 
VARIABLES-1 JSTAPT1N6  CONTOUR 

2)  CONTOUR  INTERVAL 

3)  FINAL  (HIGHEST)  CONTOUR 

ADDL  CARDS-   NONE 

16  OPERATION-REDEFINE  CONTOUR  INTERVALS  FOR  SAME  DATA, 
VARIABLES-  SAME  AS  FOR  NC0N»15 

17  OPERATION-DEFINE  IRREGULARLY  SPACED  CONTOURS* 
VARIABLES-1 5  NCONT-NU.MBER  OF  CONTOURS. 

ADDL  CARDS-1 3  FORMAT  OF  CONTOURS 

2) CONTOUR  LEVELS  AS  DEFINED  ABOVE. 

18  OPERATION-  SPECIFY  WIDTH  OF  PAPER <  11  OR  28  INCHES) 
VARIABLES-1)  WIDTH  OF  PAPER 

ADDL  CARDS-   NONE 

19  OPERATION-  SPECIFY  FORMAT  OF  DATA 
VARIABLES-  NONE 

ADDL  CARDS-FORMAT  OF  DATA  MATRIX  IN  COLUMNS  1-72 

20  OPERATION-PLOT  MAP  ACCORDING  TO  ALL  CONTROLS  GIVEN  SINCE  LAST 

CALL  TO  PLOT  (LAST  NCON»20) 
VARIABLES-  NONE 
ADDL  CARDS-  NONE 

21  OPERATION-  CHANGE  wiDTH  OF  PLOTTED  MAP  (CHANGE  SCALE) 
VARIABLES-1)  NEW  WIDTH  IN  INCHES- 

ADDL  CARDS-  NONE 

COMMENTS  AND  RESTRICTIONS 

PRIOR  TO  READING  IN  DATA  MATRIX*  THE  FORMATCNCON-19 ) AND  CONTOUR 
INTERVALS(NCON=15  OR  17?  MUST  BF  DEFINED. 

BEFORE  CALL  TO  PLOT  (NCON=20)  PAPER  WIDTH  (NC0N»1B)  AND  ACCOUNT 
CODE  (NCON*06)  MUST  3E  DEFINED.  PLOTTING  THE  MAP  RESETS  ALL  CONTROLS 

TO  PLOT  A  NEW  SET  OF  DATA  USE  NCON*l  WITH  OR  WITHOUT  A  NEW 
FORMAT  OR  CONTOURING  CONTROLS. 


SIBSYS 

^ATTACH         CDOO/0 

SAS  SYSUT1»05C*05C»OC 

$*     CSC  PLOT  TAPE  REQD  ON  Al 

$$$    CSC  PLOT  TAPE  REOD  ON  Al 

SEXECUTE        IBJ08 

SlBJOb 

SIBFTC    COtolVft    DFCK 

DIMENSION    AM  (  100  »  100)  ,REC{800)  VX(1500)  <>Y  (  1500  )  t  IPT  (  3*3  )  *  INXC  8  )  > 

1  INY(8)»IOP{20)  tXL(lO)  tYL(lO)  fZL(lO)  *LAB(l0s>6)  *DX ( 2 ) 9DY( 2 ) ,DZ  (2  )  , 

2  HT(l0)9CTHETA<4) .STHETA(A) »TEMP(6) »FMT( 12) »Z< 1500) »TIT< 8) » 

3  CFMT(l2)»VCON(20)fWORK(1024) 
DIMENSION  USER(8) 

COMMON  MT»NTtNI»lX»lYt IDX>IDY»ISS»IT#lV»NP»N,JT»PYtRECtCVtX»YtIPTi 

1  lNXtlNY»DLtAM»IOP,DPtDOt ILABtXL»YL,ZL»LAB»DX»DY»DZtXMAX»HT»SCALE9 

2  YMAX»CTHETA»STHETA»FM2»FNN»NN»MtCL»D»NLlNES»NCuRV»Z*TlTtW0RK 
COMMON  MAXX»MAXYtMINX»MINY*SlZE 

COMMON  ZSCALE 

REAL  MAXX.MINX.MAXY9MINY 

NPL0T=1 

I0P(1)=U 

NLINFS*0 

ISTAP=5 

D«0« 
9999   ILAB=0 

REWIND  3 

DO  9981  1=2.20 
9981   IOP(I)«0 

NBOUND=0 
100U   READ  (5»luO)  NC0N»(TEMP(I)tIsl,6) 
100    FORMAT  <I2t6F10.0) 

WRITE  (6»1001)  NCON*(TEMP{ !) »I«1»6) 
1001  FORMAT  (3X9I596F10.3) 

IF(NCON.LE„  0  .OR.  NCON  .6T.  21)  GO  TO  30 

GO  TO  (I»2»2»2»5»6i30»8t8»10»ll»12»13»14»15»16»17»18tl9t20»21)» 
1  NCON 
C 

C  NC0N=1 
1      M=TFMP{1) 

NN=TEMP(2) 

TEMP(1)=TEMP(2)-1« 

SCALE=TEMP ( 1 ) /TEMP ( 3 ) 

IF  (TEMP(4).EQ«0»)  TEMP<4  5«1. 

ZSCALE=TEMP«4) 

FM2=FLOAT(M)/2. 

FNN=FL0AT(NN)/2. 

YMAX=(FM2/SCALE)+2.5 

XMAX=lFNN/SCALE>+2.5 

MAXX=o„J 

MAXY»0.0 

MlNX*l.E15 

MINY=1.E15 


IF  (IOP(l)  .EO.  0)  GO  TO  5027 

I0P(15)=1 

IOP(19)=l 

5027 

I0P<1)=1 

IF(IOP(19)  .EQ.  0)  GO  TO  902 

DO  2026  I=1»M 

2026 

READ  CISTAP.FMT)  ( AM( JVI ) «J»1 VNN> 

WRITE(6»FMT)  ( (AM(J»I),J=1,NN),I=1,M) 

IF(IOP<15)  .EQ.  0)  GO  TO  903 

C 

C  VAR 

IF  (IOPV  .EQ.  0)  GO  TO  2065 

IABLF  CONTOURS 

2027 

DO  2069  I=l»NVCON 

CL=VCON( I) 

CV=CL 

CALL  SCAN 

2069 

CONTINUE 

C 

C  REG 

GO  TO  200 

ULAR  CONTOURS 

2065 

CV*CBGN 

NLINESeO 

2U66 

CL=CV 

CALL  SCAN 

CVeCV+ClNC 

IF  (CV  .LE.  CEND)  GO  TO  2066 

200 

END  FILE  3 

2068 

REWIND  3 

C 

C  NCOI 

GO  TO  1000 

M=2»3»4 

C  DEFINE  THETA  FOR  ROTATION 

2 

I0P(NC0N)=1 

NB0UND=1 

THETA=-TEMP(l)/57. 295795 
STHETA(NCON)=SIN( THETA) 
CTHFTA(NCON)=COS( THETA) 
GO  TO  loOO 

C 

C  NCON=5 

5      IOP(NCON)=l 
GO  TO  1-000 


6      I0PCNC0N1»1 

READ  (5*6000)  <USCr!(  I)  *!*1»8) 
600u   FORMAT  (8A6) 

GO  TO  100(J 
C 

C  NCON=6,9 

C  DEFINE  PERSPECTIVE  CONSTANTS 
8      IOP(NCON)=l 

IF  (NCON-8)   7083,7084*7083 
7084   NBOUND-"=l 
7083   CONTINUE 

[)p=TFMP(  1} 

DO=TEMP(2) 

GO  TO  10o0 
C 

C  NCON=10 
C  DEFINE  TRANSLATION 

10  DX(1)=TEMP(1) 
DYU)=TEMP{2> 
DZ(1)=TEMP(3) 
IOP(NCON)=l 
GO    TO    1000 

C 

C    NC0N=11 

11  ISTAP=TEMPU) 
GO  TO  1000 

C 

C  NC0N=12 

12  IOD(12)=l 

READ  C5»556>  ( T IT ( I ) , 1=1 ,8 ) 
WRITE  <6»1008)  «TIT(I)»I=1*8) 
556   FORMAT  (8A6) 

GO  TO  10o0 
10^8  FORMAT  (1H0v12A6) 
C 

C  NC0N=13 
C  FNP>  PLOT  TAPE 

13  CALL  PLOT  (0,0*999) 
STOP 

C 


^ 


C  NC0N=14 

C  READ  LABELS  AND  THEIR  CORRESPONDING  COORDINATES 

14     I0PU4)  =  1 

ILAB=ILAB+1 

XL(IL/a)=TEMP(I) 

YL(IL/8)=TEMP{2/ 

IF    (ZiCALE)    92,96,92 
96  ZL(ILAB)=TEMP(?) 

60  TC  94 
92     ZU  J)«TEMP(3)/ZSCAL-" 
94     I  Til  LAB) =TEMP (4) 

-*EAD  (5,102)  (LAB(I»ILAB)tI*l,6) 
102    FORMAT  (6A6? 

GO  TO  1000 
C 

C  aC0N«15 

C  DEFINE  CONTOUR  INCREMENTS 
13     CBGN=TFMP<1) 

CINC«TEM°(2) 

CEND=TF'S.P(3) 

I0P(15}=1 

IOPV=C 

GO  TO  1000 
C 
C  NC0N=16 

16  C6G!^=TFMP(1) 
CINJ=TEMP<2) 
CEhD=TFMP(3> 
ICPV=0 
TdWIND  3 

GO  TO  2065 
C 
C  NC0N=17 

17  NVCON=TEMP(' 
I0P(15)=1 
I0PV=1 

READ    (5*1-5)     ("FMTt I)  ,1=1,12) 
READ    <5*-FMT)     ( VCON < I ) , 1=1 ,NVCON) 
GO    TO    1J0u 
C         NC0N=18 

18  IOprJ)  =  l 
SIZ'=TFMP(1) 

go  ro  looo 


® 


C  NC0N=19 

C  READ  FORMAT  FOR  MATRIX 

19     READ  (5»1Q5)  { FMT U ) » 1=1 »12 ) 

105    FORMAT  (12A6) 

IOP(19)=l 

GO  TO  I^oj 
C 

C  NCON=20 
C  BEGIN  PLOTTING 
20    IFt I0P{2j) .EQ.l)  GO  TO  953 

IF  (IOP(6).EQ.O)  GO  TO  6001 

CALL  PLOTS  (WORK(l) ,1024»C> 

CALL  SYMBOL{0»,0«9<»2»USER,90»»48  3 

CALL  PLOT(l.,0*»-3) 

IOP(20}=1 
953   IF  <I0P{8}.EQ.0?  I0P{9)*1 

IF  (IOP(18)  *EG»  0)  GO  TO  954 

XCEN={MAXX+MlNX>/2. 

YCEN=(MAXY+MINY)/2« 

YCEN1=YCEN 

XCFNl=XCEN 

IF    U0PU7)    .EQ.    0}    GO   TO    955 

C8GN=VCON(l) 

CEND=VCONCNVCON) 

955  IF  (IOPC2)  .EQ.  OJ  GO  TO  956 
YCEN1-CTHETA(2)*YCEN-CBGN*STHFTA(2)/ZSCALE 
YCEN=CTHETA(2)*YCEN-CFND*STHETA(2)/ZSCALE 

956  IF    UOP<3)    .EQ.    3)    GO    TO    957 
XCFN1=CTHETM3)*XCEN+CSGN*5THETA(3)/ZSCALE 
XCEN=CTHETA(3)#XCEN+CEND*5THETA£3)/ZSCALE 

957  IF    (IOP«3)     .EQ.    0)    GO    TO    958 
XC£N1=CTHETAU)*XCEN1~STHETAU>*YCEN2 
XCEN=CTHETA?4)*XCEN-STHETA(4)*YCEN 

958  XCEN=«XCENl+XCEN)/2c 
XCFN=XCFN/SCALF+XMAX 
XMAX=XMAX-XCEN+SIZE/2e 

954    IF  (IOP(l)  .EQ.  0)  GO  TO  904 

C 

C  CALL  DRAFT  TO  PLOT 

C5GN=-53139.E15 

NCURV=0 

NCOUNT*0 

NCMAX=0 

IF  (NLINES)   23*23»25 
25     READ<3)  N»CV 

IF  (CBGN.EQ.  -53139. E15)  CBGN=CV/ZSCALE~. 5 

READ  (3)  (X(I  ),Y(I) »I=1.N) 

NCURV*NCURV+1 

NCOUNT=iMCOUNT+N 


7* 


IF  (N-NCMAX)   7337*7338,7338 
7338   NCMAX=N 

CLMAX=CV 
7337   CV=CV/ZSCALE 

CALL  DRAFT 

IF  (NCURV-NLINES)   25,22,22 
23   CRr,N=Oe 
22     N=5 

CVsCBGM 

X ( 1 ) «FNN 

X(2)s-X(D+1» 

X(3)*X{2) 

X(4)=X«1) 

X(5)=X(D 

Y(1)=FM2 

Y(2)=Y<1) 

YC3)*-Y(2)+1. 

Y(4)=Y(3) 

Y(5)=Y(D 

CALL  DRAFT 

XNEW=YMAX*2, 
28    CALL  PLOT(XNEW,0. .-3) 

GO  TO  loOG 
C 

C  NCON=21 
C  CHANGE  SCALE 
21     SCALE  =(FLOAT(M)-l,)/TEMP(l) 

GO  TO  1000 

6001  WRITF  (6,60023 

6002  FORMAT  < 4lH0****PLOT  USER  IDENTIFICATION  NOT  DEFINED) 
GO  TO  920 

9C2    WRITE  (6,402) 

402  FORMAT  ( 32H0****FORMAT  NOT  DEFINED  FOR  DATA) 
GO  TO  920 

903  WRITE  (6,403) 

403  FORMAT  ( 3lHQ*#**C0NT0UR  LEVELS  NOT  DEFINED) 
GO  TO  920 

904  WRITE  (6,40-*) 

404  FORMAT  ( 12H0»***NO  DATA) 
GO  TO  920 

30     WRITE  (6,531)  NCON 

531    FORMAT  ( 27H0****NCON  NCT  LEGAL*  NCON  =  ,15) 

920    WRITE  (6,6003) 

6003  FORMAT  ( 2CH0****PLQT  TERMINATED) 
STOP 

END 
SIBFTC  SCANS    DECK 
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su**nouTiwr  sc*w 

3lu£V5IO*J    Au(  1<~"">,1"0)  ,rrC(BOo)  »y(l«5Art)  »Y(15™)  tI°T(3i3J>K!X(81t 

1  I»!Y (  e  )  »  I0n  (  2" )  txi. (  in  )  ,YL  <  1C< )  «7L  ( 1  w  >  •LA°  C  10«6 )  »0X I  2  )  »DV(  2 >  »DZ  (  2 )  » 

2  HT(10)  »CTHETA(4)  »«?TH~TM4I  »Z<  15«-"M  *TITf  8  )  »wORK  1 10241 

COMMON  MT»NTtNI»IX»lY*IDr»IDYfI?S,IT»lV,KP.v,JT»PY»RE:C»CV»X»v,lPT9 

1  INX»lNY»DLfAMtIOP»DPtDO»ILA3»XLtYL»ZL»LAR»DX»DY,*7tY"AX»HT»r.CA!_F» 

2  YMAX»CTHETAtSTHETA9FM2»FNN»NN»M»CL.DtNLINES»NCunVtZtTIT,».ORK 
I2X  =  1 

NP=0 

DL=D 

NT=M 

MT=NN 

IF  (IZX)   3,3,1 
1  !PT(ltl»«8 

IPT(1»2)=1 

IPTUt3)»2 

IPT(2»1)=7 

IPT<2t3>=3 

IPT(3.1)=6 

IPT(3»2)=5 

IPT13»3)»4 

INX(1)=~1 

INX(2>=-1 

INX«3)=U 

INX(4)=1 

INX(5)=1 

INX<6)=1 

INX(7)*w 

INX<8)=-1 

INY<1)=^ 

INY(2)=1 

INY(3)=1 

INY(4>*1 

INY(5»=G 

INY(6)=-1 

INY(7)=-1 

INY(8)=-1 

1ZX=C 
3  XT=MT 

DO  58  Jsl»800 
58  REC(J)=0. 

ISS=0 

MT1=MT-1 

DO  110  1=1, MT1 

IF  <AM(I,1)-CV>  559110,113 
55  IF  (AM(I+1»1}-CV)  110,57,57 
57  IX«I+1 


IY=1 

IDX=-1 

IDY=0 

CALL  TRACE 
lie    CONTINUE 

NT1=NT-1 

DO  20  I=1»NT1 

IF  <AM(MT»I)-CV>  15,20,20 
15  IF  (AM<MT,I+1)-CV>  20,17,17 
17  IX=MT 

IY=I+1 

IDX  =  0 

IDY*-1 

CALL  TRACE 
20  CONTINUE 
22  DO  30  1=1, MT1 

MT2=MT->1-I 

IF  (AMCMT2fNT)-CV)  25,30,30 
25  IF  <AM(MT2-1,NT)-CV>   30*27,27 
27  IX=MT2-1 

IY=NT 

IDX  =  1 

IDY»0 

CALL  TRACE 
30  CONTINUE 

DO  40  I=1»NT1 

NT2=NT+1-I 

IF  <AM<1,NT2)-CV)  35,40,40 
35  IF  <AM(1»NT2-1)-CV)   40,37,37 
37  IX=1 

IY=NT2-1 

IDX»0 

IDY=1 

CALL  TRACE 
40  CONTINUE 

ISS=1 

NT1=NT-1 

MT1=MT~1 

DO  10  J=2»NT1 

DO  10  I=1»MT1 

IF  (AM(I»J)-CV)   5,10,10 
5  IF  JAM(I+1,J)-CV)   10,7,7 
7  C0M=100*U  +  1)+J 

IF  (NP)   12,11,12 
12  DO  9  10=1, NP 

IF  (REC( IDJ-COM)   9,10,9 
9  CONTINUE 
11  IX=I+1 

IY=J 

IDX*-1 

IDY=0 

CALL  TRACE 
10  CONTINUE 

RETURN 

END 


\ 


£*£» 


SIBFTC  TPACFS   OFCK 

SUBROUTINE  TRACE 

DIMENSION  AM ( 100 , 100 ) ,REC { BOO )  ,X  < 1500 > »Y ( 1500 )  ,  ?PT ( 3 ,3 ) , iNX  t  8 ) » 

1  1NY(8)  tIOP(20)  tXL(10)  tYL(10)  tZL(lC)  »LA3(10»6)  ,DX  {  2  )  }>DY(  2  )  „DZ«2  )  , 

2  HT(10)»CTHETA(4) ,STHETA{4 J ,Z( 15U0 ) *T IT* 8 ) ,WORK U024) 

COMMON   MT»NT,NI»lX,lY,  I  DX»  IDY,  ISS,  IT  ,  I  V»NP»N,JT  j>PY*REC,CV»X»  ,'»IPT» 

1  INX,INY,DL»AM,IOP,OP,CO» I  LAB, XL »YL,ZL»LAB*DX »DY»DZ tXMAX •HT .SCALE , 

2  YMAX»CTHETA,STHETA»FM2,FI\!N»NN,M»CL,D»NLINES,NCURV,Z, TIT,  WORK 
COMMON  MAXX»MAXY»MINX»MINY»SIZE 

REAL  MAXX»MAXY.MINX,MINY 
PY=0.0 
501  JT»0 
N*0 

IX0*IX 
IYO=IY 
ISX=IDX+2 
ISY=IDY+2 
IS=IPT(ISX,I5Y) 
JTB»C 
ISOMS 
IF(IS0-8)   18»18»17 

17  IS0=IS0-8 

18  IT»0 

5      CALL  CALC 

NZ=N 

N»NZ 

IF  (IT+JT-1)  49,49,47 
47  XS»X(N~1) 

YS=Y(N»1) 

X(N-l)aX(N) 

Y(N-D=Y(N) 

X(N)=XS 

Y(N)=YS 
49  IS=IS+1 

JT=IT 
9   IF  US-9)   8.797 

7  IS*IS-8 

8  IDXslNX(IS) 
I0Y~INY( IS) 
IX2=IX+I0X 
IY2=IY+IDY 

JT^sJTP+1 

IF   (JTB-1500)  51,51,308 
308    WRITE  (6,103)  CV 

103    FORMAT  <21H0  A  CONTOUR  AT  LEVEL  »F6*2»16H  WAS  TERMINATED..) 
RETURN 
51  IF  (ISS)  10,10,20 

20  IF  (IX-IXG)   12,21,12 

21  IF  (IY-IYO)   12,22,12 

22  IF  (IS-ISC)   12,23,12 

23  CALL  CALC 
GO  TO  73 


7fg 


10  IF(IX2)       13.50,13 

13    IF    (IX2-MT)    19*19,50 
19    IF(IY2)    11,50,11 

11  IF    (IY2-NT)       12,12,50 

12       IF    (CV-AM(IX2,IY2))    206,206,5 
206         IF    (IDX**2+IDY**2-1)    213,6,213 

213  DCP*{AMC IX*IY)+AM(IX2, IY)+AM( IX, I Y2 }+AM( 1X2, I Y2 ) )/4.0 
IF    (DCP-CV)    5,217,217 

217         IF    (INXCIS-in    214,215,214 

214  IX»IX+1DX 
IDX»~IDX 
PY«2.0 
CALL   CALC 
IXMX+IDX 
GO    TO   6 

215  IY=IY+IDY 
IDY=-IDV 
PY=2.0 
CALL  CALC 
IY=IY+IDY 
6   IF  {AM(IX-1 »IY)-CV)    306,16,16 
306  NP-NP+1 

RECCNPl^lQO^lX+lY 
16  IS=IS+5 
IX=IX2 
IY=IY2 
GO  TO  9 
50   XT=MT 

IF  (AM(lX-liIY)-CV)  307,73,73 
307   NP=NP+1 

REC(NP)=100*IX+IY 
73     NLINES*NLJNES+1 
DO  135  K«lfN 
X(K>=X(K}-FNN 
Y(K)=Y(K)~FM2 
135    CONTINUE 
C 

C  STORE  CURVE  ON  TAPE  3  $A4> 
IF  (N)  2072,2072,2071 

2071  WRITE  (3)  NtCL 

071   WRITE  (3)  <X< I) »Y(I),I=1,N) 
WRITE  (6,104)  CV 

104  FORMAT  < 16HOCONTOUR  LEVEL  =,F6*2) 
WRITE  (6*105)  (X(I) ,Y(  15 r i  =1 #N> 

105  FORMAT  (2X96F10.3) 

IF  (IOPC18)  .£0*  0)  GO  TO  2072 
DO  72  I- UN 

IF  (XU>  .GT.MAXXI  MAXX  =  X(I? 
F  (Yd)  sGT.MAXY)  MAXY«Y(I) 
IF  (X(I)«LT»MINX)  KINX-X(I) 
72     IF    (Y(I)  .LT.MINYJ  MINY«YCI) 

2072  N«-l 
RETUPN 
END 


SIBFTC  CALCS    DFCK 

SUBROUTINE  CALC 

DIMENSION  AM < 100,100) ,REC(80U J  ,X(15O0)  »Y(1500J  «.  I  PT  (  3»3  )  »  INX(  8  )  » 

1  INY(8)»IOP<20>  »XL(10) ,YL(10J ,ZL<10S  *LAB{ 10 s6 ) ,DXC 2 ) ,DYC 2 ) ,D2 < 2 ) , 

2  HT(10).CTHETA{4),STHETA(4),Z{ 15«0 » pTIT I  8 ) t WORK (102*) 

COMMON  MT»NT»NI*lX»lY»IDX»IDYtlSS»IT»lV»NPtN»JT»PY»REC»CV»XfYtlPT» 

1  INX  »  INY  »DL  »  AM  ♦  I  OP  »DP  9  DO.  I  LAB,  XL  s.YL,2L*LA3»DX.DY»0Z*XMAX*HT»  SCALE  9 

2  YMAX»CTHETAtSTHETA»FM2>FNNtNN».V»CL*DtNLINES»NCU^V»2.TIT»W0RK 
IT=0 

N=N+1 

IF  (IDX**2+lDY**2~li   20.1»20 

1  IF  (IDX)   10»2.10 

2  X(N)=IX 
Z=IY 

IY2»IY+IDY 
DY=IDY 

41     Y(N)  =  ((AM(IX*IY)-CV)/(AM< | x» I Y)-AM { I  X  * IY2 J ))*DY+2 

RETURN 
10     Y(N)=IY 

W»IX 

DX«IDX 

IX2=IX+IDX 
44    X(N)  =  ({AM(IX»IY)-CV)/(AM{ I  X ♦ I Y) -AM( 1X2  * i Y ) ) )*DX+W 

RETURN 
20  IX2»IX+IDX 

IY2»IY+IDY 

W*IX 

Z  =  IY 

DX=IDX 

DY=IDY 

DCP*<AM( IX»IY)+AM(!X2* IYJ+AMC IX»IY2 )+AM( I X2 • I Y2 ) ) /4«G 

IF  (PY-2-0)   24»21»24 

24  IF  (DCP-CV5  21*21.25 
21   AL=AMUX»IY)-DCP 

23   V=.5*(AL+DCP-CV)/AL 

27  X(N)»V*DX+W 
Y(N)»V*DY+Z 
PY*0.u 
RETURN 

25  IT*1 
AL=AM(TX2»TY?>-DCP 

33   V=.5*(AL-*-DCP-CV)/AL 

28  X(N)*-V*DX+W+DX 
Y<N)=-V*DY+Z+DY 
RETURN 

END 


SIBFTC  DRAFTS   DFCK 

SUBROUTINE  DRAFT 

DIMENSION  AM ( 100.100)  »REC ( 800)  .X< 1500) »Y < 1500 ). IPT( 3.3) . INXC8) » 

1  INY(8).IOP(20).XL(10)  ,YU10>  »ZL(10) »LA6( 1C.6) »DX(2 ) »DY(2) ,DZ(2) . 

2  HT(IO) .CTHETAC4) »STHETA(4)  »Z(  15o0)  »TIT(8>  .WORM  1024) 
DIMENSION  °APTK2) 

COMMON  MT»NT»NI»IX.IY»IDX»IDY»ISS»IT»JV»NP»N.JT»PY.REC.CV.X.Y.IPT. 

1  INX»INY»DL»AM»I0P.DP»D0.ILAB,XL»YL»ZL»LA8,DX»DY»DZ.XMAX»HT.SCALE. 

2  YMAX»CTHETA.STHETA,FM2tFNNtNN»M»CL»DtNLINEStNCURVtZ«TIT»W0RK 
COMMON  MAXX»MAXYtMINX»KINY»SIZE 

COMMON  ZSCALE 
REAL  MAXX.MAXY.MINX.MINY 
303  KTYPE=3 

DO  80    I=2»14 
IF   (XOP(D)    80.80,81 
81  GO  TO  {101tl01»102»l03.80»3C»80»8»8»10»80.80»80»l4)»I 
C 
C  ROTATION  ABOUT  X 

101  01*1. 
D2«0, 
D3=0. 
D4=0. 

D5=CTHETA(2) 
D8=STHETA(2) 
D7=C. 
D6«-D8 
D9=D5 

GO  TO  106 
C 
C  ROTATION  ABOUT  »' 

102  D1*CTHFTA(3) 
D2=0. 

D3=STHETA(3) 
D4=0. 
D5*l. 
D6=0„ 
D7=-D3 
D8=0» 

D9  =  D1 
GO  TO  106 
C  ROTATION   ABOUT   Z 
103  Dl*CTHETA(4) 
D4=STHFTA(45 
D3=0. 
D2=-D4 
D5=D1 
D6=0» 
D7=0S 
D8=0» 
D9=l. 


C  ROTATE 

106    CVD3=Cv*D3 

CVD6=CV*D6 

CVD9=CV*D9 

DO  20  J=1»N 

XX=D1*X( J)+D2*Y( J)+CVD3 

YY=D4*X( J)+D5*Y(J)+CVD6 

Z ( J ) =D7*X ( J ) +D8#Y  <  J  J  +CVD9 

X(J)=XX 
20     Y(J)*YY 

GO  TO  80 
C 

C  PRINT  HEADING  INFORMATION 
8      TITC0Ds,9 

IF  (1-8)  401,401*301 
401    NN1=2 

WRITE  (6.1010)  DP, DO 
1010  FORMAT  (24H0     DISTANCE  TO  PLANE  =jF9*2»22H,  DISTANCE  TO  OBJECT  - 
1»F9.2) 

GO  TO  86 
301   NN1=1 
86    IF  (IOP(12))  565,566,565 

565  CALL  SYMBOL  (TITC0D»3«»»14»TlTt90.»48 ) 

566  IOP(8)«0 
IOP(9)=0 
GO  TO  80 

C   PLOT   LABFL   ON   GRAPH 
14   DO  935   J=1,ILA8 

GO  TO  (751,752), NN1 
752  P»(DO-ZL«J) )/DP*SCALF 

F*(DO/DP)*SCALE 

C1=(XL(J)-F  INJ/E+XMAX 

C2*  ( YL  ( J  ) -Fi-12  )  /E+YMAX 

GO  TO  934 
751    C1=XL(J)/SCALE*XMAX 

C2=YL(  J)  /SCALE-s-YMAX 

F=SCALE 

934  HT(J)=HT(J)/F 

935  IOP(14)=0 
893    GO  TO  80 

C  TRANSLATE   DATA 
10   DO  602  J=1,N 

X(J)=X(J)+DX(1) 

Y(J)=Y(J)+DY(1) 
602  Z(J)=Z(J)+DZ(1) 


Ba 


C  PLOT  PLANE  OR  PERSPECTIVE  VIEW  OF  SURFACE 
80  CONTINUE 

I OPC= IOP ( 2 ) + I OP ( 3 ) + I OP { 4 ) 
405  DO  300   J=1,N 

IF  (IOPC)   894,894.892 
894    Z(J)=CV 
892    GO  TO  {402*403) »NN1 

402  Ci=X< Jj/SCALF+XMAX 
C2=Y(J)/SCALE+YMAX 
GO  TO  406 

403  E=(DO-Z(J> )/DP*SCALE 
C1  =  X(JJ./E+XMAX 
C2=Y( J)/E+YMAX 

406    IF  (I0P(5))  933*933*930 

930    CALL  PR0JfCl»C2) 

933    IF  (Cl.GT.SIZE.OR.Cl.LT.-l.)  GO  TO  936 

CALL  PLOT  (C2,C1»KTYPE) 

IF{KTYPE.E0.2)  GO  TO  304 

CVV»CV*ZSCALE 

CALL  NUMBER  { C2 »C1 , «07,CVV990. ,2 > 

CALL  PLOT  {C2.C1.3) 

CB1«C1 

CB2=C2 

304  KTYPE=2 
GO  TO  3u0 

936    KTYPE=3 

300    CONTINUE 

IFICB2oECUC2.AND.CBl.EQ.Cl)    GO    TO    305 
CALL    NUMBER    ( C2 »C1 • .07»CVV,90. ,2 ) 

305  CONTINUE 
RETURN 
END 


SIBFTC  PROJS    DECK 

SUBROUTI NE  PRO- i  X «  Y  5 

X«l» 

Y*l. 

RETURN 
END 


SIBMAP 


INPUT 

ENTRY 

p;:e 


•  UNU8 
UNI  TO 8  FILE 
END 


«UN08. 
UNIT08 
»  A  »  READY  *  I NPUT  i  3L  ;>22  *Hl GH ,BCD 


PROGRAM  CONTUR  "  SAM! 


m 


595.      5,         720. 

H   DEGREE  TREND  SURFACE«APR !{_  1967 
28.0 


23. 

672.534 
548.723 
664.519 

665.634 
647.929 
663*969 
659*002 
646.836 
663.090 
63?. 666 
645.479 
661.912 
646.6r>5 
643.387 
660.464 
640.994 
642.038 
658.777 
635.713 
640.111 
656,878 
630.838 
637.982 
654.797 
626.394 
635.729 
652.562 
622.409 
633.379 
650.199 
618*907 
630,958 
647,738 
615.914 
628*493 
645.203 
613,453 
626.008 
642.623 
611.551 
623.529 
640.022 
610.229 
621.081 
637.426 
609.51'! 
old. 687 
634.861 
609.424 


11.0 
670.131 
648.240 
666.728 
664.417 
647.730 
665.949 
658.373 
646.393 
664.8  73 
653.525 
645.759 
663.530 
648.404 
644*356 
661.952 
643,535 
642.714 
660.165 
638.947 
640.860 
658. 2C1 
634.666 
638.823 
656.088 
630.719 
636.629 
653.853 
627.132 
634. 3C 5 
651.524 
62  3.930 
631.878 
649.130 
621,139 
629.375 
646.696 
618.783 
626.319 
644.250 
616.888 
624.238 
641.817 
615.476 
621.655 
639.423 
614.371 
619.096 
637.093 
614.196 


666.802 

648.719 
667.694 
662.342 
648.412 
666.648 
6  57.560 
647.75* 
665.345 
653*065 
646.773 
663.817 
643.745 
645.500 
662.093 
644.568 
643.961 
660.203 
640. 581 
642.187 
653.177 
636,811 
640*205 
656.04? 
633.265 
638*04  2 
653*827 
630.029 
6  35*724 
651.560 
627.068 
633.280 
649.268 
6  24.429 
630.735 
546.979 
622.135 
623.114 
644.719 
620.212 
625.444 
642.514 
61P.683 
622.748 
640.390 
617.572 
620-052 
638.373 
616.902 


663*009 
650.114 
666. 953 
659.275 
649-3  931 
665.50? 
655.538 
649.380 
663.954 
651.324 
648.490 
662*224 
648*163 
647.290 
660.348 
644*583 
645.809 
658*355 
641.111 
644.076 
656.275 
637.775 
642.118 
654,136 
634.600 
639.963 
651 #966 
631*614 
637*639 
649. 3  794 
628.842 
635.171 
647.647 
626.309 
632.587 
646. 552 
624.040 
629.911 
643.537 
622,061 
627. 170 
641.627 
620.393 
624.389 
639.848 
619 .066 
621.591 
638.226 
618*097 


659.149 

652.309 
:i  6  3*575 
656.221 

6  5  2  : 
661.904 
653*214 
651.666 
660* 
650,157 
650*309 
658*141 
647.079 
649.634 
656»iil 
644.008 
648-.  170 
654*022 
540.972 
6 -A  6*  44  5 
651.903 
637.997 
644*488 
649.782 
635.110 
642*325 
647.689 
632,338 
639*936 
645*.  651 
629.706 
637.495 
643.697 
»241 
634.380 
641.853 
624.967 
632. 166 
640.146 
622*910 
629.360 
638.603 
621.092 
626.546 
637.250 
619.538 
623.688 
636.111 
618..27? 


655.354 
655.123 
657,168 

653*316 
654*973 

655.151 
650*933 

6  54  «  4»'t  T: 
653.049 
6  ft  8  .  £  3  5 
653=570 
650.893 
6  4  5  .  8  5  0 
652*376 
648*713 
643.205 
650.894 
646 * 540 
640.529 
649.151 
644s 402 
637.850 
647.175 
642.3  30 
6  3  5.193 
644.995 
640.350 
632.586 
6^2 o638 
638.492 
630.055 
540.131 
636.784 
627.624 
637.499 
635.252 
625.320 
6  34.770 
633*924 
623*167 
631.969 
632.827 
621.139 
629.121 
631.983 
619.411 
626.251 
631.425 
617.856 


616.372 
632. 350 
609*985 
614*139 
629.920 
611.219 
612.071 
627.  50? 
613, 146 
510.131 
625,391 
615.79C 
608.361 
623,340 
619,169 
606.782 
621*461 
623.305 
605,416 
629.777 
628.218 
604.284 
618.309 
633.923 
6W3.406 
617.079 
640.453 
602.803 
616.107 
647.012 
602.493 
615,415 
656.024 
602.497 
615.02! 
665.106 
602.832 
6X4.947 
675.077 
6o3.51S 
615.210 
6B5.953 
604.573 
615.831 
697.751 
606.013 
616.826 
710.487 
607*856 
618,215 


616*183 
634. «5 2 
614,375 
614.141 
632.725 
615.128 
611.792 
630.735 
616*479 
609*560 
62  8.907 
618*448 
607,466 
627,262 
621*057 
605*533 
625*325 
624.326 
603*782 
624,616 
628*275 
602.233 
623«659 
632*923 
600.9C8 
622*974 
638*291 
599*827 
622.563 
644.397 
599. 0C8 
622*506 
651.?60 
598*473 
622.763 
658.897 
593*239 
623,375 
667,326 
593.325 
624,360 
676,565 
59 d. 750 
625,737 
686.630 
599.530 
627.526 
697,537 
600. 683 
629.743 


617.379 
636.487 
616*696 
614.753 
634*757 
616*976 
612.198 
633.207 
617.765 
609.736 
631,860 
619*083 
607*3b9 
630*740 
620*952 
605*180 

62°. 870 

623*393 
603.130 
629*2  72 
626*426 
601*260 
628*968 
630.070 
599.591 
628*979 
6 34 c 345 
598.143 
629*327 
639,270 
596*937 
630«033 
644.863 
595.990 
631*117 
651,144 
595*323 
632*59B 
658*128 
594.954 
6  34  «,49  6 
665*835 
594*901 
636*831 
674.280 
595*182 
639*620 
683*480 
595.815 
642.882 


618.802 
636.786 
617.511 
616.045 
635*552 
617*330 
613.343 
634*548 
617. 57P 
610.719 
633*799 
618*274 
608,196 
633.328 
619.441 
605*796 
6339157 
621*099 
603.540 
633.309 
623*269 
601,450 
633.307 
625.970 
599*547 
634*671 
629.222 
597. S50 
6  35.924 
633.045 
596.330 
637»586 
637.456 
595,156 
639.677 
642.474 
594.198 
642«,217 
648*117 
593,524 
645.227 
654*403 
592*152 
648.72* 
661,347 
593.101 
652*729 
668*963 
593.388 
657*258 


620.832 
635.214 
617.31? 
618*001 
634.581 
616.694 
615*218 
634.238 
616.427 
612.507 
634*208 
616*536 
609*890 
634*515 
617*044 
607*389 
635,181 
617.971 
605*027 
636*2  30 
619*337 
602*823 
637*684 
621*163 
600.800 
639*564 
623*469 
598.978 
641*892 
626*273 
597*376 
644*689 
629*594 
596*015 
647.975 
633*451 
5949914 
651.770 
637*861 
594*091 
656,094 
642.842 
593*566 
660*96? 
648*412 
593*355 
666.406 
654*586 
593.476 
672*431 


623.383 
631.173 
616.547 
620.542 
63U253 
615.507 
617*751 
631.689 
614.75? 
615*032 
632.506 
614.321 
612*409 
633.727 
614*218 
609.904 
635.375 
614*471 
607.539 
637*473 
615.100 
605,336 
640*043 
616,125 
603,314 
643*107 
617,565 
601*496 
646*687 
619*44! 
599.900 
650*804 
621.770 
598*548 
655*478 
624.572 
597.457 
660.729 
627.864 
596.648 
666.578 
631*665 
596*138 
673.043 
635*990 
593*946 
680.144 
640.358 
596*090 
687.398 


626.307 

615.571 
623.524 

614.143 
620.802 

612.951 
618.162 

612*015 
615.628 

611*358 
613.221 

611.000 

610*965 

610*962 

608*879 

611*265 
606.987 

611.929 
605.30? 

612*972 

603^861 

614.414 
602.669 

616.273 
601*749 

618.568 
601*122 

621*316 
600.804 

624.535 
6OO08I6 

628*241 
601*174 


xvo.sia 

sp?.rs* 

rri.^s 

a*r,na 

.TC 

IEI, 

SAMPLE    OUTPUT    FROM     PROGRAM     CONTUR 
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